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INTRODUCTION 

Th«  work  performed  and  the  results  obtained  under  Contract  • 
FH-11-7904  to  the  Federal  Highway  Administration,  Structures 
and  Applied  Mechanics  Division,  for  a  study  on  "Prediction  of 
Loadings  on  Highway  Bridges--Phase  II",  is  described  in  four 
separate  reports,  FHWA-RD-73-42,  43,  44  and  45. 

A  primary  objective  of  this  study  was  to  extend  the  traffic 
simulator  computer  program  developed  by  the  Kelly  Scientific 
Corporation  to  a  useable  engineering  tool  capable  of  generating 
bridge  loads.   It  was  further  the  purpose  of  the  study  to  develop 
a  finite  element  stress  analysis  program,  including  dynamic 
effects  due  to  the  live  load,  which  directly  interfaced  with 
the  load  generator,  produce  several  stress  histograms  for 
various  bridge  and  traffic  configurations,  and  to  develop  and 
implement  an  alternative  method  based  upon  analytic  methods 
rather  than  traffic  simulation. 

Report  FHWA-RD-73-42  describes  the  work  performed  in  the 
study.   This  includes  the  background  which  led  to  the  work 
performed,  a  description  of  the  work  performed  and  problems 
encountered  in  revising  BRIGLD1.   The  results  of  sensitivity 
testing  of  BRIGLDl,  development  of  the  stress  program,  genera- 
ted histograms,  description  of  the  analytic  methods  investigation, 
and  the  conclusions  and  recommendations  are  also  included. 

Report  FHWA-RD-73- 43  is  the  Users  Manual  which  provides 
utilization  instructions,  data  preparation  instructions,  output 
and  variables  definitions  and  a  description  of  the  BRIGLDl 
computer  program.   The  description  includes  a  narrative  description 
section,  flow  charts  and  program  listings  on  the  main  line 
program  and  each  subroutine.   The. use  of  this  program  for  stress 
range  prediction  is  illustrated  in  Figure  1. 
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Figure  1.   Traffic  Simulation  Based  System 
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Report  FHWA-RD-73-44   is  the  Users  Manual  for  the  dynamic 
stress  analysis  computer  program,  BRGSTRS.   This  report 
contains  the  same  descriptive  type  matter  as  indicated  above 
for  BRIGLD1. 

Report  FHWA-RD-73-45   contains  the  Users  Manuals  for  two 
computer  programs  which  operate  as  a  system  with  BRGSTRS. 
The  first  is  the  Synthetic  Load  Generator,  SYNGEN,  which 
generates  single  axle  loads  for  the  dynamic  stress  analysis 
program,  which  in  turn  generates  a  stress  signature  curve,  trace, 
for  each  defined  axle  load.   The  second  is  the  histogram 
computer  program,  HISGEN,  which  generates  long-term  stress 
range  histograms  from  the  synthetic  single  axle  stress  trace 
data  generated  by  the  dynamic  stress  analysis  program,  as  shown 
in  Figure  2.   This  is  accomplished  by  first  forming  composite 
truck  stress  traces  for  a  given  truck  population  from  the 
single  axle  data.   Then,  forming  composite  truck  platoon  stress 
traces  for  a  given  platoon  population.   Long-term  effects  are 
estimated  from  traffic  density  estimates  and  the  estimated 
incidence  of  each  platoon  configuration.   The  information  contained 
in  this  report,  for  both  SYNGEN  and  HISGEN,  is  of  the  same  form 
as  described  above  for  BRIGLD1. 
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Figure  2.   Axle  Signature  Based  System 


BACKGROUND 

Reference  (1)  describes  the  work  performed  under  contract- 
FH-ll-7314  in  developing  a  traffic  simulator  computer  program 
(BRIGLD1)  for  the  purpose  of  Forecasting  of  Heavy  Loading 
Patterns  on  Highway  Bridges.   The  purpose  of  this  work  was 
toward  developing  better  techniques  for  predicting  long-term 
fatigue  of  bridges,  specifically  the  main  beams.   This  work 
was  completed  in  June  1970. 

The  initial  form  of  this  program  was  developed  on  a  CDC6600  in 
Run  FORTRAN.   Subsequent  to  delivery  of  the  program  to  the 
Federal  Highway  Administration,  the  program  was  modified  to 
be  compatible  to  IBM  (360/65)  FORTRAN  IV. 

In  February  of  1971,  the  Federal  Highway  Administration  solic- 
ited research  and  development  study  proposals  for  Phase  II  of 
the  work.   The  problem,  which  the  research  and  development 
study  was  to  attack,  was  defined  as  follows: 

"The  bridge  loading  simulation  model  in  its  original  state 
was  able  to  forecast  certain  combinations  of  loadings  on 
bridges.   Parameters  accounted  for  in  the  model  were  total 
traffic,  percent  of  truck  traffic,  separation  of  trucks  by 
types ,  separation  of  types  into  several  axle  groupings  and 
weight  distributions,  geometric  highway  features,  vehicle 
speed,  passing  maneuvers,  traffic  direction,  and  number  of 
lanes.   Of  interest  was  the  summation  of  certain  critical 
loading  combinations  on  various  lengths  of  bridges. 

"The  model  was  not  thoroughly  tested  and  stopped  short 
of  converting  loading  combinations  into  stresses.   The 
model,  in  its  original  configuration,  was  only  an  inter- 
mediate research  tool  which  was  not  directly  useable  by 
bridge  engineers. 


"Accompanying  the  development  of  the  model  was  a  limited 
amount  of  analytical  development  of  bridge  loading  fore- 
casting.  The  analytical  work  was  intended  to  supplement 
results  from  the  model. 

In  June  1972 ,  this  firm  was  awarded  the  Phase  II  study  by  the 
Federal  Highway  Administration. 

The  objectives  of  the  study  were  to: 

1.  Perform  a  sensitivity  analysis  on  the  simulation 
model  (BRIGLD1)  to  evaluate  the  effects  of  changes  in 
all  the  important  parameters. 

2 .  Extend  the  scope  of  the  model  to  include  the  calcu- 
lation of  stress  in  representative  bridge  members. 

3.  Prepare  charts  or  tables  of  stress  histograms  based 
on  several  bridge  and  traffic  configurations  which  are 
suitable  for  use  as  design  guides. 

4.  Extend  the  analytical  approach  to  load  estimation, 
and  use  the  resulting  analysis  for  comparison  and  veri- 
fication of  the  simulation  results. 

In  order  to  satisfy  these  objectives,  the  work  was  defined  in 
terms  of  four  major  tasks,  i.e., 

1.  Sensitivity  Analysis 

2.  Extend  Model  to  Stress  Calculations 

3.  Prepare  Stress  Histograms/Tables 

4 .  Extend  the  analytical  Approach 

The  specifications  of  the  work  required  under  each  of  these 
tasks  are  presented  in  the  following  subsections . 


Task  I  -  Sensitivity  Analysis 

Vary  all  the  input  data  over  the  allowed  range  of  each  element 
of  data,  from  extrema  to  extrema,  and  all  Input/Output  control 
data  and  option  elements  for  the  purpose  of  testing  and  validat- 
ing BRIGLD1. 

Establish  guidelines  for  user  selection  of  parameters  and  verify 
the  remaining  part  of  the  functional  performance  of  the  simula- 
tor and  accomplish  accuracy  testing  at  the  same  time.   This 
would  be  accomplished  through  the  use  of  a  selected  set  of  test 
cases.   Each  case  would  be  run  with  a  particular  parameter  vary- 
ing over  its  allowable  interval;  wherever  possible,  multiple 
parameters  would  be  varied  during  a  given  run,  as  long  as  their 
effect  was  independent.   Checks  of  the  results  of  each  calcu- 
lation were  to  be  performed  to  verify  the  accuracy  of  the  cor- 
responding programmed  calculation. 

Where  inaccuracies  were  detected,  reasons  would  be  investigated 
and  solutions  sought  and  incorporated  into  BRIGLD1  in 
consultation  with  the  contract  manager. 

Once  valid  functioning  and  accuracy  of  the  simulator  was  achieved, 
a  parametric  variation  on  the  significant  parameters  utilized 
by  the  simulator  would  be  performed.  .The  choice  of  parameter 
variations  would  be  guided  by  an  analysis  of  the  equations  and 
of  the  dependent  and  independent  relationships  of  the  simulator 
in  order  to  accomplish  the  desired  results  with  a  minimum  num- 
ber of  computer  runs . 

Task  II  -  Extend  Model  to  Stress  Calculations 

The  original  simulator  program  provided  time-dependent  traffic 
loadings  for  bridges.   It  accumulated  the  vehicle  weights  that 
pass  off  the  roadway  and  onto  the  bridge  within  specified  time 
intervals.   The  contractor  was  to  transform  this  data  into  a 
suitable  form  for  distribution  to  the  bridge  deck.   At  each 


time  increment  (determined  by  the  user)  the  spatial  distribu- 
tion and  associated  magnitudes  would  be  sampled  and  the  live- 
load  '  calculated  and  appropriately  distributed.   From  the  dis- 
tributed loads  the  corresponding  bridge  stresses  would  be 
calculated  at  critical  points  to  be  specified  by  the  user. 

In  calculating  the  stresses,  dynamic  effects  were  also  to  be 
considered.   The  added  dynamic  load  would  be  superimposed  on 
the  bridge  span.   The  approach  utilized  would  be  similar  to 
that  used  in  the  computer  program  from  the  Illinois  highway 
research  program.   The  actual  vehicle  axle  loadings  would  be 
distributed  to  the  lateral  and  longitudinal  beams,  and  then  to 
the  girders  via  a  2-dimensional  mesh.   Also  to  be  included 
would  be  a  user  option  that  would  allow  the  retention  of  the 
load  data,  as  a  function  of  time  and  location,  on  magnetic 
tape  for  future  use.   Continuous  as  well  as  simple  spans  would 
be  considered,  along  with  reinforced  or  prestressed  concrete 
or  steel  beam  bridges.   The  stress  routine  would  be  modular  in 
that  other  types  of  bridges,  such  as  box  beam  or  orthotropic 
plate,  might  be  added  by  the  user  if  so  desired.   The  materials 
would  be  assumed  to  be  fatiguing  metals. 

The  stress  calculation  portions  would  directly  interface  with 
the  load  simulator  at  each  increment  of  time  during  the  simu- 
lation. 

The  finite  element  method  would  be  utilized  to  perform  the 
stress  calculation.   The  assembly  of  the  stiffness  matrix  of 
the  bridge  member  with  the  nodal  forcing  function  would  enable 
the  calculation  of  nodal  displacements. 

Subsequently,  the  strain  matrix  was  to  be  obtained  and  finally 
the  stresses  would  be  computed.  The  various  existing  programs 
or  subprograms  would  be  considered  and  adapted  with  some  modi- 
fication for  tailoring  to  the  bridge  structural  members. 


The  simulation  program  would  be  verified  through  reasonable 
agreement  between  comparisons  of  the  strain  or  stress  output 
of  the  simulator  with  values  obtained  from  available  observed 
data  (data  was  to  be  obtained  through  consultation  with  FHWA) . 

From  the  developed  simulator,  the  pattern  of  loadings  on  the 
bridge  member  due  to  a  loading  would  be  utilized  in  the  stress 
histogram  would  be  determined  in  this  portion  of  the  work  and 
output  in  tabular  and  graphical  hard-copy  forms.   The  precise 
format  of  these  outputs  would  be  submitted  to  the  FHWA  for  re- 
view and  approval  prior  to  the  performance  of  Task  III. 

Task  III  -  Prepare  Stress  Histograms/Tables 

The  stress  portion  of  the  program  would  output  the  stress  his- 
tory data  as  stress-histograms  in  both  tabular  and  graphical 
forms.   A  user  option  was  also  to  be  included  to  allow  output 
of  the  ordered  time-dependent  stress  history  generated  by  the 
simulator. 

The  actual  number  and  types  of  configurations  would  be  sub- 
mitted to  the  FHWA  for  review  and  approval.   The  approval 
cycle  was  to  be  the  normal  30  days  from  time  of  submission. 
The  stress  histograms  for  approximately  ten  bridges  and  traf- 
fic configurations  would  be  prepared.   It  was  envisioned  that 
this  would  consist  of  six  different  bridge  configurations  of 
which  four  should  have  two  sets  of  traffic  configurations. 
The  configurations  would  allow  for  a  comparison  of  types, 
lengths,  widths,  and  loads. 

Upon  acceptance  by  the  FHWA  of  the  prescribed  bridge  and  traf- 
fic configurations  to  be  used  as  the  basis  for  developing  the 
tabular  and  graphical  stress  histograms,  the  contractor  would 
generate  the  necessary  input  data,  verify  and  validate  it,  and 
operate  the  load  simulator-stres^s  calculation  system  on  the 
data.   The  output  from  this  would  be  submitted  in  draft  form 
to  the  FHWA  for  review  and  approval  for  inclusion  in  the  final 


report.   This  would  be  done  substantially  before  the  submission 
of  the  draft  version  of  the  final  report  in  order  to  allow 
adequate  time  for  any  reruns. 

Task  IV  -  Extend  the  Analytical  Approach 

The  extension  of  the  loading  prediction  capability  through 
"analytical  methodology"  would  be  accomplished  through  the 
development  of  a  probabilistic-based  load  prediction  technique. 
This  technique  would  be  implemented  as  a  computer  program  to 
provide  load  histograms  for  test  cases.   These  test  case  load 
histograms  would  be  compared  to  the  load  histograms  produced  by 
the  BRIGLD1  Simulator  for  the  same  test  cases. 

The  probabilistic  techniques  would  utilize  such  significant  traffic 
behavioral  characteristics  as  flow  rates,  propagation/density,  flow 
patterns,  lane  changing,  vehicle  types,  platoon  lengths,  platoon 
composition,  platoon  size,  load  propagation,  and  load  distribution. 

Deficiencies  in  the  BRIGLDl  Simulator 

As  work  progressed  in  familiarization  and  investigation  of  the 
original  BRIGLDl  simulator,  it  became  evident  that  the  program 
was  almost  totally  unsuited  for  the  objectives  of  the  study. 
Further,  coding  errors  and  poor  representations  were  discovered 
in  the  simulation  model.   These  factors  caused  the  planned  attack, 
to  satisfy  the  study  objectives,  to  be  completely  revised.   The 
primary  and  priority  effort  became  one  of  salvaging  the  traffic 
simulator  for  use  as  a  bridge  load  generator.   All  other  efforts 
were  relegated  to  subordinate  positions  until  this  problem  was 
solved  adequately.   Hence,  the  sensitivity  analysis  actually 
became  a  subtask  in  the  major  task  of  salvaging  the  simulator. 

Parallel  efforts  were  carried  on  for  the  stress  analysis  program 
and  the  analytic  method  effort.   However,  the  impact  of  the 
original  simulator  forced  performance  of  the  histogram 

l  n 


generation  task  to  the  very  end  of  the  project,  due  to  the 
shifting  of  the  completion  of  the  stress  analysis  program.. 
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DESCRIPTION  OF  WORK  ON  BRIGLDl  SIMULATOR 

Nearly  a  100%  rework  of  the  BRIGLDl  simulator  was  required. 
Only  30%  of  the  original  program  was  retained  in  the  final 
version  produced  in  this  effort.   However,  that  30%  required 
significant  modification  and  correction. 

Some  general  observations  on  the  BRIGLDl  simulator  to  be  used 
for  the  generation  of  bridge  deck  loads  are  presented  herein, 
and  the  work  performed.   Specific  comments  on  programming  prob- 
lems are  also  discussed  in  this  section.   In  general,  the  simu- 
lator was  not  designed  for  the  purposes  of  structural  analysis 
of  bridges.   The  originally  generated  statistics  and  output 
data  did  not  satisfy  the  needs  of  deck  load  data  for  use  in  struc- 
tural analysis.   Mostly,  the  original  output  data  reflected  the 
traffic  generating  capability  of  the  simulator.   It  was  recommend- 
ed that  all  of  the  statistical  programs  be  eliminated  and  simple, 
user  controlled,  statistical  output  be  provided  that  would  allow 
a  user  to  verify  the  proper  generation  of  truck  traffic  for  his 
defined  bridge.   This  significantly  reduced  the  program  size  and 
had  some  effect  on  running  time,  in  that  substantial  data 
handling  was  eliminated. 

The  manner  in  which  the  data  representing  the  motion  and 
positions  of  the  vehicles  being  simulated  was  being  handled 
in  the  original  program  was  extremely  inefficient.   State-of-the- 
art  methods  in  the  use  of  pointers  and  buffers  had  not  been 
utilized.   It  was  this  particular  problem  which  caused  the 
excessive  use  of  computer  time,  e.g.,  only  a  2  to  1  compression 
ratio  of  real  time  to  computer  time. 

The  arbitrary  class  interval  size  on  truck  loads  was  not  real- 
istic in  terms  of  bridge  response  properties  and  required  change. 
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Another  significant  cause  of  excessive  computer  utilization 
was  the  motion  simulation  and  handling  of  non-truck  traffic. 
Approximately  5/6  of  all  of  the  vehicles  generated  in  the 
traffic  data  contained  in  the  program,  was  non-truck  traffic. 
In  this  case,  for  a  U.S.  301  bridge  in  Md. ,  the  motion  of  these 
vehicles  was  simulated  over  7000  feet  of  roadway,  per  lane.   On 
the  basis  of  the  data  available  for  this  particular  bridge,  truck 
traffic  occurred  on  an  average  of  one  per  two  minutes.   (It 
should  be  noted  that  the  simulator  erroneously  generated  at  three 
per  minute  in  the  originally  reported  results.)   At  the  rate 
of  one  per  two  minutes,  and  a  rough  assumption  of  two  cars  per 
truck,  only  1-1/2  vehicles  per  minute  need  be  simulated.   Further, 
only  those  non-truck  vehicles  which  would  occupy  a  portion  of 
the  bridge  deck,  at  some  time  during  the  passage  of  a  truck  over 
the  deck,  need  be  simulated.   Hence,  for  the  given  test  data 
only  approximately  50  vehicles  need  be  generated  over  a 
simulated  33  minutes.   Be  reducing  the  simulator  to  an  event 
type  simulator,  rather  than  a  continuous  time  simulator,  where 
the  event  would  be  the  generation  of  a  truck (s),  approximately 
50%  of  the  computing  time  could  be  saved,  i.e.,  based  upon  using 
a  7000  foot  roadway.   As  the  roadway  becomes  longer,  the 
savings  in  computing  time  increases. 

As  a  result  of  this  evaluation  it  was  recommended  and  agreed 
upon  that  the  program  be  modified  to  use  an  event  type 
approach,  each  event  being  the  generation  of  a  truck,  or 
platoon.   This  was  to  be  accomplished  by  forming  an  envelope 
about  each  truck,  or  platoon,  which  would  essentially  contain 
all  of  the  non-truck  vehicles  which  would  be  on  the  bridge 
deck  with  the  truck  or  platoon,  as  an  example: 

Simple  logic  could  be  used  if  a  user  inputs  maximum 
and  minimum  truck  and  non-truck  traffic  speeds.   The 
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program  need  only  generate  autos  for  approximately 
35  00  feet  of  roadway.   Any  non-truck  traffic  exceeding 
this  distance  down  the  roadway  prior  to  a  truck 
event  may  be  discarded.   Any  non-truck  traffic 
generated  after  approximately  3500  feet  behind  a  truck 
event  would  depend  on  the  subsequent  occurrence  of  a 
truck  event.   An  approach  of  this  nature  would  provide 
substantial  savings  in  computer  time. 

While  the  above  modification  was  initially  planned  for 
implementation  subsequent  developments  indicated  that  it 
was  unnecessary  in  the  program  as  it  was  being  revised.   This 
was  primarily  due  to  the  improved  vehicle  data  handling 
methodology  and  determination  that  a  shorter  roadway  provided 
more  realistic  behavior  than  the  longer  7000  foot  roadway. 

It  should  be  noted  that  the  only  traffic  data  available  in  and 
established  for  the  program  was  for  a  Md.  bridge  on  U.S.  301. 
This  data  was  not  necessarily  valid  for  use  on  any  other  bridge 
and  theoretically  applies  only  to  the  one  bridge. 

An  option  existed  in  the  original  program,  the  parametric 
bridge  length  option,  that  served  no  current  need.   It 
provided  variations  of  deck  sampling  lengths  for  the  specific 
parametric  and  distribution  data  defining  a  given  section 
of  highway.   Since  the  traffic  data  collected  for  a  given 
bridge  is  essentially  unique  to  that  bridge,  it  can  be 
established  that  it  is  not,  in  general,  suitable  to  any 
other  bridge.   Also  since  the  structural  configuration 
of  a  bridge  is  essentially  unique,  the  incremental  length- 
ening of  a  specific  bridge  has  no  meaning.   Hence,  it  was 
recommended  that  this  option  be  eliminated.   Subsequently 
the  design  of  the  revised  program  and  the  interface  philos- 
ophy to  the  structural  program  obviated  this  option  anyway. 
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rt  was  further  determined  that  with  multiple  bridge  lengths 
the  first  bridge  must  be  the  longest  when  specified  in  this  . 
option.   Otherwise,  subsequent  ends  of  bridges  would  be 
identical  to  the  first  bridge.   This  was  not  specified  in  the 
user  information  of  the  original  program  documentation. 

The  original  program  did  not  directly  generate  truck  platoons, 
nor  did  .it  attempt  to  preserve  the  actual  platoon  statistics 
for  a  given  bridge.   Any  platooning  that  may  have  occurred  was 
incidental  to  the  simulation.   A  more  realistic  basis  for 
the  generation  of  a  truck  platoon  should  be  established  in  the 
original  program  by  the  introduction  of  realistic  grades  of 
the  bridge  approaches  and  realistic  truck  traffic  behavior  on 
the  grades.   Also,  a  higher  correlation  of  weight  and  speed 
relationships  should  be  established.   However,  this  would 
amount  to  a  substantial  effort  with  questionable  results. 

A  more  realistic  approach  was  to  superimpose  a  truck  platoon 
routine  over  the  original  vehicle  generator  which  preserves 
platoon  statistics.   This  was  recommended  and  implemented. 

It  should  be  noted  that  for  the  purposes  of  this  program  a 
truck  platoon  is  defined  as  a  load  event  where  a  truck  axle  is 
on  the  bridge  or  will  go  on  the  bridge  during  an  integration 
interval.   The  event  lasts  until  there  are  no  truck  axles  on 
the  bridge. 

The  original  simulator  had  a  fairly  complex  passing  logic 
contained  in  it  for  bidirectional  and  unidirectional  two- 
lane  highways .   It  was  recommended  that  it  be  evaluated  in 
terms  of  its  significance  to  the  objectives  of  bridge 
load  prediction  and  structural  analysis.   Relative  to  this 
problem  the  whole  concept  of  roadway/ traffic  simulation 
for  the  purposes  of  bridge  structural  analysis  becomes 

15 


questionable,  in  terms  of  its  analytic  contribution.   It 
was  the  opinion  of  the  investigators  that  a  straight  forward . 
truck  traffic  generator,  which  preserves  observed  statistics 
on  a  bridge  deck,  would  have  been  far  more  meaningful  and 
certainly  more  economical  than  the  BRIGLDl  simulator. 

There  were,  besides  programming  deficiencies  and  those  indicated 
in  the  above  discussion,  certain  utilization  deficiencies  in  the 
simulator.   Some  of  these  are  defined  below: 

1.  To  satisfy  the  needs  of  modern  highway  designs  and  the 
ever  present  increasing  highway  widths,  it  was  recommended 
that  the  program  be  modified  to  handle  multiple  unidirection- 
al bridges  of  up  to  six  lanes.   However,  as  work  continued 
and  the  questionable  value  of  the  use  of  the  BRIGLDl 
simulator  increased,  it  was  decided  that  no  immediate  effort 
be  expended  to  provide  this  capability. 

2.  It  was  also  recommended  that  the  maintenance  of  the 
bidirectional  two-lane  highway  simulation  whould  be  evaluated 
as  to  its  value  to  the  Federal  Highway  Administration.   If  it 
would  be  of  little  value,  it  whould  be  eliminated.   This  would 
significantly  reduce  the  size  of  the  simulator  and  the  com- 
plexity of  its  logic. 

3.  There  was  no  minimum  vehicle  speed  provided  in  the 
original  simulator.   This  was  unrealistic  and  could  be  a  source 
of  unnecessary  computer  utilization,  and  logic  problems.   A 
minimum  vehicle  speed,  as  a  user  input  and  a  control  on  the 
vehicle  motion,  was  recommended  and  implemented. 

4.  The  traffic  generated  did  seem  to  preserve  the  dis- 
tribution statistics  which  were  used  as  input  to  the  program 
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for  short  runs.   There  was,  however,  some  evidence  of  degrada- 
tion of  the  randomness  of  the  "random"  numbers  generated  for 
longer  runs.   This  appeared  to  be  due  to  one  single  seed  used 
to  generate  all  random  data.   If  serious  use  is  to  be  made  of 
the  BRIGLD1  simulator,  especially  for  long  runs,  the  program 
should  be  modified  to  use  separate  random  number  seeds  to 
generate  distinct  functions,  i.e.,  one  each  for  vehicle  type, 
speed,  weight,  platoon  size,  etc.   Also,  for  long  runs  each 
of  these  functions  should  be  reseeded  periodically. 

5.  Additional  consideration  should  be  given  to  making 
vehicle  performance  relate  better  to  the  vehicle's  weight 
and  power.   Consideration  is  given  to  this  relationship  only 
in  the  calculation  of  acceleration  in  the  present  form  of 
BRIGLD1 . 

6.  Further,  a  significant  number  of  programming  errors 

were  initially  uncovered  in  the  preliminary  work  and  more  un- 
covered during  debug  testing  and  sensitivity  testing.   Some 
of  the  originally  discovered  errors  are  shown  below: 

PROGRAM  ERRORS 

Program     Line     Error 

INDATA      528      Read  restricted  zone  data  from  unit  6  - 

must  be  5,  6  is  the  line  printer.   Note: 
that  REZONE  has  never  been  used  or  tested. 
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PROGRAM  ERRORS  (Continued) 


Program     Line     Error 


PASPOS 


571/575   Repeat  of  same  statement,  SPDIFF=ABS 
(SPD  (MU,1)-(ILV,1)) 


PASPOS 


589 


Equation  of  the  form:  (D2/AM-.5D2/AM) -C, 
excessive  coding  here 


SORPOS 


1043 


IF (LANE (18))  3,3,2   should  be:  IF (LANE (IB)) 
3,3,2 


SORPOS 


1044      2  IBB=TB  should  be  2  IBB=IB(TB  does  not 
exist) 


SORPOS 


1078 


22   IBB=TBB+3  should  be:  22  IBB=IBB+3, 
(TBB  does  not  exist) 


SUMCK 


1228 


EQUIVALENCE  (DIFF(51),  HOUR)  50  values  of 
HOUR  are  overwritten  by  DIFF  values  51 
to  100  in  SUMLD.   However,  HOUR  is  accum- 
ulated from  call  to  call  and  must  be  saved 


BRGLD 


CONTRO 


20/48    Plot  tape  number  read  in  under  label 
NPLOT  and  used  in  program  as  IPLOT 

206      NR  not  initialized  in  the  program  before 
it  is  used. 


Specific  recommendations  that  were  made  with  regard  to  im- 
proving the  coding  of  the  simulator  included: 

1.   Subroutines  used  in  the  program  GEN,  vehicle  generation/ 
HDWAY,  TYPE,  SPEED,  WEIGHT  functions  were  extremely  trivial 
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and  should  have  been  done  in  line.   In  addition,  the  HDWAY., 
SPEED  AND  WEIGHT  funcitons  included  a  linear  interpolation 
which  was  being  solved  by  use  of  an  iterative  process  when 
it  should  have  been  a  simple  explicit  solution. 

2.  The  DTNRD  subroutine  used  by  READ  should  have  been 
done  in  line. 

3.  All  moving  of  vehicle  data  -  RENUM,  PASTES,  SORPOS, 
should  be  eliminated  and  replaced  by  pointers  for  each 
vehicle  to  vehicle  ahead  and  vehicle  behind,  and  a  buffer 
allocation  scheme  to  assign  vehicle  data  location  in  the 
tables.   All  searching  for  vehicle  ahead  and  behind  should 
be  eliminated  and  pointers  used  to  simplify  logic  in  UPDATE 
routines . 

4.  The  P0S(IV,1)  array  should  be  reduced  to  a  singly 
indexed  variable  POS(IV). 

The  details  of  the  work  performed  on  BRIGLD1  is  described 
in  the  following  subsections: 

User  Data  Input 
Traffic  Generation 

Vehicle  Data  Buffer  Handling 

Statistical  Output 

Bridge  Loading  Output 

Debug  Output 

Other  Considerations 


The  program  description,  including  listings  and  flowcharts, 
appear  in  FHWA-RD-73-43 ,  Users  Manual  for  BRIGLD1. 
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The  list  of  variables  was  corrected,  updated  and  expanded 
to  reference  all  subroutines  using  the  variable,  and  to 
include  units  of  measurement  where  applicable. 


User  Data  Input 

The  entire  data  input  portion  of  the  BRIGLD1  Simulator 
was  rewritten  with  the  user  in  mind.  The  data  input  required 
by  the  original  program  placed  an  extraordinary  burden  on  the 
engineer  using  the  program.   Also,  in  some  cases,  formats 
were  specified  such  that  significant  data  was  lost  and 
zeros  were  read.    All  single  data  variables  were  placed  in 
a  namelist  and  given  default  values  within  the  program. 

Default  values  were  assigned  on  the  basis  of  normal  design 
use  of  the  program  and  on  the  basis  of   results  from  the 
sensitivity  testing.   These  variables,  their  default  values, 
units,  and  recommended  use  are  described  in  the  Users  Manual. 
Most  of  the  tabular  data  input  was  also  redesigned  and 
data  input  forms  generated  for  this  data.   Default  values 
were  assigned  again  to  all  tabular  data.   Default  values  for  these 
variables  is  currently  set  to  the  values  used  in  the  original 
simulator  and  the  traffic  distribution  is  that  of  U.S.  301  in  Md. 
Truck  types  are  basically  the  same  as  that  used  by  the  original 
simulator  except  for  corrections  required  for  some  lengths 
and  weight  distributions.   All  double  axles  were  compressed 
to  single  axle  values  at  the  midpoints.   Default  values 
of  the  'tabular  data  were  coded  in  -the  Block  Data  Frogram 
and  can  be  easily  changed.   A  platoon  distribution  table 
was  added  to  allow  the  user  to  generate  platoons  according 
to  a  specific  distribution.   The  BRIGLD1  simulator  did  not, 
specifically  generate  or  maintain  truck  platoon  statistics. 
It  relied  totally  upon  the  adequacy  of  its  traffic  simulation. 
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A  printout  of  all  data  input  and  default  values  used  have, 
also  been  added  to  the  program. 

Traffic  Generation 

The  traffic  generation  is  functionally  the  same  as  the  original 
However,   many   small  subroutines  were  eliminated  and  coded 
"in  line"  in  the  GEN  subroutine.   Also,  several  program 
interpolation  loops  were  replaced  by  a  simple  solution  of  an 
interpolation  equation  saving  many  lines  of  code,  and  computer 
time.   The  vehicle  type  generation  was  expanded  to  include 
the  truck  platoon  generation,  previously  mentioned.   If 
a  platoon  distribution  is  not  input,  all  trucks  are  generated 
as  a  single  truck  event  and  platooning  will  take  place 
as  it  did  in  the  original  version. 

A  deficiency  in  the  time  of  vehicle  generation 
was  corrected.   The  original  version  generated 
a  time,  HDWY,  for  the  next  vehicle  to  be  generated.   This 
time  was  decremented  by  the  integration  interval  At  until 
it  reached  zero  or  a  negative  number.   The  vehicle  was  then 
generated  at  position  0  and  a  new  HDWY  generated.   This 
method  adds  an  average  of  1/2  At  to  each  HDWY  time  de- 
creasing the  amount  of  traffic  generated.   This  was  corrected 
by  adding  the  new  HDWY  to  the  residual   (0,  or  negative 
less  than  At)  and  by  placing  the  vehicle  generated  on  the 
roadway  at  the  distance  it  would  have  travelled,  i.e.,  SPD  X 
(-residual  time) . 

Currently  traffic  is  generated  only  in  the  right  lane. 
Consideration  should  be  given  to  generating  traffic  in  both 
lanes  according  to  an  input  distribution.   If  the  simulator 
is  expanded  to  more  than  two  lanes  or  if  an  extremely  short 
roadway  is  used  this  expansion  of  the  program  would  be  necessary 
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Vehicle  Data  Buffer  Handling 

The  BRIGLD1  simulator  used  excessive  computer  time  in  the 
continuous  shifting  of  all  data  related  to  a  vehicle  whenever 
its  relative  position  changed.   Also,  excessive  time  was  being 
used  in  searching  for  vehicles  ahead  and  behind.  All  of  the 
searching  and  data  shifting  was  eliminated  and  replaced  by 
a  buffer  allocation  scheme  which  kept  track  of  vehicles  ahead, 
behind,  and  the  location  of  the  related  data.   Three  tables  were 
added,  IFWD,  IBAK,  and  INDX.   These  tables  must  be  the  same 
length  as  the  data  buffers,  ITYPE,  WGT,  POS,  LANE,  etc.   The 
length  of  these  buffers  was  shortened  from  1200  to  400  words 
and  could  have  been  easily  shortened  to  100  for  the  recommended 
approach  roadway  of  1000  ft.   Some  minor  program  modification 
would  be  required  to  change  the  buffer  length.   The  program 
was  also  modified,  such  that  in  the  event  buffer  space  was  not 
available,  a  vehicle  was  generated  as  soon  as  space  was  available 
without  disturbing  the  traffic  statistics.  A  detailed  descrip- 
tion of  the  vehicle  data  buffer  allocation  methodology  is 
presented  below. 

Three  buffer  allocation  tables  were  generated  -  IFWD,  IBAK, 
and  INDX.   The  INDX  table  contains  the  pointer  to  the  vehicle 
data,  that  is,  SPD,  ACC,  POS,  WGT,  KSTAT,  ITYPE,  and  LANE. 
ACC  (acceleration  data  buffer)  has  been  added  to  the  original 
program.   Initially,  the  INDX  table  was  ordered  sequentially 
starting  with  the  second  element  in  the  table  referencing  the 
first  vehicle  data  elements.   The  first  position  of  this  table 
is  not  used. 


22 


All  vehicles  are  positionally  ordered  through  use  of  the 
IFWD  and  IBAK  pointers.   The  first  IFWD  element  points  to 
the  INDX  table  position  referencing  the  forward  most  vehicle 
on  the  road.   The  IFWD  element  referenced  by  the  first 
IFWD  element  points  to  the  INDX  table  position  referencing 
the  vehicle  directly  behind  the  most  forward  vehicle  on  the 
road,  etc. 

The  IBAK  table  strings  the  vehicles  together  in  the  other 

direction,  that  is,  from  last  to  first.   The  first  element 

of  the  IBAK  table  points  to  the  INDX  table  position  referencing 

the  last  behicle  generated  or  last  on  the  road.   The  IBAK 

element  referenced  by  the  first  IBAK  element  points  to  the 

INDX  table  position  of  the  second  last  car  on  the  road, 

and  so  forth. 

The  tables  are  initialized  as  follows: 

.IFWD (I)  =1+1,  I  =  1,  399 

IFWD (1200)  =  -1 

IBAK (I)  =1-1,  1=2,  400 

IBAK(l)  =  1 

INDX (I)  =  i  -  1,  i  =  1,  400 
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When  a  vehicle  is  generated,  its  data  is  assigned  to  the  position 
referenced  by  the  INDX  value  at  the  forward   (IFWD)   reference 
of  the  last  vehicle  generated,  i.e.,  at  INDX  (IFWD  (IBAK(1)l  )  .   The 
last  vehicle  generated  reference  is  changed: 

BAK(l)  =  IFWD (IBAK (1)) 

When  a  vehicle  at  POS(INDX(J))  is  passed  by  its  IFWD  referenced 
vehicle,  POS (INDX (IFWD (J) ) ,  the  pointers  in  the  INDX  table 
at  position  J  and  at  IFWD (J)  are  interchanged. 

Any  number  of  vehicles  in  a  string  from  J  to  K  may  be  elim- 
inated by  inserting  the  string  at  the  position  of  the  next 
vehicle  to  be  generated  and  connecting  the  J-l  and  K-l  vehicles 
by  changing  the  following  forward  and  backward  references: 

IFWD(K)  — *•  TEMP 
IFWD(IBAKU)  — *■   IFWD(K) 
J  —►  IFWD  (IBAK  (1)) 
TEMP  —►IFWD  (IBAK  (J)) 
IBAK(TEMP)  — *  IBAK (IFWD (K) ) 

IBAK (J)  *IBAK(TEM|>) 

IBAK(l)  ►IBAK  (J) 

J  may  be  equal  to  K   (one  vehicle  only)   and  the  string  may 
appear  anywhere  including  the  forward  or  backward  ends  of  the 
traffic. 

The  following  examples  are  offered  to  illustrate  the  use  of 
the  IFWD,  IBAK,  and  INDX  tables   (limited  here  to  a  length 
of  16).   The  positional  and  lane  vehicle  data  tables  are 
included  in  the  examples,  see  Figures  3  through  6. 
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Initialized  tables: 


I 

IFWD 

IBAK 

INDX 

POS 

LANE 

1 

2 

1 

0 

- 

- 

2 

3 

1 

1 

- 

- 

•3 

4 

2 

2 

- 

- 

4 

5 

3 

3 

- 

- 

5 

6 

4 

4 

- 

- 

6 

7 

5 

5 

- 

- 

7 

8 

6 

6 

- 

- 

8 

9 

7 

7 

- 

- 

9 

10 

8 

8 

- 

"- 

10 

11 

9 

9 

- 

- 

11 

12 

10 

10 

- 

- 

12 

13 

11 

11 

- 

- 

13 

14 

12 

12 

- 

- 

14 

15 

13 

13 

- 

- 

15 

16 

14 

14 

- 

- 

16 

-1 

15 

15 

- 

- 

■ 

Figure  3.   Vehicle  Buffer  Allocation  Initial 
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Now  we  will  generate  12  vehicles  on  the  road: 


[10] 


[  5] 


U2T  TUT 


[  9] 


T~8T    [  7]  T~6T 


HI  T^T  T^T  [   H 


100 


2Q0 


300 


400 


500 


600 


700 


Updated  Tables: 


I 

IFWD 

IBAK 

INDX 

POS 

LANE 

1 

2 

13 

0 

700 

1 

2 

3 

1 

1 

650 

1 

3 

4 

2 

2 

600 

1 

4 

5 

3 

3 

525 

1 

5 

6 

4 

4 

500 

2 

6 

7 

5 

5 

400 

1 

7 

8 

6 

6 

350 

1 

8 

9 

7 

7 

275 

1 

9 

10 

8 

8 

200 

1 

10 

11 

9 

9 

150 

2 

11 

12 

10 

10 

100 

1 

12 

13 

11 

11 

25 

1 

13 

14 

12 

12 

- 

- 

14 

15 

13 

13 

- 

- 

15 

16 

14 

14 

- 

- 

16 

-1 

15 

15 

™ 

*■ 

Figure  4.  Vehicle  Buffer  Allocation  Vehicles  Added 
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Next  vehicle  No.  10  will  pass  vehicle  No.  9  and  vehicles  1  and  2 
will  be  deleted  as  they  leave  the  roadway: 


[8]       [5] 

[11] 

[12] 

E  4] 

[  9]  [10]  [7]  [6] 

[  3] 

100 


200 


300 


400 


500 


600 


700 


Updated  Tables: 


I 

IF  WD 

IBAK 

INDX 

POS 

LANE 

1 

4 

13 

0 

- 

- 

2 

3 

13 

1 

- 

- 

3 

14 

2 

2 

650 

1 

4 

5 

1 

3 

575 

1 

5 

6 

4 

4 

525 

2 

6 

7 

5 

5 

500 

1 

7 

8 

6 

6 

450 

1 

8 

9 

7 

7 

425 

2 

9 

10 

8 

8 

350 

1 

10 

11 

9 

10 

400 

1 

11 

12 

10 

9 

275 

2 

12 

13 

11 

11 

200 

1 

13 

2 

12 

12 

- 

- 

14 

15 

3 

13 

- 

- 

15 

16 

14 

14 

- 

- 

16 

-1 

15 

15 

- 

- 

Figure  5.  Vehicle  Buffers  Modified 
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Three  more  vehicles  are  generated,  (1,  2,  and  13),  vehicle  5  will 

pass  vehicle  4,  and  vehicle  11  will  pass  vehicle  9.   Then  vehicles 

* 
10,  8,  and  7  are  eliminated  as  non  truck  influencing: 


[11] 


[  5; 


[13] 


[2]  [  1] 


[12] 


9] 


e     100 

Updated  tables 


200 


300 


400 


6] 


500 


4] 


600 


Figure  5.  Vehicle  Buffers  Modified (Cont. ) 


[  3] 


700 


I 

IF  WD 

IBAK 

INDX 

POS 

LANE 

1 

4 

14 

0 

175 
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2 

3 

13 

1 

125 

1 

3 

14 

2 

2 

700 

1 

4 

5 

1 

3 

600 

1 

5 

6 

4 

5 

650 
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6 

7 

5 

4 

550 

1 

7 

11 

6 

6 

- 

- 

8 

9 

14 

7 

- 

- 

9 

10 

8 

8 

350 

1 

10 

15 

9 

10 

- 

- 

11 

12 

7 

11 

375 

2 

12 

13 

11 

9 

250 

1 

13 

2 

12 

12 

50 

1 

14 
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3 

11 

- 

- 

15 

16 

10 

14 

- 

- 

16 

-1 

15 

15 

- 

- 

The  elimination  of  car  traffic  determined  not  to  be  in  the  range 
of  influence  of  truck  traffic  has  not  been  implemented  in  view 
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Bidirectional  traffic  is  handled  by  splitting  the  400  word 
buffer  allocation  tables  into  two,  200  word  tables: 

One  Direction  Bidirectional  Traffic 

Forward  Forward       Reverse 

IFWD(400)  IFWD(200)     JFWD(200) 

IBAK(400)  IBAK(200)     JBAK(200) 

INDX(400)  INDX(200)     JNDX(200) 

The  first  half  of  the  data  buffers  (1-200)  are  then  allocated 
to  the  forward  direction  traffic  and  the  second  half  (201-400) 
are  allocated  to  the  reverse  direction  traffic.   Reverse  traffic 
direction  tables  operate  then  in  the  same  manner  as  the  forward 
traffic  direction  tables. 

Vehicle  Motion 

Vehicle  motion  is  essentially  the  same  as  is  described 
in  (1)  as  far  as  passing  and  acceleration  decision 
logic  is  concerned.   However,  major  changes  were  made  to  this 
portion  of  the  program.  The  original  version  started  with  the 
foremost  vehicle,  deciding  its  passing  status,  acceleration, 
and  changing  its  position  immediately.   Status  and  acceleration 
decisions  for  all  following  vehicles  were  made  on  the  basis  of 
all  vehicles  forward  of  it  integrated  to  the  next  time  step, 
while  all  vehicles  behind  were  in  the  present  time  step.   This 
was  changed  in  the  "Bridge  Load  Generator"  to  a  two  step  process, 
decision  and  integration,  where,  starting  with  the  foremost 
vehicle,  status  and  acceleration  decisions  were  made  for  all 
vehicles  at  positions  in  the  same  time  step.   Lane  changes  are, 
however,  still  made  during  the  decision  phase.   All  vehicles 
are  then  integrated  forward  on  the  basis  of  acceleration  and 
the  old  speed.   This  required  the  addition  of  an  acceleration 
table,  ACC.   Speed  is  then  changed  only  by  acceleration. 


of  the  current  compression  of  computer  running  time  relative  to 
simulation  time  and  the  recommendation  of  a  short  roadway 
(approximately  1000  ft) ,  as  indicated  earlier.   The  example  does 
illustrate  how  to  eliminate  one  or  more  vehicles  from  the 
middle  of  the  roadway. 
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Significant  coding  errors,  in  addition  to  those  described 
earlier,  were  found  in  the  vehicle  motion  portion  of  the 
program  and  were  eliminated.   A  debug  routine,  GRAPH,  was 
generated  in  order  to  visualize  the  motion  of  vehicles  on 
the  roadway.   This  printout  is  obtained  along  with  allocation 
buffers  and  position,  acceleration,  status  and  speed  data  for 
each  vehicle  by  setting  the  DBUG  parameter  "True"  during  pro- 
gram data  input.   There  was  no  way  to  "see"  or  verify  traffic 
motion  and  behavior  in  the  original  simulator. 

There  are  still  some  problems  with  bidirectional  traffic  with 
occasional  "crashes"  of  passing  cars  not  completing  the  pass 
in  time  where  extremely  dense  traffic  is  encountered,  as  in 
the  originally  defined  traffic  data.   In  general,  however,  the 
bidirectional  traffic  behaves  reasonably  well,  but  was  of 
questionable  value  to  the  objectives  of  this  effort. 

Statistical  Output 

The  statistical  output  and  statistic  gathering  portions  of 
the  original  simulator  were  entirely  replaced. 

The  following  statistics  and  data  are  gathered  by  the  revised 
"Bridge  Load  Generator  (BRIGLD1)": 

1.  Total  vehicles  generated. 

2.  Simulated  time  -  starting  from  the  time  the  first 
vehicle  enters  the  bridge. 

3.  Platoon  distribution  of  1  through  9  and  10  and 
above.   This  statistic  is  gathered  on  the  bridge  and  at  the 
time  the  platoon  is  generated,  both  forward  and  reverse 
directions.   A  truck  is  part  of  a  platoon  on  the  bridge 

if  it  enters  the  bridge  during  a  time  frame  when  there  is 
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•still  all  or  part  of  another  truck  on  the  bridge  or  during 
the  next  following  time  step. 

4.  Type  distribution  -  this  statistic  is  gathered 
as  the  vehicle  enters  the  bridge. 

5.  Load  distribution  -  this  statistic  is  gathered 

as  the  vehicle  enters  the  bridge  and  is  distributed  according 
to  class  intervals  specified  by  the  user. 

6.  Statistical  output  is  provided  at  the  conclusion 
of  a  run,  and  at  every  simulated  hour  of  the  run. 


Bridge  Loading  Output 

The  Bridge  Loading  Output  portions  of  the  program  are 
additions  to  the  original  Simulator.   Bridge  loading  data  is 
output  whenever  there  is  either  a  truck  on  the  bridge  or  one 
is  entering  the  bridge  during  the  next  time  frame.   This  data 
is  associated  with  a  loading  event  number.   Three  types  of 
data  are  output: 

1.   General 


Simulation  time 

Load  event  number 

Integration  time  increment 

Load  Data  -  by  axle 

Number  of  axles  on/going  on  the  bridge 

Position  relative  to  the  start  of  the  bridge 

Lane 

Axle  weight  -  may  be  a  double  axle 

Distance  the  axle  will  travel  during  this 

integration  step 
Acceleration 
Vehicle  identification 
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Number  of  trucks  on/going  on  the  bridge  since  start 

of  this  event 
Truck  Type 
Truck  Weight 
Truck  Speed 
Lane 
Time  of  truck  entering  the  bridge 


The  position  of  a  vehicle  axle  is  referenced  to  zero  at 
the  beginning  of  the  pseudo  bridge. 

Debug  Output 

A  printout  of  the  following  data  is  available  by  setting 
the  DBUG  parameter   to  TRUE  during  program  data  input. 

-Vehicle  generation  data 

Printout  of  buffer  allocation  tables  and 

vehicle  data  buffers 
Printout  of  the  roadway  and  bridge  up  to  6000  ft 

and  two  lanes  with  vehicle  numbers  showing 

positionally  on  the  roadway 
Printout  of  bridge  loading  data 

An  example  of  the  debug  output  is  provided  in  Appendix  A. 

Other  Considerations 

The  use  of  multiple  bridge  lengths  was  eliminated  as  being  un- 
necessary. However,  bridge  loading  data  should  be  generated  for  the 
longest  anticipated  bridge  and  is  directly  useable  for  any 
shorter  length  bridge. 
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SENSITIVITY  TESTING  OF  LOAD  GENERATOR 

The  first  task  required  by  the  contract  was  to  perform 
sensitivity  testing  of  the  input  parameters  of  the  Government 
furnished  bridge  loading  program.   As  was  indicated  in  the  fore- 
going, it  was  not  possible  to  either  use  this  program  as  a 
load  generator  to  a  stress  program  or  to  generate  traffic  data. 
Further,  most  of  the  parameters  required  in  the  input  were 
either  unnecessary  or  related  to  the  control  of  meaningless 
statistics  generated  by  the  original  program.   Subsequent 
to  performing  the  major  modifications  discussed  previously, 
it  was  then  possible  to  identify  meaningful  parameters  for 
the  performance  of  sensitivity  testing,  or  at  least  to  be 
considered  for  sensitivity  testing.   It  should  be  noted  that 
the  technique  of  sensitivity  testing  is  a  process  of  evaluating 
the  effects  of  variables  in  a  mathematical  or  logical  model 
whose  relationships  in  the  model  cannot  be  explained  analytically 
This  is  the  case  with  representations  that  utilize  engineering 
fudge  factors  or  stochastic  processes. 

Potential  Input  Parameter  Candidates 

The  input  parameters  which  remained  in  the  load  generator 
subsequent  to  modification,  were  not  traffic  statistics, 
and  which  potentially  could  affect  the  output  of  the  generator, 
were  defined  as  follows: 

NTH  Number  of  load  samples,  i.e.,  number  of  input 

cases  to  be  run;  it  defaults  to  1. 

MD  Maximum  number  of  vehicle  types  over  a 

simulation  run,  i.e.,  all  NTH  spaces  per 
run;  it  defaults  to  11,  with  a  maximum  of  20 
allowed. 
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LT 


Number  of  class  intervals  for  weight  histo- 
grams; defaults  to  12  equal  intervals,  with 
a  maximum  of  99  intervals. 


NL 


Number  of  lanes,  i.e.,  bridge  width  in  lanes. 
(No  meaning  of  a  1  lane  utilization  presently. 
However,  should  be  retained  for  future  ex- 
pansion to  more  than  one  lane.)   It  defaults 
to  2. 


ND 


Number  of  directions  of  traffic  flow,  if  1 
then  2  lane  unidirectional  flow;  if  2  then 
2  lane  bidirectional  flow;  defaults  to  1. 


NZ 


Number  of  restricted  or  graded  zones,  NZ<5 , 
each  are  defined  below;  defaults  to  0. 


SIMTIM 


Length  of  simulation  time  for  given  simu- 
lation case,  minutes;  defaults  to  program 
abort. 


DELTIM 


Interval  of  integration  for  solution  of  traffic 
motion,  sees. 


BRPOS 


Center  of  bridge  relative  to  the  beginning 
of  the  roadway,  establishes  roadway  approach 
length  to  the  bridge.   This  value  defaults 
to  1100  ft. 


TALINC 


Class  interval  for  weight  histograms  collected 
over  simulation  times,  lbs;  defaults  to 
8000  lb.  intervals. 
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ACCEL 


Maximum  acceleration  allowed  for  any  vehicle, 

2  2 

fps  ;  defaults  to  15  f/sec  . 


SPDLIM 


Highway  speed  limit,  mph;  defaults  to  65 
mph. 


EXSPD 


Excess  speed  permitted  for  passing,  mph; 
defaults  to  15  mph. 


BRLEN 


Bridge  length  for  case  to  be  run,  normally 
should  be  of  maximum  length  of  interest. 
The  structural  program  will  allow  the  input 
of  the  actual  bridge  lengths  to  be  analyzed 
which  will  be  less  than  or  equal  to  this 
value  for  the  section  of  highway  of  interest; 
defaults  to  200  feet. 


SDFAC 


SAFDIS 


Speed  factor  to  determine  maximum  decision 

point  for  passing,  f t-hr/mi ,  defaults  to  15  which 
provides  approximately  one  vehicle  length  per 
10  mph  speed. 
Minimum  distance  between  vehicles,  ft.; 

defaults  to  10  ft. 


TRKLIM 


Truck  speed  limit,  mph.;  defaults  to  55  mph. 


SPDMIN 


Minimum  vehicle  speed,  mph.;  defaults  to  40 
mph. 


In  evaluating  the  input  parameters,  two  basic  considerations 
were  used  as  guidelines: 
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1.  Preservation  of  the  input  traffic  population 
statistics,  on  the  bridge  deck,  and 

2.  Minimization  of   computer   time  for  generating 
synthetic  deck  loads. 

With  these  factors  in  mind,  an  initial  analysis  of  all  input 
parameters  was  performed  to  minimize  the  sensitivity  testing 
A  subsequent  computer  based  sensitivity  test  was  performed 
on  those  parameters  which  remained  in  doubt  as  to  their 
affects  upon  generated  loads. 

Analysis  of  the  Candidate  Parameters 

The  initial  analysis  of  the  above  specified  input  parameters 
indicated  the  following: 

1.  NTH  is  a  case  count  parameter  that  does  not  af- 
fect generated  loads  and  did  not  require  sensitivity 
testing. 

2.  MD  is  basically  a  traffic  parameter  based  upon 
the  categorization  of  truck  types,  all  autos  form 

one  category,  and  is  totally  dependent  upon  the  field 
collected  population  sample  or  an  hypothesized  sample. 
It  did  not  require  sensitivity  testing. 

3.  LT  is  an  output  control  value,  i.e.,  it  controls 
the  number  of  class  intervals  in  the  weight  histograms, 
but  in  no  way  affects  the  generated  load  data  to  be  used 

by  the  stress  program.   It  did  not  require  sensitivity 
testing. 
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4.  NL  is  the  bridge  width  in  lanes.   The  load 
generator  is  limited  to  a  maximum  of  2  lanes  of  either 
unidirectional  traffic  flow  or  bidirectional  flow.   At 
present  a  one-lane  bridge  has  no  meaning  to  the  generator. 
Hence,  it  did  not  require  sensitivity  testing. 

5 .  ND  is  the  number  of  traffic  flow  directions. 
It  did  require  testing  for  functional  performance  and  a 
gross  form  of  sensitivity  testing ,  i.e., was  there  a  sig- 
nificant difference  between  unidirectional  flow  simulation 
and  bidirectional  flow  simulation. 

6.  NZ  is  the  restricted  zone  control  number.   This 
parameter  definitely  required  evaluation.   However,  the 
related  parameters  e.g.,  restricted  zone  length,  grade,  etc., 
had  to  be  simultaneously  evaluated  in  terms  of  their  impact 
upon  the  generated  data  and  proper  functional  operation. 

7.  SIMTIM  is  the  length  of  simulated  time  over  which 
loads  are  continually  generated.   This  variable  required 
evaluation  in  terms  of  its  impact  upon  the  preservation  of 
the  stochastic  parameters  defining  the  traffic  for  any  given 
highway.   Desirably, the  best  value  was  that  which  adequately 
preserved  the  stochastic  parameters  and  was  the  smallest 

in  magnitude.   Assuming  normal  distribution  of  each  type  of 
truck,  in  order  to  achieve  a  90%  confidence  interval  of  the 
occurrance  of  the  lowest  incidence  vehicle  type,  approximately 
a  5-hour  simulation  time  is  necessary  for  the  stochastic 
parameters  defined  for  the  MD  301  bridge  data  contained  in  the 
original  program. 
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8.  DELTIM  is  the  integration  internal  used  in  the 
solution  of  the  motion  of  the  vehicles  by  the  load  generator. 
A  conservatively  safe  value  of  DELTIM  for  the  purpose  of 
maintaining  proper  vehicle  dynamics  is  .5  seconds. 

A  one  second  interval  appears  to  provide  a  nominally  safe 
value.   However,  the  error  in  the  simple  numerical  integration 
utilized  in  the  modified  program,  Euler ' s  and  Trapezoidal, 
due  to  internal  size  required  evaluation.   The  largest  value 
which  maintained  the  accuracy  necessary  and  consistent  with 
the  smaller  values  should  be  selected.   Two  effects 
were  considered: 

a.  Error  in  calculating  vehicle  motion  dynamics 

b.  Variation  in  vehicle  logical  behavior  which 
affected  the  generated  bridge  loads. 

9.  BRPOS  establishes  the  length  of  the  roadway  used 
for  traffic  simulation  by  the  load  generator.   It  is  actually 
measured  from  the  beginning  of  the  roadway  to  the  center 

of  the  bridge.   This  parameter  required  evaluation  to  deter- 
mine its  affect  upon  the  generated  loads.   The  smallest 
value  which  best  preserves  the  statistical  para- 
meters defining  the  truck  and  traffic  population/    input 
to  the  program  should,  was  selected  for  use. 

10.   TALINC  is  an  output  control  parameter  which 
establishes  the  size  of  the  weight  histogram  class  intervals 
which  are  used  to  collect  the  weights  of  the  vehicles 
generated  on  the  bridge  deck  over  the  period  of  simulation. 
This  parameter  in  no  way  affected  the  generated  load  data 
for  the  stress  program  and  is  merely  a  user  preference  for 
weight  histogram  output. 
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11.   ACCEL  is  a  limiting  value  on  the  acceleration 

allowed  any  vehicle  synthesized.   A  maximum  default  value 

2 

of  15  fps   was  established  based  upon  realistic  short  term 

accelerations  developed  by  high  performance  automobiles. 
It  was  not  considered  necessary  to  test  this  parameter 
since  a  realistic  limit  can  be  established  from  real  ve- 
hicle performance.   There  is  an  assumption  inherent  in  this 
decision,  i.e.,  that  the  generator  adequately  synthesizes 
the  behavior  and  motion  of  highway  traffic. 

Many  vehicle  behavior  characteristics  have  been  simplified 
or  ignored  in  BRIGLD1 . 


12.  SPDLIM,  this  value  should  be  established  by 

a  user  as  a  function  of  the  highway  of  which  his  bridge  is 
a  part.   A  default  value  of  65  mph .  has  been  established  for 
this  parameter.   No  sensitivity  testing  was  considered 
necessary  for  the  same  reasons  as  stated  in  11.,  above. 

13.  EXSPD  is  the  amount  of  excess  speed  allowed 
above  SPDLIM  for  passing.   Again,  this  parameter  is  a  function 
of  the  particular  highway  and  should  be  established  by  a 
user.   A  default  value  of  15  mph.  was  established  in  the 
program  and  no  sensitivity  testing  was  considered  necessary 
for  the  same  reasons  as  stated  in  11.,  above. 

14.  BRLEN  is  the  length  of  the  pseudo  bridge  deck 
on  which  the  generated  traffic  loads  are  captured,  struc- 
tured and  passed  to  the  stress  program.   This  value  should 
be  selected  by  a  user  on  the  basis  of  potential  re-use 
for  the  section  of  highway  which  is  represented  £>y  the 
statistical  traffic  data  input  to  the  bridge.   It  should  be 
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at  least  as  long  as  the  longest  possible  span,  simple  or 
continuous,  whichever  is  greater,  which  may  be  needed  on 
the  represented  section  of  highway.   The  data  generated 
by  the  load  generator  for  one  span  of  a  multi-span  bridge 
may  be  re-used  by  the  stress  program  for  all  spans  of  a 
given  bridge  and  for  spans  of  nearly  all  bridges  on  the 
same  highway,  provided  the  input  statistical  traffic  data 
is  valid,  since  the  stress  program  separately  discerns 
its  own  load  events  for  any  given,  and  specific,  bridge 
length.   A  default  value  of  200  feet  was  established  for 
this  parameter.   Because  of  the  foregoing  no  sensitivity 
testing  was  necessary,  i.e.,  it  in  no  way  affected  the  loads 
input  to  the  stress  program. 

15.  SDFAC  is  the  safe  distance  factor  between 
following  vehicles.   It  is  based  upon  the  use  of  a  gap 
between  vehicles  as  a  function  of  a   following  vehicle's 
speed..  Using  the  rule  of  thumb  of  one  vehicle  length  per 
10  mph  of  speed,  a  default  value  of  15  was  established  for 
this  parameter.   The  actual  value  for  this  rule  of  thumb 

is  14-2/3/sec.   However,  it  is  fairly  well  established  that 
the  desired  safe  following  rule  of  thumb  is  not,  in  general, 
maintained.   Hence,  it  was  decided  that  the  effect  of  this 
parameter  upon  the  generated  loadings  should  be  determined 
so  that,  if  any  significant  dependency  existed,  a  user 
would  have  a  basis  of  selecting  an  input  value  for  this 
parameter  as  a  function  of  traffic  behavior  on  any  given 
highway . 

16.  SAFDIS  is  the  absolute  minimum  distance  that  a 
following  vehicle  may  approach  a  leading  vehicle.   At 
reaching  this  distance  with  no  other  actions  having  occurred 
the  program  forces  the  following  vehicle  to  either  pass 
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or  slow  to  the  same  speed  as  the  leading  vehicle.   This 
parameter,  is  a  lower  limit  and  a  function  of  traffic 
behavior  on  the  given  highway.   A  default  value  of  10 
feet  was  established  and  no  sensitivity  testing  was  felt 
necessary  for  this  parameter. 

17.  TRKLIM  is  the  speed  limit  imposed  upon  trucks, 
actual  or  theoretical,  for  a  given  highway.   If  this 
value  is  set  at  the  legal  limit,  then,  EXSPD  should  be 
set  realistically,  i.e.,  on  the  basis  of  actual  behavior 
of  truck  traffic  speedwise.   A  default  value  of  55  mph 
was  established  for  this  parameter  and  sensitivity  testing 
was  not  felt  necessary. 

18.  SPDMIN  is  the  nominal  minimum  speed  anticipated 
or  recorded  for  the  given  highway,  adjacent  to  the  bridge. 
It  is  established  within  the  load  generator  program  as  an 
absolute  minimum  speed  at  which  vehicles  must  travel.   This 
was  necessary  because  the  BRIGLDl  simulator  used  an  al- 
gorithm for  vehicle  speed  which  appeared  to  converge 

toward  zero  rather  than  the  statistical  nominal  speed  assigned 
to  a vehicle  as  a  steady  state  speed.   Sufficient  inter- 
action of  vehicles  over  a  long,  synthetic,  roadway  fails  to 
continue  at  a  high  enough  level  for  the  interaction  portion 
of  the  algorithm  to  force  a  maintenance  of  speed  and  since 
the  algorithm  does  not  converge  upon  the  steady  state 
speeds  assigned,  a  trend  toward  zero  speed  occurs.   It  is 
essential  that  a  value  of  SPDMIM  be  input  to  prevent  this 
problem.   Inhibition  of  this  tendency  is  also  accomplished 
by  using  a  minimum  acceptable  roadway  approach  to  the  bridge. 
A  default  value  of  40  mph.  was  established  in  the  program. 
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Parameters  to  be  Tested 

As  a  result  of  the  foregoing  analysis, the  set  of  para- 
meters requiring  sensitivity  testing  was  established  as: 

1 .  ND  -  Number  of  flow  directions 

2.  NZ  -  Number  of  restricted  zones 

3.  SIMTIM  -  Simulation  Time  (Sec) 

4.  DELTIM  -  Time  increment  (Sec) 

5.  BRPOS  -  Length  of  approach  roadway  (plus  1/2 
bridge  length) (feet) 

6.  SDFAC  -  Safe  distance  factor 

The  implications  of  testing  the  above  parameters  was  as 
follows : 

1.  The  ND  test  was  primarily  a  functional  test  and 
provided  evidence  of  the  ability  of  the  BRIGLD1  generator 
to  function  in  either  a  unidirectional  or  bidirectional 
manner . 

2.  NZ  is  merely  a  control  parameter  and  did  not 
directly  imply  anything  of  consequence.   However,  grade 
variation  testing  would  indicate  the  ability  of  the 
generator  to  impose  the  effects  of  grades  adjacent  to  a 
bridge.   This  test  would  merely  indicate  whether  the  load 
generator  provided  realistic  effects  from  adjacent  grades. 
The  no-passing  zone  effect  was  of  small  consequence  and 
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had  significant  meaning  only  in  terms  of  bidirectional 
traffic  flow. 

3.  SIMTIM  required  evaluation  to  determine  the 
nominal  minimum  stable  simulation  time,  which  was  defined 
as  that  simulation  time  sufficient  to  preserve  the  input 
statistical  parameters. 

4.  The  testing  of  DELTIM  would  indicate  the  most 
economic  integration  interval  for  calculating  traffic 
motion. 

5.  The  most  economic,  and  stabilizing,  highway  approach 
length  which  preserved  prescribed  traffic  statistics  would 

be  indicated  by  the  testing  of  BRPOS . 

6 .  The  testing  of  SDFAC  would  provide  indication 
of  the  sensitivity  of  this  parameter  and  establish  an 
estimate  of  a  worst  case  value  to  utilize  the  default  input 
purposes . 

Parameter  Test  Order 

In  addition,  there  was  a  very  definite  interdependence 
of  the  above  parameters.   Hence,  those  parameters  most 
dependent  upon  the  others  in  the  set  would  be  tested  last. 
Those  least  dependent  upon  the  remainder  of  the  set  were 
tested  first.   The  order  of  the  testing  performed  was 
as  follows : 

1.   BRPOS  to  determine  minimum  stable  load  generation 
approach  length. 
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2.  SIMTIM  to  determine  minimum  characterizing  sim- 
ulation time. 

3.  DELTIM  to  determine  maximum  practical  interval 
of  motion  integration. 

4.  SDFAC  to  determine  worst  case  value  for  default 
purposes  and  significance  of  this  variable. 

5.  NZ  to  determine  significance  and  validity  of 
synthesized  grade  effects. 

6.  Nr>    to  determine  the  sensitivity  of  this  vari- 
able in  the  load  generator,  and  its  functional  behavior. 


Sensitivity  Test  Specification 

As  a  result  of  all  of  the  analysis  previously  described, 
a  specification  of  the  "sensitivity"  testing  to  be  per- 
formed was  developed  at  the  point  in  the  project  where  the 
BRIGLDl  simulator  was  as  well  debugged  as  could  be  deter- 
mined for  unidirectional  flow  traffic. 

The  tests  were  specified  as  follows: 

1.   ND  -  Number  of  directions  of  traffic  flow 

a.  Test  bidirectional  flow  to  assure  valid 
simulation. 

b.  Run  one  bridge  configuration  for  uni-  and 
bidirectional  traffic  flow,  based  upon  the  same  highway 
traffic  data  base.   This  was  necessary  due  to  the 
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form  of  BRIGLD1.   It  was  anticipated  that  the  bi- 
directional flow  would  produce  a  higher  incidence  of 
loading  events  and  tend  toward  a  worst  case  loading, 
as  compared  to  the  unidirectional  flow.   Thus,  a  synthetic 
bidirectional  flow  would  produce  more  significant  input 
for  the  structural  program.   Constant  approach  lengths 
and  the  best  set  of  parameters  from  previous  tests  were 
utilized  in  this  test. 

2.   NZ  -  NZ  would  be  tested  for  one  bridge  con- 
figuration, for  a  minimum  valid  simulation  period,  and  for 
the  following  cases,  with  the  best  set  of  parameters  from 
previous  tests: 


a.   For  functional  and  accuracy  testing 
(1)  Unidirectional  2  lane  flow  for 


(i)     NZ=0 

(ii)    NZ=5  grades 


(2)  Bidirectional  flow  for  2  lanes  for 

(i)     NZ=0 

(ii)    NZ=5  (as  above) 

(iii)   NZ=1  for  z(l)  =0 
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(a) 


1950' 


50' 


0%  grade 


(b) 


1200 


k750' 


50 


10%  grade 


(c) 


1900 


50 


5%  grade 


50' 


(d) 


1200' 


750' 


10%  grade 


(e) 


30 


191 


50 


5%  grade 


(f) 


2000 


control  core 
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3.  SIMTIM  -  SIMTIM  would   be  either  tested  or 
analytically  determined  such  that  the  minimum  simulation 
time  preserved  the  input  traffic  statistics.   Probability 
analysis  of  available  data  indicated  approximately  5  hours 
to  provide  a  90%  confidence. 

4.  DELTIM  -  DELTIM  would  be  tested  as  follows: 

For  maximum  anticipated  bridge  length,  run, 
250  ft., 

d  O        •  D   S6C • 

b.  1.0  sec. 

c.  1.5  sec. 

d.  2.0  sec. 

It  was  anticipated  that  1.0  sec.  integration 
intervals  will  be  near  optimum  for  all  conditions 
Use  a  constant  stable  approach  length  and  a 
valid  simulation  time,  as  indicated  from  prior 
testing  would  be  used. 

5.  BRPOS  -  BRPOS  would  be  tested  for  one,  nominal, 
bridge  length  and  allowed  to  vary  to  determine  statistical 
variations  of  truck  loads  approaching  on  the  bridge.   This 
would   be  done  as  follows: 

For  a  fixed  SIMTIM,  At  =  1.0  sec,  and  uni- 
directional traffic  flow,  as: 

a.  26,400  ft. ,  5  miles 

b.  10,560  ft. ,  2  miles 

c.  4,000  ft. 
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d.  1,000  ft. 

e.  200  ft. 

6.   SDFAC  -  SDFAC  would  be  tested  for  unidirectional 
traffic  flow,  constant  approach  length,  constant  bridge 
length,  for  constant  highway,  for  values  of 

a.  5 

b.  10 

c.  15  (default  value) 

d.  20 

e.  25 

using  the  best  set  of  parameters  from  prior  testing. 

Conduct  and  Results  of  Tests 

In  accordance  with  the  previously  presented  test  specification, 
the  sensitivity  tests  were  conducted,  after  review  and 
concurrence  by  the  Government's  Contract  Monitor  of  the 
specification.   The  testing  and  the  results  are  described 
within  this  section  of  the  report. 

It  should  be  noted  that  this  testing  effort  was  the  first, 
essentially,  production  type  use  of  the  program  and  numerous 
difficulties  were  encountered  in  achieving  a  "production" 
mode.   These  problems  were  of  the  nature  of  structuring  the 
JCL  to  accommodate  "production"  runs  of  the  program.   Two 
problems,  bugs,  were  encountered  during  the  testing,  i.e., 

1.   Failure  of  the  program  to  run  with  extremely 
short  roadway  approaches.   This  was  quickly  and  easily 
corrected. 
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2.   Difficulty  in  getting  reasonably  valid  bidirectional 
traffic  flow.   This  was  an  unanticipated  problem.   Since 
the  original  BRIGLDl  simulator  was  supposed  to  provide 
valid  synthesis  of  bidirectional  flow.   However,  this  was 
consistent  with  the  difficulties  encountered  in  making  the 
simulator  simulate,  as  discussed  earlier  in  this  report. 

BRPOS  Testing  -  Much  emphasis  was  laid  upon  the  need  to 
provide  a  sufficient  roadway  length,  i.e.,  the  approach 
(es)  to  the  bridge  deck,  in  the  original  BRIGLDl  report  (1) 
The  basis  of  this  emphasis  was  the  need  to  provide  adequate 
simulation  time,  prior  to  bridge  entry,  to  develop  a  stable 
traffic  behavior,  i.e.,  to  overcome  the  effects  of  the 
pseudo-random  generator  algorithm.   However,  a  detailed 
evaluation  of  the  pseudo-traffic  behavior  indicates  that 
for  the  traffic  population,  and  characteristics  utilized, 
i.e.,  the  same  as  the  original  work,  Md .  301,  the 
traffic  synthesis  appeared  realistic  after  a  very  short 
distance.   The  use  of  long  roadway  approaches  earlier 
evidenced  the  tendency  of  the  simulator  to  converge  to  a 
zero  speed  and  resulted  in  the  establishment  of  a  minimum 
speed  to  inhibit  this  tendency.   However,  the  simulator 
still  behaves  in  this  fashion  and  with  long  roadway  approaches 
will  tend  for  all  vehicles  to  converge  on  first  the  steady 
state  speed,  and  second,  on  the  minimum  speed.   Further, 
dynamic  interaction  of  the  vehicles  appears  to  damp  out. 
This  is  definitely  a  handicap  because  of  the  established 
relationship  of  speed  and  weight  to  bridge  structural 
responses.   Further,  the  preservation  of  the  given  highway 
traffic  statistics  fails  to  be  preserved  for  long  roadway 
approaches. 

The  evaluation  criterion  was: 
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1.  Preservation  of  original  traffic  statistics 
on  the  bridge  deck,  i.e.,  captured  loading  data, 

2.  Minimum  roadway  length,  for  economization  of 
computer  time ,  and 

3.  Normalization  of  traffic  behavior  from  initial 
generation  state,  especially  representative  distribution  into 
the  second  lane  for  the  unidirectional  flow. 

Summaries  of  the  results  of  the  BRPOS  tests  are  shown  in 
Tables  1.  through  7.   Tables  1.  through  5.  reflect  the 
main  test  parameters,  for  each  roadway  approach  length,  and 
comparisons  of  prescribed  truck  incidence,  i.e.,  the  input 
distribution  function,  versus  the  truck  population  sample 
captured  on  the  bridge  deck.   These  tables  also  contain  the 
comparisons  of  the  given  platoon  distribution  function 
versus  the  platoon  population  generated  and  versus  the  platoon 
population  captured  on  the  bridge  deck.   The  platoon  distri- 
bution function  utilized  was  purely  theoretical  and  was 
a  binary  distribution  function.   No  data  on  the  incidence 
of  various  platoon  sizes  over  a  given  length  of  roadway 
was  found  to  provide  an  experimental  basis.   The  binary 
distribution  function  was  used  for  testing  purposes.   How- 
ever, a  theoretical  basis  for  its  use  can  be  developed 
to  support  interim  use  of  the  function. 

Table  6  provides  a  comparison  of  the  variations,  in  percent, 
of  the  truck  population  captured  on  the  deck  versus  the  given 
truck  population.   The  mean  variation,  in  percent,  of  each 
roadway  length,  extracted  from  Table  6,  is  as  follows: 
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Approach  Length  Mean  Variation 

200'  .12% 

1,000'  .12% 

4,000'  .13% 

10,560'  .14% 

26,400'  .20% 

while  the  mean  absolute  variation,  in  percent,  from  Table 
6,  is  .14%.   The  above  data  implies  that  the  two  shorter 
roadway  lengths  better  preserve  prescribed  platoon  distri- 
bution functions  than  the  longer  roadway  approach  lengths. 

The  mean  variation,  in  percent,  of  each  roadway  length, 
extracted  from  Table  7.,  is  as  follows: 

Approach  Length  Mean  Variation 

200'  3.9% 

1,000'  2.0% 

4,000'  3.25% 

10,560'  3.25% 

26,400'  5.45% 

while  the  mean  absolute  variation,  in  percent,  from  Table 
7.,  is  3.56%.   The  above  data  implies  that  the  1,000  foot 
roadway  approach  preserves  specified  platoon  distributions 
for  a  given  highway  better  than  the  other  tested  roadway 
approach  lengths.   The  next  two  best  roadway  lengths  were 
4,000  feet  and  10,560  feet. 
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As  a  result  of  the  aforementioned  tests,  it  is  obvious 
that  short  roadway  approach  lengths,  BRPOS,  are  better 
for  use  with  the  simulator.   Of  the  above 
lengths  listed,  1,000  feet  gave  the  best  results. 

The  reason  for  the  poor  showing  of  the  longer  approach 

lengths,  while  intuitively  they  would  appear  to  be  better 

candidates  for  providing  stable  traffic  behavior,  is  due 

to  the  degradation  of  the  simulation  algorithm  over  long 

roadways.   This  problem  has  been  discussed  elsewhere  in  this 

report.   Briefly,  the  algorithm  first  converges  to  the 

vehicle  steady  state  speed,  initially  generated  in  the 

pseudo-random  process  and  second  converges  on  a  zero  speed, 

without  a  lower  limit  imposed.   In  order  to  utilize  the 

program  a  value  of  BRPOS  must  be  selected  which  allows 

sufficient  simulation  time,  to  the  bridge,  for  dynamic 

traffic  behavior  to  propagate  and  for  a  distribution  into 

the  second  lane  to  occur,  as  a  lower  limit,  prior  to  any  de- 
gradation of  the  performance  of  the  vehicles. 

The  tests  described  herein,  then,  indicate  that  an  upper  limit 
on  BRPOS  exists  wherein  the  simulator  fails  to  simulate  properly 

Some  dependency  may  exist  on  the  given  highway  statistical 
data.   However,  this  should  not  be  a  strong  tendency  unless 
the  population  is  extremely  small.   In  general,  select 
BRPOS  as  follows: 


200  feet   <  (BRPOS  -  BRLEN)  <   4,000  feet 
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An. initial  value  of  1,000  feet,  plus  one-half  the  bridge 
length,  should  be  used  for  any  new  application  until 
evidence  indicates  that  a  better  value  exists. 

It  should  be  noted  that  this  tendency  of  the  simulator 
is  a  fortunate  one,  since  shorter  computer  execution 
times  are  implied  by  lower  values  of  BRPOS,  e.g.,  from 
Tables  1  through  5. 

Approach  Length         Computer  Execution  Time 

2001  32  sec. 

1,000'  37  sec. 

4,000'  69  sec. 

10,560'  137  sec. 

26,400'  315  sec. 

SIMTIM  Testing  -  Upon  completion  and  quick  analysis  of  the 
results   of  the.  BRPOS  tests,  a  value  of  1,000  feet  for 
roadway  external  to  the  bridge  was  selected  for  this  test. 
All  other  parameters  were  held  at  the  same  values  as  utilized 
in  the  BRPOS  tests. 

Prior  to  specifying  the  SIMTIM  test  an  estimate  of  the 
length  of  simulation  time  necessary  to  provide  a  90%  con- 
fidence   interval  that  at  least  one  occurrance  of  each 
type  vehicle  would  happen,  based  upon  the  given  truck  popu- 
lation distribution  function.   This  led  to  a  value  of 
approximately  five  hours  of  simulated  time.   Also,  the 
previous  test  had  used  one-half  hour  as  a  simulation  period 
and  achieved  reasonable  results.   Hence,  the  test  was  specified 
as  previously  defined  and  conducted  as  specified. 
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Similar  to  the  previous  test,  this  test  should  be  evaluated 
upon  the  basis  of  preservation  of  given  input,  traffic  statistics 
Tables  8  through  12  present  the  variations  in  percent, 
from  the  given  statistics  for  total  traffic  population  and 
for  distribution  of  platoon  size.   Theoretically,  the  longer 
the  simulation  tiiue,  1 . .& . ;    the  larger  the  sample,  the 
closer  the  generated  statistics  converge  upon  the  given 
statistics.   It  was  not  felt  necessary  to  compare  the 
generated  traffic  population  data  since  its  content  remained 
invariant.   However,  the  generated  platoon  size  data  was 
included  in  order  to  evaluate  maintenance  of  the  given  platoon 
statistics  at  the  point  of  generation.   Further,  it  was 
seen  that  the  traffic  simulation  superimposed  a  platoon 
stringing  effect  which  was  a  function  of  the  pseudo  bridge 
length.   This  was  evidenced  by  the  occurrence  of  platoon 
sizes  larger  than  those  generated.   It  was  due  to  the  definition 
of  a  bridge  platoon  used  in  the  simulator,  as  previously 
discussed.   As  the  pseudo  bridge  length  grew   longer, 
the  size  of  a  platoon  increased.   Conversely,  as  the  pseudo 
bridge  length  grew   shorter,  the  size  of  a  bridge  platoon 
decreased.  This  was  further  impacted  by  the  simulator  logic 
and  data  interface  requirements  to  the  stress  analysis 
program.   The  simulator  continued  a  truck  platoon  event 
beyond  the  actual  existence  of  a  truck  axle  within  the  pseudo 
bridge  length.   Any  truck  axle  that  would  appear  on  the 
pseudo  bridge  within  the  next  time  increment  acted  to  con- 
tinue the  event.   Hence,  at  maximum  acceptable  speed  the 
pseudo  bridge  length  was  artifically  extended  approximately 
100  feet.   In  the  test  cases  reported  on,  the  platoon 
defining  distance  became,  in  the  worst  case,  approximately 
300  feet.   Thus  the  platoon  statistics  determined  within 
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the  pseudo  bridge  length  reflected  the  existence  of  truck 
separation  distances,  i.e.,  between  rear  and  front  axles 
of  following  trucks,  ranging  from  approximately  2  0  feet  to 
approximately  300  feet  for  the  test  cases.   This  had  no 
material  effect  upon  the  prediction  of  the  stress  dynamics 
The  stress  analysis  program  defines  its  own  truck 
platoon  event  based  upon  a  steady-state  through  load- 
state  to  steady-state  behavior  of  all  points  evaluated 
on  the  bridge.   Thus  a  platoon  event  specified  by  the 
simulator  was  purposely  defined,  conservatively,  to 
contain  structural  dynamic  platoon  events  such  that  they 
were  not  arbitrarily  destroyed  or  modified  due  to  any 
simulator  mechanisms. 

The  above,  then,  tended  to  disqualify  the  bridge  captured 
platoon  size  statistics  as  a  valid  means  of  assessing  a 
"best"  SIMTIM  value.   This   left   the  overall,  gross,  popu- 
lation statistics  and  the  generated  platoon  statics  preser- 
vation on  the  bridge  to  be  evaluated  against  the  given  input 
statistical  data.   Table  13  reflects   a   comparison  of  the 
statistics  for  the  values  of  SIMTIM  used  in  the  testing. 
As  indicated  above,  the  longer  the  simulation  time, theoretically, 
the  better  the  representation  of  the  given  statistics  should 
become.   However,  a  rather  strange  behavior  of  the  simulator 
was  evidenced  in  the  data  shown  in  Tables  8  through  13. 
The  theoretically  anticipated  improvement  was  evidenced  in  the 
decrease    of  variation  from  the  given  distributions 
from  a  one  hour  value  through  the  five  hour  value,  but  the 
ten  hour  value  suddenly  provided  a  poorer  representation  of 
the  given  distributions.   This  behavior  was  evidenced  in 
both  the  generated  traffic  population  distribution  and  the 
generated  platoon  size  distribution  data.   The  cause  of  this 
degradation  at  a  ten  hour  simulated  time  was  not  investi- 
gated in  depth  due  to  the  constraint  imposed  by  other  work 
remaining  to  be  accomplished  and  the  limitations  of  funds 

and  time. 
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This  could  represent  another  very  serious  handicap 
in  the  use  of  the  BRIGLDl  simulator  as  a  synthetic 
bridge  load  generator.   Since  future  use 

of  this  generator  was  predicated  on  taking  time  slice  samples 
of  approximately  2  weeks,  336  hours  of  simulated  time,  as 
a  representative  sample  of  a  longer  period,  e.g.,  ten 
years,  and  repetition  of  this  for  a  longer  total  bridge 
life  span,  e.g.,  fifty  years,  the  above  encountered  problem 
bears  significant  investigation  prior  to  such  use. 

The  authors  believe  that  the  basis  of  this  problem  arose 
from  the  misuse  of  the  pseudo-random  number  generating 
subroutine.   It  is  used  up  to  a  maximum  of  six  calls  in  the 
generation  of  a  truck  event,  i.e.,  to  propagate  up  to  six 
separate  and  independent  random  functions.   In  the  case  of 
the  ten  hour  test,  this  implied  the  generation  of  a  maximum 
of  110,760  pseudo  random  numbers.   This  in  itself  is  not 
a  serious  problem.   However,  all  of  these  numbers,  for  six 
separate  and  independent  functions ,  are  generated  from  the 
same  series  and  initiated  by  a  single  seed.   It  is  strongly 
recommended,  prior  to  extended  use  of  this  load  generator, 
that  six  separate  seeds  be  utilized,  and  maintained,  one 
for  each  function.   Further,  each  of  the  six  series  should 
be  periodically  re-initiated  with  new  seeds. 

i 
It  should  be  noted  that  the  original  version  of  this  pro- 
gram did  not  provide  for  user  initiation  of  even  the  single 
seed.   This  caused  the  program  to  always  generate  the  same 
vehicle  events  in  the  same  order. 

In  terms  of  preserving  the  given  traffic  population  statistics 
and  the  platoon  size  statistics  the  five  hour  simulation  time 
provided  the  best  results.   However,  a  five  hour  run  is  not 
an  adequate  sample  for  long  term  purposes. 
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DELTIM  Testing  -  Since  the  variation  of  DELTIM  (At)  re- 
lates primarily  to  the  integration  of  the  differential  equations 
of  motion,  it  could  only  be  evaluated  in  terms  of  the  accuracy 
of  vehicle  motion  synthesis.   However,  in  this  particular 
problem,  i.e.,  bridge  loading,  the  motion  was  of  lesser 
consequence  than  the  representation  of  the  loadings  on  a 
bridge  deck.   Since  no  control  data  exists,  i.e.,  field 
data  on  how  any  given  bridge  deck  is  loaded  over  any  signifi- 
cant period  of  time ,  there  was  no  direct  comparative  test 
available . 

Further,  analyzing  the  motion  behavior  of  a  large  set 
of  vehicles  in  a  detailed  accuracy  evaluation  of  motion 
behavior  is  beyond  the  scope  of  this  effort.   A  compromise 
evaluation  approach  for  this  parametric  test  was  utilized. 
For  the  four  cases  run  a  comparison  at  simulation  times 
common  to  all  cases  was  made.   The  comparison  was  in  terms  of 

1.  Commonality  of  deck  loading  vehicle  set 

2.  Sequence  retention  of  the  vehicles  forming  the 
set 

3.  Relative  position  of  the  specific  vehicles  in 
the  set. 

This  analysis  was  then,  essentially,  based  upon  determining 

1.  Similarity  of  the  truck  load  events  generated 

2.  Occurrence  of  the  events  at  the  same  time. 
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It  was  necessary  to  assume  the  At  =  .5  sec.  as  a  reference, 
or  control  case.   This  value  of  At  was  small  enough  to 
prevent  the  overshooting  of  any  logical  operation  per- 
formed by  the  simulator  in  synthesizing  traffic  behavior. 
Hence/  all  represented  traffic  behavior  functions,  within 
the  simulator,  should  have  occurred  properly  with  this  value 

The  test  parameters,  and  the  test  cases,  utilized  were 
as  follows: 

1.  Traffic  distribution  data  for  Md.  301 

2.  A  binary  platoon  distribution  at  generation 

3.  SIMTIM  =5.0  hours 

4.  Print-out  of  deck  loads  =  1.0  hour 

5.  Bridge  length  =  250  ft. 

6.  Roadway  length  =  1000  ft. 

7.  At  variations  were: 

Case  #1  =  .5  sec. 

Case  #2  =  1.0  sec. 

Case  #3  =  1.5  sec. 

Case  #4  =  2.0  sec. 

8.  Data  evaluation  limited  to  1/2  hour  of  bridge 
loading  data  at  comparison  values  of  time. 
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The  reason  for  limiting  the  detailed  bridge  load  data  out- 
put to  one  hour  was  due  to  the  extremely  large  quantity 
of  output  and  the  difficulties  in  arranging  for  the  generation 
of  such  a  large  amount  of  printing.   System  constraints 
further  created  difficulties  in  the  generation  of  these 
large  quantities  of  output. 

Basic  data  generated  by  the  test  cases  is  shown  in  Table 
14.   It  should  be  noted  that  the  characteristic  of  stringing 
smaller  platoon  sizes  together  to  form  a  larger  platoon 
also  was  true  of  larger  integrating  intervals,  as  is  in- 
dicated in  Table  14.   This  again  was  directly  due  to  the 
definition  of  a  platoon  event  utilized  within  the  simulator, 
as  discussed  previously. 

The  analysis  of  the  results  are  ordered  in  the  following 
manner: 

1.  DELTIM  -  1.5  sec.  vs.  .5  sec. 

2.  DELTIM  -  2.0  sec.  vs.  .5  sec. 

3.  DELTIM  =  1.0  sec.  vs.  .5  sec. 

The  above  order  was  due  to  the  order  in  which  the  test 
case  results  were  produced,  from  the  computer,  and  for  no 
other   reason. 

Calculationally  there  is  no  difference,  i.e.,  at  most 
+  .01  feet,  in  the  use  of  At  =  1.5  versus  0.5  where  the 
event  is  exactly  the  same  on  the  bridge  deck.   This  implies 
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that  the  simulated  vehicles  have  undergone  exactly  the  same 
logical  process  in  both  cases.   It  further  implies  that 
the  numerical  integration  error  propagated  by  the  choice 
of  At  was  satisfactory  for  1.5  seconds,  for  the  dynamics 
involved. 

However,  there  may  be  significant  differences  in  the  structure 
of  a  load  event,  at  any  instant  of  time  on  a   bridge  deck. 
An  analysis  of  the  At  =  1.5  case  against  the  At  =  0.5  case 
indicated 

1.   A  significant  difference  of  the  total  number 
of  load  events  defined,  i.e., 

a.  140  events  for  At  =  0.5  sec. 

b.  120  events  for  At  =  1.5  sec. 

This,  in  itself,  could  cause  no  significant  affect  upon  the 
structural  program,  if  the  spatial  relationships  of  the  vehicles 
defining  each  event  were  preserved.   This  is  true  because 
the  structural  program  disaggregates  the  simulator  defined 
events  into  a  new  set  of  load  events  that  are  defined  by 
the  dynamic  behavior  of  any  given  bridge,  i.e.,  from  steady 
state  through  the  transient  load  state  and  back  to  steady 
state . 

However,  the  observed  behavior  of  events,  see  Table  15., 
between  the  At  =  0.5  and  the  At  =  1.5  cases,  indicated the 
following  behavior: 
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1.  Exact  duplication  of  an  event. 

2.  Minor  spatial  variations  but  retaining  the 
same  event  content  and  relationships. 

3.  Large  spatial  variations,  between  10  feet  and  100 
feet  positionally ,  but  retaining  the  same  event  content 

and  relationships . 

4 .  Large  spatial  variations  and  variation  in 
relationships  but  retaining  the  same  content. 

5.  Some  vehicles  retaining  same  positions  and 
relationships,  but  total  vehicle  content  not  the  same. 

6.  Same  or  similar  event,  but  an  unexplainable  time 
of  occurrence  variation. 

At  this  point  in  the  research,  the  above  variations  are  not 
necessarily  significant.   The  economies  saved  by  using  a 
larger  At  may  be  worth  far  more  than  attempting  to  preserve 
unknown  or  doubtful  event  statistics.   The  variations  in 
the  definitions  of  the  events  due  to  choice  of  At  were  due 
to  the  reasons  discussed  below. 

The  cause  of  the  variation  in  the  number  of  events  as  a 
function  of  At  was  due  to  the  effective  bridge  length  used 

by  the  program  to  define  an  event,  i.e.,  if  a  truck  would 

have  an  axle  on  the  bridge  at  the  next  time  step  and  a  truck 

still  had  an  axle  on  the  deck  in   a  current  time  step, 

then  the  event  continued  as  the  same  event.  The  implied 
bridge  lengths  as  a  function  of  At  were: 
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At  Max .  Speed        Minimum  Speed 

.  .5  279  ft.  301  ft. 

1.0  308  ft.  353  ft. 

1.5  337  ft.  404  ft. 

2.0  366  ft.  456  ft. 

where  the  specified  bridge  length  was  250  ft. 

When  the  logical  occurrences  were  identical  a  At  of  1.5  seconds 
gave   exactly  the  same  positional  values, at  the  same  time, 
as  At  of . 5  seconds. 

Automobile  motion  appeared  to  be  logically  manipulated  better 
at  At  =  .5  than  at  At  =  1.5  which  also  causes  a  variation 
in  event  structure.  Typically  more  autos  were  present 
in  the  At  =  .5  samples  than  the  At  =  1.5  and  tended  to  be 
more  intermingled  with  the  trucks  than  the  At  =  1.5  case. 
This  was  valid  if  the  manner  in  which  the  passing  logic 
operates  is  considered. 

At  maximum  allowed  speed,  an  auto  could  move~117  feet  for  a 
At  of  1.0  seconds.   However,  it  needed  to  only  move  60  feet 
to  satisfy  the  minimum  passing  criteria,  three  logical 
steps  were  involved, 

1.  Move  from  lane  1  to  lane  2 

2.  Pass  automobile  in  lane  1 

3 .  Move  from  lane  2  to  lane  1 . 

In  order  to  accomplish  these  three  steps  a  minimum  of  two 
time  steps  were  required.  To  pass  a  maximum  length  truck  an 
automobile  needed  to  only  move  9  4  feet  in  one  time  step. 
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A  comparison  of  the  At  =  2.0  seconds  with  the  At  =  0.5 
seconds  indicated  the  same  numerical  error.   The  incidence 
of  agreement  in  event  content,  spatial  accuracy  and  re- 
lationships was  obviously  less  than  for  At  =  1.5,  see  Table 
16.,  8.8%  for  At  =  2.0  sec.  vs.  16%  forAt  =  1.5  for  exact 
agreement.   It  should  be  noted  that  these  figures  did  not 
represent  all  time  points,  only  those  points  sampled. 
The  use  of  a  At  =  2.0  seconds  appeared  to  induce  very 
serious  degradation  in  the  definition  of  deck  load  events. 
However,  calculational  accuracy  of  the  numerical  integration 
was  retained  for  those  events  having  exactly  the  same  logical 
operations  imposed  upon  them. 

Again,  the  use  of  At  =  2.0  seconds  is  of  consequence  only 
if  valid  event  statistics  are  to  be  preserved. 

The  comparison  of  At  =  1.0  sec.  case  with  the  At  =  0.5 
sec.  case  indicated  a  much  higher  correlation  of  the  total 
number  of  event  samples,  see  Table  17,  than  the  At  =  1.5 
sec.  case.   The  At  =  1.5  sec.  case  had  an  incidence  of  21% 
of  the  event  samples  varying  from  the  At  =  . 5  case,  while 
the  At  =  1.0  sec.  case  had  an  incidence  of  8.1%.   The 
At  =  1.0  sec.  case  in  general  provided  a  much  better  com- 
parison to  the  At  =  0.5  case  than  either  the  At  =  1.5 
or  At  =  2.0  sec.  cases. 

Analytically,  At  =  . 5  seconds  ensures  that  passing  logic 
will  be  executed  properly,  since  it  prevents  an  overshoot 
during  one  time  step.   However,  unless  valid  statistics 
are  provided  as  parameters  to  this  program  the  events 
generated  on  the  deck  for  At  =  1.0  seconds  are  quite  as 
likely  as  those  generated  at  At  =  0.5  seconds.   Hence,  for 
the  sake  of  economy  and  a  compromise  with  the  complete 
preservation  of  event  behavior  a  At  =  1.0  second  is  recommended, 
at  this  time. 
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If  valid  load  event  statistics  are  provided,  via  the  para- 
metric input  to  the  program,  a  At  =  .5  seconds  should  be  used 

A  comparison  of  the  three  test  case  variations  is  shown 
in  Table  18.   Significant  factors  were 

1.  As  the  time  increment  increases  fewer  comparison 
event  points  existed 

2 .  Existence  of  the  exact  event  in  position  and  time 
on  the  pseudo-bridge  showed  a  non-decisive  behavior. 

3.  The  existence  of  the  same  event  structure  and 
content,  with  minor  positional  variances  of  a  few  feet, were 
about  the  same  for  increments  of  1.0  and  1.5  seconds. 
However,  these  provided  twice  the  similarity  as  does  an 
increment  of  2.0  seconds . 

4.  Similar  event  structures  and  content  with  large 
positional  variances  showed  a  non-decisive  behavior. 

5.  There  appears  to  be  about  a  10%  existence,  for 
all  cases,  of  the  same  event  content  but  different  spatial 
relationships  between  the  vehicles,  i.e.,  leading  vs. 
following,  and  lane  occupancy. 

6 .  There  was  clearly  a  decrease  in  the  incidence  of 
dissimilar  events  at  the  same  sample  time  points  as  the 
time  increment  approaches  .5  seconds. 

7.  There  was  also  an  apparent  decrease  in  the  occur- 
rence of  unexplainable  time  variations  for  the  same  or 
similar  events,  as  the  time  increment  approached  .5  seconds. 
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There  also  appeared  to  be  a  divergence  trend,  as  simulation 
time  increases,  of  event  equality  and  similarity  between 
the  different  time  increment  cases.   The  larger  time  in- 
crement cases  diverged  more  rapidly  than  the  smaller  time 
increment  cases.   There  was  also  a  definite  decrease  in  the 
number  of  unexplainable  deviations  in  time  as  the  DELTIM 
value  approaches  .5  seconds. 

In  general,  the  results  shown  in  Table  18.,  as  should  be 
expected,  indicate  better  correspondence  with  the  .5  second 
case  as  the  time  increment  approaches  that  value.   As  pre- 
viously indicated,  a  value  of  1.0  seconds  is  borderline 
and  usable  unless  specific  platoon  event  statistics  are 
to  be  preserved.   In  this  instance,  a  value  of  .5  seconds 
should  be  used.   For  the  purposes  of  economy  the  larger 
time  increments  are  better,  obviously. 

Restricted  Zone  Test  -  The  testing  of  the  restricted  zone 
capability  of  the  simulator  was  more  of  an  analysis  of  its 
proper  functioning  than  a  sensitivity  analysis.   However, 
a  sensitivity  analysis  was  performed  on  the  resulting  data, 
in  terms  of  affect  upon  generated  truck  load  events.   A 
control  case  having  all  of  the  same  parameters,  excepting 
restricted  zones,  was  run  to  provide  a  basis  for  evaluating 
the  effects  to  the  generated  bridge  loads  by  various  zoning 
options  and  configurations,  as  previously  described  in  the 
test  specification.   Any  other  form  of  analysis  would  be 
meaningless  in  the  context  of  bridge  loading  synthesis. 

The  basic  parameters  were  as  follows: 

1.  Md.  301  traffic  distribution  data. 

2.  A  binary  platoon  distribution  function. 
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3.  SIMTIM  =  1  hour. 

4.  SDFAC  =  15. 

5.  An  approach  roadway  of  2000  feet. 

6 .  t  =  1 . 0  sec. 

Traffic  flow  was  basically  unidirectional  for  the  sensitivity 
testing  with  functional  testing  of  bidirectional  flow. 

The  analysis  of  the  effects  was  performed  in  the  same 
manner  as  for  the  DELTIM  testing,  i.e.,  in  terms  of  the 
similarity  of  the  load  events  on  the  pseudo-bridge  at  equiva- 
lent simulation  times. 

The  effects  reflected  in  the  unidirectional  flow  cases 
were  first  evaluated,  in  the  following  order: 

Unidirectional 

1.  Case  f  (Control  Case)  vs.  Case  2a  (50  ft.  re- 
stricted with  no  grades). 

2.  Case  f  (Control  Case)  vs.  Case  2b  (750  ft.  re- 
stricted with  10%  grade) . 

3.  Case  f  (Control  Case)  vs.  Case  2c  (1950'  re- 
stricted with  -5%  grade) . 

4.  Case  f  (Control  Case)  vs.  Case  2d  (750*  restricted 
with  10%  grade) . 

5.  Case  f  (Control  Case)  vs.  Case  2e  (1950  *  restricted 
with  5%  grade) . 
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6.  Case  f  (Control  Case)  vs.  Case  la  (5  zone  hilly 
type  approach) . 

Bidirectional 

7.  Case  lb  (i)  (Control  Case)  vs.  Case  lb  (iii) 
(50  ft.  restricted  with  no  grades) . 

8.  Case  lb(i)  (Control  Case)  vs.  Case  lb(ii) 
(5  zone  hilly  type  approach) . 

The  analysis  of  Case  2a,  the  50'  restricted  zone  with 
unidirectional  flow,  adjacent  to  the  bridge  without  grade 
indicated  that  variations  were  induced  when  compared  to  the 
control  case.   The  choice  of  a  50'  zone  length  was  almost 
too  short  to  provide  a  measurable  variation.   However,  the 
following  characteristics  were  observed,  due  to  this  parameter 
variation,  in  a  sample  of  80  events  in  a  simulation  time 
of  1/2  hour: 

1.  No  change  in  event  occurrence  was  observed.   A 
one-to-one  correspondence  existed. 

2.  All  events  occurred  at  exactly  the  same  time 
as  for  the  control  case. 

3.  Sparce  vehicle  events  gave  no  evidence  of  variation 
Variations  tended  to  occur  for  dense  vehicle  events. 

4.  Exact  duplication  of  events,  and  their  time  of 
occurrence,  occurred  in  76.25%  of  the  events. 

5.  Minor  positional  variations  were  observed  in 
17.5%  of  the  events.   These  variations  were  typically  due 

to  auto  position  variations,  as  opposed  to  truck  variations. 
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6.  Large  positional  variations  occurred  in  3.75%' 
of  the  events ,  and  were  typically  again  observed  in  auto 
position  change. 

7.  Only  2.5%  of  the  events  evidenced  vehicle  con- 
tent variation.   In  both  cases,  the  variation  in  the  vehicle 
content  was  one  auto. 

The  use  of  a  restricted  zone,  without  grade,  as  short  as 
50  feet  will  produce  variations  in  the  behavior  of  truck 
load  events  for  a  pseudo-bridge  length.   However,  such 
variations,  based  on  this  test  do  not  materially  affect 
the  vehicle  content  nor  the  timing  of  a  load  event. 

It  should  also  be  noted  that  the  control  case  utilized  1 
minute  and  29  seconds  of  computer  time  in  the  "go"  step, 
while  this  case  utilized  1  minute  and  32  seconds  in  the  go 
step.,  an  increase  of  3  seconds. 

The  analysis  of  case  2b,  the  750'  -10%  down  grade,  test 
case,  as  described  in  the  specification;  showed  absolutely 
no  effect  of  the  imposed  down  grade  restricted  zone  for  uni- 
directional traffic.   The  analysis  was  performed  on  an  80 
load  event  sample  over  1/2  hour  of  simulated  time  by  com- 
parison with  the  control  case,  2f.   The  generated  loads 
were  identical  in  space  and  time  to  those  generated  by  the 
control  case.   This  parameter  variation  capability  within 
the  simulator  either  has  no  meaningful  effect  because  of 
its  logical  definition  to  the  form  of  use  applied  in  this 
test  case,  or  it  has  never  worked.   It  is  the  opinion  of  the 
investigator  that  this  use  of  the  BRIGLD1  simulator  had  not 
been  either  utilized  or  debugged,  previously.   No  additional 
effort  was  expended  on  this  apparent  problem. 
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It  should  be  noted  that  all  of  the  restricted  zone  cases 
required  a  longer  approach  roadway  length,  approximately 
2000  feet,  and  because  of  this  the  accelerations  had  largely 
converged  to  zero  and  the  speeds  had  largely  converged  to  the 
minimum  allowed  speed.   CAse  2b  utilized  exactly  the  same 
amount  of  computer  time,  in  the  "go"step,  as  did  the  control 
case,  i.  e.,  1  minute  and  29  seconds.   This  further  reinforced 
the  above  conclusion  that  this  variation  of  the  restricted 
zone  capability  had  no  affect  whatsoever  on   the  behavior 
of  the  synthesized  traffic. 

The  analysis  of  case  2c,  the  1950',  -5%  downgrade,  test 
case  resulted,  again,  in  determining  that  the  generated 
load  data  was  exactly  the  same  as  the  control  case,  case 
2f.   No  effect  whatever  from  this  parameter  variation  was 
reflected-  in  the  output  for  1/2  hour  of  simulated  time  spanning 
80  load  events.   The  conclusions  drawn  on  this  case  were 
exactly  the  same  as  those  presented  above  on  case  2b. 

One  interesting  variation  in  the  case  was  found.   While  the 
exact  replica  of  output  was  generated  by  this  case  for  the 
control  case,  the  computer  time  utilized  in  the  "go"  step 
for  this  case  was  1  minute  and  21  seconds  as  compared  to  1 
minute  and  29  seconds  for  the  control  case,  a  decrease  of 
8  seconds  over  the  control  case  time. 

In  general  ,  this  case  further  reinforced  the  belief  that 
the  restricted  zone  downgrade  capability  is  not  functioning 
in  any  manner  for  unidirectional  traffic  flow,  i.  e.,   it 
had  no  meaning  in  the  simulator. 
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The  results  of  case  2d,  the  750 '  10%  upgrade  case,  were 
heartening  after  the  two  previous  test  cases.   When  compared 
to  the  control  case,  the  following  was  observed. 

1.  The  first  event  for  this  case  lagged  the  control 
case's  first  event  by  2  seconds.   This  implies  the  existence 
of  slower  moving  traffic  than  for  the  control  case. 

2.  This  case,  while  generating  exactly  the  same 
set  of  vehicles  as  the  control  case,  compressed  the  number 
of  separate  load  events,  as  previously  defined,  into  135 

as  compared  to  151  for  the  control  case;  implying  a  stringing 
affect  due  to  slower  vehicle  action. 

3.  Significantly  slower  speeds  for  truck  traffic 
were  noted  at  the  point  the  trucks  entered  the  pseudo- 
bridge  . 

4 .  Longer  times  on  the  pseudo-bridge  were  also 
noted  for  the  turck  traffic,  which  was  completely  consistent 
with  the  other  noted  behavior. 

5.  The  structure  and  content  of  the  generated  load 
events  on  the  pseudo-bridge  were  radically  changed,  almost 
totally,  from  those  generated  by  the  control  case. 

6.  A  few  instantaneous  points,  especially  at  the 
commencement  of  a  load  event,  were  found  where  the  load 
event  was  exactly  the  same  as  for  the  control  case. 

7.  The  same  strong  tendency  for  the  accelerations  to 
converge  to  zero  and  the  speeds  to  converge  to  the  minimum 
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allowed  speed  was  noted.   However,  velocities  and  acceler- 
ations, where  different  from  the  minimum  speed  and  non- 
zero respectively,  were  almost  totally  different  from  those 
output  by  the  control  case.   Convergence  occurred  more 
rapidly  than  for  the  control  case  due  to  the  grade  slowing 
effect. 

In  general,  conclusions  drawn  from  this  test  indicated  that 
the  upgrade  effect,  at  least  for  a  significant  distance, 
did  significantly  alter  the  generated  truck  load  events 
from  a  non-grade  form.   Further,  the  simulation  did  appear 
to  reflect  the  physics  and  dynamics  expected  from  the  im- 
position of  a  significant  grade  in  the  approach  to  a  bridge. 
However,  no  indication  of  the  restricted  portion  of  the 
restricted  zone  capability  could  be  determined  separate 
from  the  effects  of  the  upgrade.   This  could  have  been  due 
to  the  fact  that  the  effects  of  both  variations  should  cause 
similar  effects  in  the  behavior  of  the  traffic.   However,  it 
is  the  opinion  of  the  investigators  that  the  restricted 
operation  did  not  function  and  that  all  of  the  effects  ob- 
served for  this  case  were  due  to  the  upgrade  variation.   This 
opinion  was  based  primarily  upon  the  results  from  the  two 
previous  test  cases. 

For  the  same  period  of  simulation  time,  this  case  utilized 
1  minute  and  27  seconds  of  computer  time  in  the  "go"  step 
as  compared  to  1  minute  and  29  seconds  for  the  control  case. 

The  results  of  case  2e ,  the  1950'  5%  upgrade  case,  were 
very  similar  to  those  of  the  previous  case,  case  2d,  in 
behavior.  In  general,  all  of  the  observed  effects  from  the 
previous  case  held  for  this  case.  Although,  each  event 
was  essentially  different  from  the  previous  case.  Some 
specific  results  were : 
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1.  The  first  event  for  this  case  lagged  the  control 
case's  first  event  by  3  seconds.   This  implies  slower 
moving  traffic  behavior  than  the  control  case. 

2.  Again,  the  load  events  were  largely  of  different 
structure  than  the  control  case  and  the  previous  case. 
Truck  behavior  was  significantly  different  than  the  previous 
case.   In  some  instances,  the  same  truck  in  this  case 
lagged  in  the  previous  case.   It  was  also  true  that  this 
case  lagged  the  previous  case  time-wise  for  some  trucks. 

3.  This  case  compressed  the  total  number  of  load 
events  to  132  as  compared  to  131  for  the  previous  case  and 
151  for  the  control  case.   In  general,  this  implied  an 
overall  slightly  faster  traffic  behavior  than  the  previous 
case,  which  was  consistent  with  the  slopes  of  the  grades 
imposed . 

4.  The  dynamics  of  the  vehicles  forming  the  load 
events  while  on  the  pseudo-bridge  were  much  different  than 
for  either  the  control  case  or  the  previous  case . 

The  overall  behavior  of  this  case  as  compared  to  the  con- 
trol case  was  consistent  with  the  anticipated  behavior 
expected  from  the  physics  of  the  problem.   Again,  no  measur- 
able effect   due  to  the  restricted  operation  could  be  dis- 
cerned from  the  effects  due  to  the  grade.   The  same  con- 
clusions were  drawn  with  respect  to  this  case  as  for  the 
previous  case. 
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Comparing  a  sampling  of  identical  trucks  between  this  case 
and  the  previous  case  indicated  an  arrival  time  on  the 
pseudo-bridge  for  this  case  of  approximately  3  seconds 
before  the  previous  case,  for  the  same  trucks.   In  a  very 
few  cases,  the  arrival  time  was  behind  the  previous  case, 
approximately  10%  of  the  sample.   The  speeds  of  the  identical 
trucks  on  arrival  at  the  pseudo-bridge  tended  to  be  slower 
or  the  same,  in  equal  proportions,  for  this  case  when 
compared  to  the  previous  case.   In  a  small  portion  of  the 
sample,  identical  trucks  arrived  at  the  bridge  with  faster 
speeds  than  in  the  previous  case.   These  results  appeared 
to  be  consistent  with  the  relative  physics  of  the  two 
cases,  i.e.,  this  case  had  a  grade  that  was  2.6  times  the 
length  of  the  grade  in  the  previous  case  although  with  only 
50%  of  the  grade. 

In  analyzing  these  two  cases,  another  deficiency  in  the 
simulator  became  obvious.   Truck  power  is  limited  to  a  single 
fixed,  input,  weight  to  horsepower  ratio  for  each  truck 
type.   No  provision  was  made  for  weight  variations.   A 
more  realistic  approach  would  be  to  input  a  nominal  horse- 
power value  for  each  type  and  calculate  the  ratio  as  a 
function  of  the  truck's  load.   As  it  is,  a  fully  loaded 
truck  of  a  given  type  has  the  same  weight  to  power  ratio, 
i.e.,  the  same  acceleration  potential,  as  an  empty  truck 
of  the  same  type.   This  poses  a  serious  question  about 
utilizing  the  load  data  generated  by  the  BRIGLD1  simulator 
for  structural  analysis  purposes.   Work  at  Le  High  has 
related  velocity-weight  relationships  with  bridge  structural 
response.   The  simulator  not  only  does  not  preserve  such 
relationships,  but  falsely  relates  the  two  variables  to  a 
probable  50%  of  the  trucks  generated. 
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In  analyzing  case  la,  the  hilly-like  approach  to  the  pseudo- 
bridge,  in  comparison  with  case  2f,  the  unidirectional  con- 
trol case,  significant  variation  in  the  structure  of  load 
events  and  the  dynamics  were  observed.   However,  the  effects 
were  much  less  in  magnitude  than  for  the  simple  grade  cases 
2d  and  2e. 

Observed  effects  were  as  follows: 

1.  The  first  event  occurred  at  exactly  the  same 

time  as  for  the  control  case,   but  was  different  in  structure 

2.  The  total  number  of  load  events  was  compressed 
to  145  events  as  compared  to  151  for  the  control  case,  and 
131  and  132  for  cases  2d  and  2e  respectively.   This  further 
substantiated  the  conclusion  that  the  hilly-like  variation 
provided  effects  in  lesser  magnitude  than  the  simpler 
upgrade  cases. 

3 .  Truck  behavior  and  ordering  was  exact  in  many 
instances,  with  a  definitely  measurable  variation  in  a 
significant  portion  of  the  truck  population. 

4.  Event  vehicle  content  variation  was  largely 
due  to  auto  variations. 

5.  Definite  lag  effects  were  observed  which  are 
consistent  with  other  observed  effects  from  this  case. 

6 .  The  initiation  of  many  events  were  the  same  or 
nearly  the  same  in  most  cases.   Variations  initially  or 
subsequently  were  usually  due  to  change  in  auto  content. 
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The  overall  observed  effects  were  qualitatively  what  should 
be  anticipated  in  terms  of  comparison  with  the  control  case. 
They  were  also  consistent  with  the  two  simple  upgrade  cases. 
While  the  grades  imposed  were  severe,  their  duration  was 
extremely  short,  with  a  total  upgrade  length  of  only  500 
feet.   Again,  no  direct  conclusions  can  be  drawn  about  the 
restricted  operation  for  the  same  reasons  as  previously 
discussed.   Similarly,  the  downgrade  effect  cannot  be 
evlauated,  but  it  is  assumed  to  be  non-existent  based 
upon  the  two  simple  downgrade  tests. 

The  results  of  case  9b(iii)  for  bidirectional  traffic 

flow  with  a  restricted  50  ft.  zone  at  one  approach  to  the 

pseudo-bridge  indicates  very  little  variation  from  the  control 

bidirectional  flow  case,  case  lb(i).   On  the  order  of 

ten  variations  were  observed  over  a  1400  second  simulation 

time  period. 

Both  cases  timed  out  on  the  computer  at  5  minutes  of  exe- 
cution.  The  control  case  utilized  4  minutes  and  40  seconds 
of  computer  time  in  the  "go"  step  for  1672  seconds  of 
simulated  time,  while  this  case  used  exactly  the  same  time, 
but  for  only  1401  seconds  of  simulation  time. 

Variations  were  in  terms  of  auto  positions  at  the  same 
sample  time.   Truck  behavior  was  exactly  the  same  as  for 
the  control  case.   No  change  in  the  number  of  events  over 
the  same  simulation  period  occurred. 

The  zone  restriction  operator  appeared  to  function  for  bi- 
directional traffic  flow,  as  opposed  to  the  unidirectional 
case . 
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Additionally,  initialization  of  truck  velocities,  in  the 
2nd  direction,  by  the  simulator's  generator  evidenced  a 
tendency  to  initialize  improper  velocities. 

The  results  of  the  5  restricted  zones  bidirectional  test 
case,  case  lb(ii),  when  compared  to  the  bidirectional 
control  case,  case  lb(i),  indicated  minor  but  definite 
variations  from  the  control  case.   This  case  represents 
a  hilly  section  of  road  adjacent  to  one  end  of  the 
pseudo-bridge.   Ordering  of  truck  load  events  appear 
consistent  with  the  control  case.   Observed  variations  were 
in  terms  of  change  of  position  at  the  same  simulation  times 

A  total  of  1401  seconds  of  simulation  time  was  synthesized 
in  a  total  of  4  minutes  and  40  seconds  of  computer  time 
in  the  "go"  step.   This  compared  very  closely  to  the  pre- 
viously described  case,  case  lb(iii).   It  condensed  the 
total  number  of  cases  to  53  as  compared  to  54  for  the  con- 
trol case. 

A  total  of  22  events  varied  from  the  control  case.   Most 
variations  were  due  to  positional  differences  for  autos . 
In  a  few  cases,  a  change  in  the  number  of  autos  contained 
in  an  event  occurred  and  in  a  few  others,  five,  truck 
positions  changed. 

This  case  evidenced  the  same  form  of  variation,  slight  but 
definite,  to  its  control  case,  as  did  the  similar  uni- 
directional case.   In  this  instance,  the  restricted  zoning 
largely  only  affected  50%  of  the  generated  vehicles. 

Its  relative  effects  appeared  to  relate  to  the  physics  of 
the  case. 
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Bidirectional  traffic  flow  simulation  increased  computer 
CP.  time  by  a  factor  of  7  over  unidirectional  flow  simulation, 
using  5  stipulated  zones  for  both  cases,  i.e.,  exactly  the 
same. 

Bidirectional  traffic  flow  with  5  zones  ran   at  an  increase 
of  .03  CP  seconds  per  second  of  CP  time  for  a  no  zone 
bidirectional  case.   This  implied  that  the  bidirectional 
simulation  is  a  very  costly  process. 

A  comparison  of  1  zone  stipulated  to  5  makes  no  significant 
change  in  running  time  (CP  time) . 

In  addition  to  the  excessive  computer  utilization  incurred 
by  this  option,  and  the  incorrect  velocity  generation  of 
trucks  in  the  2nd  direction,  there  is  a  serious  question 
about  the  validity  of  the  generated  density  for  a  two  di- 
rectional highway.   It  is  essential  that  a  nominal  one  lane 
traffic  load  be  defined  for  a  bidirectional  highway  be- 
cause the  simulator  apparently  doubles  the  given  vehicle 
density. 

Use  of  this  option  is  not  recommended  unless  significant 
improvement  on  the  two  way  synthesis  is  performed.   Further, 
the  excessive  computer  costs  make  its  use  questionable  for 
the  purpose  it  is  supposed  to  serve,  i.e.,  long  term  load 
prediction. 

Conclusions  drawn  on  the  use  of  the  restricted  zone  options 
are  summarized  as  follows: 

1.   Restricted  flow  of  unidirectional  flow  did 
not  appear  to  have  any  meaning  in  the  simulator  and  only 
qualified  effect  for  bidirectional  flow. 
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2.  Downgrade  definitions  did  not  appear  to  have 
any  meaning  in  the  simulator. 

3.  Upgrade  definitions  appeared  to  operate  properly. 

4.  As  a  result  of  analyzing  these  test  cases,  test 
results  indicate  that  the  bidirectional  flow  simulation  is 
questionable  and  certainly  used  a  large  amount  of  computer 
time,  as  compared  to  unidirectional  flow. 

5.  Weight-velocity  relationships  were  not  adequately 
preserved. 

SDFAC  Testing  -  This  variable  controls  the  value  utilized 
for  safe  following  distance  during  a  simulation.   A  value  of 
15,  as  indicated  earlier  is  approximately  the  rule  of  thumb 
of  one  car  length  per  each  ten  miles  of  speed,  e.g.,  at 
60  mph  it  establishes  117  feet  as  a  safe  following  distance. 

In  evaluating  these  test  cases,  the  SDFAC=15  case  was  utilized 
as  a  control  case.   In  all  cases,  a  simulation  time  of  one 
hour  was  utilized,  with  detailed  loading  output  at  each 
time  point  within  each  event.   The  use  of  one  hour  was 
based  upon  the  pragmatics  of  analyzing  a  large  quantity 
of  output  and  the  actual  printing  of  output.   One  hour  of 
simulation  generated  approximately  1800  vehicles  for  the  traffic 
distributions  utilized.   The  control  case  generated  131 
separate  load  events  for  this  period  of  simulation. 

The  other  important  parameters  utilized  were: 

1.   An  approach  length  of  1000  feet 
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2.  A  bridge  length  of  250  feet 

3.  Unidirectional  traffic  flow 

4  .   At  =  1.0  sec. 

In  performing  these  tests,  it  was  noted  that  the  use  of  the 
vehicle  length  speed  equation  to  establish  safe  following 
distances  may  not  be  realistic  in  the  case  of  trucks.   In 
fact,  if  this  program  is  actually  utilized  to  generate 
bridge  truck  loadings,  two  equations  should  be  established 
with  different  parameters,  i.e.,  one  for  trucks  and  one  for 
autos . 

The  results  of  the  test  case  for  SDFAC=5  compared  to  the 
control  case,  SDFAC=15,  clearly  indicated  measurable  effects 
occur  in  the  generated  truck  load  events  on  the  pseudo- 
bridge  when  this  parameter  was  varied.   Directly,  a  change 
from  15  to  5  for  a  speed  of  60  mph  changes  the  safe  following 
distance  from  117  feet  to  352  feet.    This  change 
implies  one  or  both  of  two  occurrances, 

1.  Lighter  vehicle  density  is  imposed,  and/or 

2.  Slower  event  dynamics. 

This  case,  because  of  its  tendency  to  widely  separate 
vehicles,  increased  the  number  of  load  events  generated 
to  132.   This  was  to  be  expected. 

Observations  of  the  effects  of  using  SDFAC=5  as  opposed 
to  SDFAC=15  include: 

1.   Truck  behavior  was  mostly  either  exactly  the  same 

or  relatively  minor  position  changes  occurred  at  the  same 

points  of  time. 
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2.  Auto  content  in  an  event  varied  considerably  from 
varying  in  number  and  position  to  being  exactly  the  same. 

3.  Relative  positional  relationships  were  not  pre- 
served. 

4.  Ordering  of  vehicles  was  not  preserved,  including 
occasional  shifts  in  truck  ordering. 

5.  An  initial  lag  of  one  second  occurred  for  this 
test  case,  but  was  not  consistently  maintained  due  to  the 
creation  of  a  new  event.  The  lag  was  restored  subsequent 
to  the  new  event,  i.e.,  over  the  control  case. 

6.  The  first  time  point  in  each  event  was  usually 
the  same  in  both  cases  in  terms  of  the  truck  present  and 
its  position,  i.e.,  considering  the  imposed  lag.   Total 
structure  of  the  event  initially  varied  because  of  auto 
behavior. 

7.  While  the  initial  time  point  might  provide 
duplication,  in  some  instances,  invariably  a  later  time 
point  in  the  event  would  evidence  variation  of  the  pseudo- 
bridge  load  in  terms  of  position  and/or  content. 

Decreasing  the  SDFAC  value  created  a  restructuring  of 
load  events  due  to  changing  the  safe  following  distance. 
This  may  impose  changes  in  traffic  dynamics  which  were  not 
straightforwardly  and  easily  predictable. 

The  results  of  the  test  case  for  SDFAC=10  compared  to  the 
control  case,  SDFAC=15,  were  very  similar  as  those  for 
SDFAC=5.   The  observed  effects  were  pronounced  and  were 
of  the  same  nature  as  the  previous  case.   This  value  of 
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SDFAC  establishes  a  safe  following  distance  at  60  mphf  for 
an  auto  of  234  feet  when  compared  to  the  control  case.  The  same 
number  of  generated  load  events  were  generated,  131.   However, 
the  time  structure  and  content  were  not  the  same. 

Compared  to  the  two  previous  cases,  the  bridge  load  data 
tended  to  be  more  like  the  control  case  than  the  SDFAC=5 
case.   More  events  were  identical  in  time  and  space  to 
the  control  case  than  the  previous  case.   Vehicles  tended 
to  be  more  frequently  in  the  same  position,  at  a  given  time, 
as  that  shown  in  the  control  case,  than  the  previous  case. 

The  observations  were  completely  consistent  with  increasing 
the  following  distance  between  vehicles. 

The  results  of  the  SDFAC=20  case,  compared  to  the  control 
case  also  clearly  evidencedvariations  in  load  event  number, 
vehicle  content  and  structure.   This  value  of  SDFAC  produced 
a  safe  following  distance  of  88  feet  at  60  mph  for  an  auto. 

The  sample  of  output  analyzed  indicated  the  following: 

1.  Repetition  of  exactly  the  same  load  event  as  the 
control  case  appeared   in    about  25%  of  the  events.   This 
case,  on  the  basis  of  repetition  of  exactly  the  same  event, 
was  less  similar  to  the  control  case  than  for  the  SDFAC=10 
case . 

2 .  Auto  variations  in  position  were  more  frequent 
than  truck  variations.   However,  truck  variations  were 
plentifully  present. 
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3.   Event  structure  usually  varied  to  the  inclusion 
of  an  auto  in  the  event  for  this  case  which  was  not  present 
in-  the  same  event  of  the  control  case. 

This  case  generated  133  events  as  compared  to  131  for  the 
control  case  and  was  contrary  to  the  expected  and  generally 
indicated  decrease  in  the  number  of  events  with  an  increase 
in  the  magnitude  of  SDFAC. 

The  test  results  for  SDFAC=25,  as  compared  to  the  control  case, 
and  the  SDFAC-2  0  case  were  evaluated  over  a  sample  of  the 
runs  generated.   Unlike  previous  cases,  it  was  felt  that  this 
case  should  be  compared  to  the  SDFAC=20  case  as  well  as  the 
control  case.   This  was  based  upon  the  fact  that  the  change 
in  the  safe  following  distance  was  very  small  as  compared 
to  the  earlier  cases,  e.g.,  see  Table  19. 

As  anticipated,  of  the  sample  analyzed  28.6%  of  the  time  points 
were  exactly  the  same  for  the  SDFAC=20  case  and  this  case. 
This  compares  to  11.9%  of  the  time  points  sampled  being 
exactly  the  same  as  the  control  case.   Further,  another 
28.6%  of  the  time  points  were  very  similar  between  this  case 
and  the  SDFAC=20  case.   The  same  comparison  between  the 
control  case  and  this  case  revealed  only  14.3%  of  the 
time  points  being  very  similar.   The  remainder  of  the  sample 
points  varied  mostly  due  to  changes  in  auto  position  in  the 
event  or  the  auto  content.   A  few  had  significant  truck 
position  variation  and  content. 

The  variations  noted  for  this  case  were  of  the  same  form 
observed  previously  for  the  other  cases  in  this  set.   The 
small  variation  in  the  safe  following  distance  between 
this  case  and  the  SDFAC=20  case,  e.g.,  18  feet  at  60  mph 
for  autos,  still  caused  significant  variations  in  the 
descriptions  generated  of  the  load  events. 
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As  the  following  distance  is  allowed  to  decrease  there  should 
be  a  tendency  to  have  more  vehicles  in  an  event.   Thus,  a 
reduction  in  the  total  number  of  events  should  be  anticipated. 
As  can  be  seen  in  Table  20.,  there  appears  to  be  an  overall 
weak  tendency  to  decrease  the  total  number  of  events  as 
SDFAC  increases.   This  is  due  to  shortening  of  the  safe 
following  distance  and  a  consequential  reduction  of  the 
distances  between  trucks.   This  then  has  the  tendency  to 
string  platoons  into  larger  platoons,  in  terms  of  the 
definition  utilized  in  the  program. 

The  simulation  of  vehicle  dynamics  by  the  BRIGLD1  simulator 
is  extremely  sensitive  to  Variations  in  the  value  of  SDFAC. 
Consequently,  load  events  are  similarly  sensitive  to  changes 
in  this  value.   At  equal  incremental  values  of  SDFAC  the 
change  in  the  safe  following  distance  increases  significantly 
as  the  value  of  SDFAC  decreases,  i.e.,  it  varies  as  the 
reciprocal  of  SDFAC.   Larger  values  of  SDFAC  will  cause 
proportionately  smaller  changes  in  the  safe  following  distance 

It  is  felt  Lhat  as  indicated  earlier,  a  better  representation 
of  safe  following  distance  is  needed  in  the  simulator, 
as  opposed  to  the  one  simple  equation  presently  used. 
The  choice  of  this  value  materially  affects  the  vehicle 
dynamics  and  acts  as  a  constraint  on  the  vehicle  statistics. 
The  mathematical  model  of  the  safe  following  distance  should 
not  only  consider  vehicle  differences,  but  also  traffic 
behavior  differences  as  a  function  of  a  given  highway. 

Extreme  care  must  be  used  in  the  variation  of  this  parameter. 
Studies  of  bridge  loading  effects  on  a  given  highway  should 
consistently  utilize  the  same  value,  and  it  should  relate 
to  the  traffic  behavior  of  that  highway. 
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TABLE  20  -  Platoon  Distribution  of  Bridge  Captured  Loads 


PLATOON  COUNT 


PI 
lC 

atoon 
Size 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Tot. 

sdf; 

5 

62 

28 

15 

10 

9 

4 

2 

0 

0 

1 

131 

10 

61 

29 

12 

10 

11 

4 

2 

0 

0 

1 

130 

15 

61 

23 

17 

13 

9 

5 

1 

0 

1 

0 

130 

20 

63 

24 

16 

14 

8 

5 

1 

0 

1 

0 

132 

25 

60 

25 

15 

13 

8 

3 

.de 

0 

1 

0 

129 

■"rcitcnoayc  xii^j 

5 

47. 

3 

21, 

4 

11. 

4 

7.6 

6.9 

3. 

0 

0 

0 

0 

0.8 

100 

10 

- 

46. 

9 

22. 

3 

9.: 

• 

7.7 

8.5 

3. 

1 

0 

0 

0 

0.8 

100 

15 

46. 

9 

17. 

7 

13. 

1 

10.0 

6.9 

3. 

8 

0.8 

0 

0. 

8 

0.0 

100 

20 

47. 

7 

18. 

2 

12. 

1 

10.6 

6.1 

3. 

8 

0.8 

0 

0. 

8 

0.0 

100 

25 

46. 

5 

19. 

4 

11. 

6 

10.1 

6.2 

3. 

1 

2.3 

0 

0. 

8 

0.0 

100 
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DEVELOPMENT  OF  THE  DYNAMIC  STRESS  PROGRAM 

This  section  of  the  report  presents  the  basis  of  the  dynamic 
structural  analysis  program  which  was  developed  to  utilize 
the  generated  load  data  from  the  BRIGLD1  Simulator  and  to, 
in  turn,  predict  the  stress  caused  by  the  synthetic  loads 
in  a  given  bridge  design.  The  ultimate  purpose  of  the  joint 
use  of  these  two  programs,  as  indicated  earlier,  is  to 
provide  a  basis  for  the  prediction  of  long  term  fatigue. 

The  generated  load  data  are  assumed  to  form  a  force  vector 
which  varies  with  time,[F(t)],  and  is  the  driving  function  in 
the  nonhomogeneous  differential  equation. 

[F(t)]  =  [K]  [<S]  +  [mh  • 

where  [K]  is  the  unit  stiffness  matrix  for  an  "effective" 
beam, 

[6]  is  the  displacement  vector, 

m  is  the  element  mass  which  is  accelerating 

and  6  is  the  element  acceleration. 

As  is  implied  in  the  foregoing,  the  bridge  is  subdivided  into 
a  rectangular  grid  system,  elements.  The  longitudinal  center 
line  of  each  rectangular  grid  is  the  center  line  of  a  main 
girder.  The  longitudinal  width  of  each  grid  is  nominally  10 
feet.  The  time  dependent  load,  [F(t)J,  is  applied  to  the  geo- 
metric center  of  each  grid  element. 

Two  forms  of  basic  bridge  structures  were  assumed: 

1.  A  simple  span  bridge,  and 

2.  A  continuous  span  bridge. 
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The  simple  span  bridge  is  evaluated  on  the  basis  of  the 
deflections  of  the  effective  longitudinal  beams  when 
subjected  to  bending.   The  continuous  span  bridge  is  identi- 
cally calculated,  except  that  each  supported  section  of  the 
span  is  assumed  to  be  a  free  body.   This  imposes  the  need  to 
determine  the  discontinuity  moments  at  each  support  position 
and  reintroducing  these  discontinuity  moments  in  the  simple 
span  stress  calculations.   This  re-establishes  necessary 
continuity. 

The  solution  to  the  discontinuity  moments  and  shears  at 
the  intermediate  supports  are  calculated  via  the  "3  moments 
equations"  for  the  bending  fiber  stresses.   The  strength 
of  material  approach  utilized  herein  for  the  beam  element 
stiffness  matrix  provides  the  user/engineer  a  means  to  achieve 
better  understanding  of  the  stress  fatigue  problem.   The  direct 
stiffness  method  developed  for  the  analysis  of  assembling  the 
plate  and  rib  stiffness  is  certainly  comparable  to  the  beam 
element  stiffness  process  in  which  the  anisotropic  composition 
of  deck  plate  and  beam  rib  is  taken  care  of  as  the  effective 
modulus  of  the  bending  rigidity.   The  correct  stiffness  for 
the  representation  of  the  bridge  is  only  relatively  important 
so  long  as  the  load  distribution  causes  the  proper  results 
in  deflection  values. 

Structural  Analysis 

The  study  of  truck  loadings  on  a  bridge  deck  has  been 
investigated  by  many  authors,  besides  the  general  guideline 
presented  by  AASHO  methods.   The  objective  of  this  approach 
is  to  provide  a  basic  flexural  beam-element  analysis  method, 
suitable  for  a  computerized  analysis  of  a  bridge  design, 
undergoing  dynamic  loading.   The  complete  treatment  of  the 
finite  element  method  of  an  orthotropic  slab-beam  in  the 
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principal  directions  was  not  the  intent  of  this  project.   It 
was  more  convenient  to  formulate  an  orthogonal  beam  element 
(grid  network)  on  the  bridge  deck,  so  as  to  distribute  the 
loading  in  a  reasonable  fashion.   It  was  pointed  out  by 
Scordelis,  et  al7  that  the  empirical  formula  Nr7T  =  §-  wheel 
load  for  interior  girders  may  not  be  accurate  for  design 
use.   Thus  some  other  assumptions  were  made  and  are  described 
herein,  for  the  analysis  of  loadings  on  the  girders.   The 
exterior  girders  are  designated  to  carry  a  fraction  of  wheel 
loads  N    =   i/7,   (S.  is  equal  to  the  half-flange  width  plus 
the  overhang  flange  in  feet) . 

A  complete  treatment  to  the  force  displacement  via  the 

finite  element  method  was  much  too  cumbersome  to  achieve  the 

basic  objectives  within  project  constraints.   A  more  rational 

approach  was  selected  for  the  computing  of  the  deck  stiffness 

in  the  transverse  direction  by  considering  the  beam-plate  in 

the  transverse  direction  to  be  simple  supported  at  both  ends 

which  is  generally  considered  to  be  rigid  with  the  girder 

and  overhang  flange.   In  accordance  with  the  Timoshenko  and 

Wainowsky-Rrieger ,  the  plate  is  assumed  cross  stiffened  by 

two  sets  of  equidistant  stiffeners.   The  bridge  can  then  be 

assessed  for  bending  behavior  by  calculating  the  longitudinal 

and  transverse  plate  stiffnesses  separately  for  the  effective 

modulus.   The  determination  of  the  effective  modulus  for  the 

bending  stress  calculation  is  shown  later  in  the  report.   If 

the  computerized  program  is  to  be  developed  for  the  solution 

of  the  detailed  stress  distribution  in  the  vicinity  of  a 

wheel  load,  then,  a  finer  grid  will  be  needed  and  it  will 

be  necessary  to  utilize  the  anisotropic  slab  stiffness  of 

D  ,  D   and  D    to  account  for  applied  loading  distribution, 
x    y       xy  rir  ^ 

However,  a  modular  anisotropic  routine  can  be  added  for  this 
particular  interpretation  for  the  determination  of  the  stress 
distribution,  in  the  composite  section  of  concrete  slab  and 
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girder,  as  well  as  the  bents  or  diaphragm.   This  portion  of 

the  report  will  be  restricted  to  considering  the  main  deck  stress 


The  variation  of  the  EI  along  the  beam  must  be  taken  into  ac- 
count.  The  local  average  of  the  beam-slab  combination  is  used 
in  the  calculations  as  will  be  discussed  later  in  this  report, 
for  example, 

If      X.  =  X  .  ,  then 

11 

It  should  be  noted  that  this  forms  the  principal  diagonal  in 
the  assembled  stiffness  matrix  for  the  bridge  deck,  which  con- 
sists of  multiple  beam  elements  on  a  bridge  deck.  The  lateral 
beams  on  the  bridge  are  defined  as  the  deck  plate  with  inter- 
mittent synthetic  beams,  diaphragms,    connected  between 
the  main   beams.  These  synthetic  beams  or  stiffners  account 
for  the  beam  stiffness  in  the  lateral  direction,  much  the  same 
way  as  the  longitudinal  effective  beams.  These  synthetic  stiff- 
ners are  assumed  supported  at  the  edges  of  a  deck.   The  numeri- 
cal stiffness  is  computed  by  using  the  same  equation  as  above, 
when  the  lateral  plate  and  the  intermittent  rib  stiffnesses  are 
considered  as  the  effective  cross-section  of  the  beam. 

The  continuous  span  was  considered  in  different  manner  for 
the  stiffness  matrix  formulation.   The  influence  coefficients 
can  be  evaluated  by  considering  the  bending  moment  area  for 
the  main  girder-deck  system.   The  discontinuity  moment  at 
the  support  juncture  is  considered  to  be  contributory  to  the 
deflection  of  the  beams  and  can  be  determined  by  the  integral 
of  the  area  under  moment  diagram  M/EI.   Then,  the  deflection 
between  grid  points  can  be  represented  by  the  bending  moment 
area  integral  with  respect  to  any  grid  point  under  unit  loading 
at  the  beam  on  the  grid  point. 
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JAB 


D   /Mxdx\     B   /Mxdx\ 
A      \    El    J  D   \  EI  / 

Moment  Area 
EI 


(ih 


Moment  Area 
EI 


^   x. 


Simple  Support  Span  and  Continuous  Span 

Considering  a  linear  system  of  beam  rib  elements  along  the 
longitudinal  direction  simply  supported  at  ends,  the  influ- 
ence coefficients  of  a  simply  supported  beam  will  be  defined 
as  6. .,  i.e.,  the  deflection  of  the  beam  at  node  i  due  to  a 

unit  force  at  node  i .   6  .  .  will  be  the  deflection  at  i  due  to 

J     11 

a  unit  load  at  i.   The  simple  supported  beam  can  then  be 
subdivided  into  the  various  grid  points,  both  in  the  longitudinal 
and  transverse  directions.   For  the  longitudinal  beam,  6 . . =  =- 


ii 


at  node  x.  is  written  as: 

l 


x     1     x.  (£  -  x. ) 

0  .  .  =  rr    =   1   1 

ix     ii    61   EI    ' 


{  £2  -  x2  -  U  -xi)2} 


For  the  deflection  of  any  node  i,  due  to  an  axle  load  located 

at  node  j  where  j  is  situated  on  the  right  side  of  the  node 

1# 


can  be  written  as 

«•■  -i 

3-D    kij 


x.   U  -  X  .  ) 


This  represents  the  deflection  of  any  node  i  to  the  left  of  the 
axle  load  j  on  the  beam.   In  evaluating  this  deflection 


119 


for  any  node  along  the  beam,  it  is  important  to  note  that 
the  EI  are  the  composite  beam  stiffness  (plate  and  rib) , 
in  the  longitudinal  direction.   <$  /   is  the  deflection  due  to 
a  unit  load  at  point  A  with  respect  to  point  B  at  any  grid 
point  between  the  beam,  say  6. .  .   Then,  the  inverse  of  the 
influence  coefficient  matrix  will  be  the  stiffness  matrix. 
6 . ,  is  the  influence  coefficient  of  the  grid  system  and 

[K]  =  [6]"1 

Once  the  moment  diagram  on  each  span  of  the  beam  is  calculated, 
then,  the  reaction  force  may  be  produced  in  a  similar  way 
using  the  discontinuity  bending  moment  method;  i.e., 

R   =  R'   +   R"   +   n-1     n      n     n+1 

n    n      n -  — — 

*n  *n+l 

R1  =  Reaction  force  due  to  load  on  span  n 

n  r 

R"  =  Reaction  force  due  to  load  on  span  n  4-  1 

'I       1    R" 


777/7  7777 


Rn  ■        ■    -n 


Where  R'    = 

n 


R"    = 


P   b 


n 


Pn+1  An+1 


B  Pnbn      +      Pn+l^n+l   +  Mn-r  Mn      .    -Mn   +  Mn+1 

n  '     ~1~  i    „  C    "  I 


n 


n+1  n  n+1 
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P  b     M   .  -  M 
n  n  +  n-1    n 


n+1  n+1   n    n+1 


n 


n+1 


For  a  continuous  span,  the  loading  on  bridge  deck  can  be  utilized 
to  calculate  the  moment  on  the  main  girder,  or  beam.   The  total 
summarized  moments  can  be  computed  with  the  procedure  outlined 
below. 

T> 


M 


\/ 


"5F 


^k: 


7T 


M. 


M 

3    U3 


7T 


M4 


M 
6   5 


1. 


Calculate,  M,,  M- ,  M  based  on  the  discontinuity 


moment  matrix  from  loads  P,  ,  P?/  and  P_  independently  on  each 

span.   Then,  sum  up  the  moments  M   ,  M   ,  an£  m   from  the 

Pi   P2       P3 
loads  P1/p2'  and  p3  on  the  various  spans.   The  total  discontinuity 

moment  at  any  given  time  t. ,  on  support  1   will  be  equal  to 


Dl- 

(M                  + 
Pi 

M 
P2 

+   M 
P3 

+   M       ) 
pn 

= 

P=P 

z    n 

P=l 

M 
Pn 

(t,) 

P2           1 

E  (M      + 
Pi 

M 

P2 

+   M 

P3 

+   M      ) 
P4 

If  more  than  one  load  exists  on  any  span,  the. moment  calculation 
remains  the  same  as  above,  and  is  essentially  additive. 


2.  The  moment  due  to  any  load  P.  distribution  on  any 

span  I     will  be  computed  as  if  the  span  is  simply  supported. 


M 


,  =  Moment  due  to  load  P .  =   i  1   1  3- 
0-1  l    5 

in  span  0  -  i 


X 
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^ 


X- 


p. 

1 


-* 


X. 


■> 


-&- 

% 

M 


ik. 


~~ 


M 


M 


> 

X.  I 

l         ■ 


"2^ 


M 


M 


M 


D. 


M. 
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m     =   p.x.  a  -x., 

n      J-  J-  n ij_  ,  and 


£ 


n 


,th 


is  the  moment  due  to  load  P.  of  the  n   span. 

i 

load  exists  on  span  n,  then 

n 

i  =  1 


If  more  than  one 


M  = 


Vi^n-Xi) 

*n 


3. 


The  total  moment  distribution,  on  span  in, 


including  the  discontinuity  moment,  is  determined  from 


M 


n 

=  Z 
i  =  1 


p..x..(*  -  x.>   iyL   (£n-x. ) 

l  i   n     i  '■      +      D    n   l 


n 


n 


M      (&    -  X    ) 
+  D   •  A  n-i   Ai~l' 


n"1 


n  -1 


For  example: 


M .  due  to  P . 
l         l 


i+1 


In  general,  the  discontinuity  moment  MD   is  negative  in  value/ 
except  the  complicated  case  when  there   is  more  than  one  load 
acting,  for  example /  loads  acting  on  an  adjacent  span  which  may 
cause  an  inflection  in  the  curvature  of  span  i. 
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4.  Bending  stress 


o      (stress)  =  Mi   =  Mi 

a 


Z.      EI  . 
on  steel  girder  1 


E  C0 
s  2 


aB  (concrete  girder)  =     M-jy  EIci 


EcC2 


Each  end  on  the  continuous  span  is  assumed  to  be  simply 

supported,  if  it  is  a  built-in  end,  one  more  equation 

expressing  the  condition  that  no  end  rotation  will  take  place 

at  the  support.   The  angular  rotation  vanishes  at  the  built-in 

end.  n  =  0 

6=0=  angular  slope  at  built-in  end 

Left  Span,      e'  =  Vn  +  ^-l31  n  +  AnCn 

n    3EI     6EI  ET~ 

n      n       n  o 

p.  ..  c         „■■_  Vn+1    Mn+l*n+l    An+ldn+l 
Right  Span,      0n  -  ^r~   +  + 

n+1        n+1      n+1   n+1 


To  compute  the  total  bending  moment,  the  moments  due  to  each 
load  must  first  be  computed  for  the  span,  that  is, 

M^ (x)  =  moment  due  to  axle  load  1, 


M2 (x)  =  moment  due  to  axle  load  2, 


M^ (x)  =  moment  due  to  axle  load  3, 


v  m    =   total  bending  moment 
B 


=   Mx(x)+M2(x)+M3(x) 
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Where , 


P. a,b  x 

M1(X)  =  ^~ 


Plblx 


,  for  x,  <a1 


or 


PlalblX 


for  a,  <  x  <  (a,+b,)  «  PlalX 


ll  "1 


M2(x)  =   P2a2b2X 


ha2 


or   P2a2X 


2  2     for  x  <  a. 


for  a2<x<l1 


M  (x  )  +  M2(x2) 


plbi 

P2b2X 

1 

x    i        c 
*1 

P,a,  „  ,  P~b~x 
or   1  1  x  +   2  2 


I 


U 


(valid  up  to  axle  load  P  ) 


(valid  for  a,<x  <a_ 


Plal     P2a2 

— x  +       x  (valid  for  a0<  x    ) 


Ul         ili  /        £i 

fPlal+P2a2  \ 
V   %      /x, 


+P2b2 


x   ,  for  x<a. 


for  a, <x<a_ 


for  a~<x 


To  determine  the  centroid,  C  ,  of  the  bending  moment 

n 

diagram, 


E  Area  of  Moment  diagram  x  centroid  distance  to 


C   =   point  "0 

n     c 


Total  area  of  moment  diagram 


To  compute  A  =  Area  of  moment  diagram,  due  to  loads  on 
span  n, 
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A„  =  IE  v 

5.  To  compute  the  discontinuity  moments  M  ,  M  ,  , 

on  span  n. 

at  n  =  -1,  M  =0  for  the  simple  supported  case, 

and  for  the  built-in  end  at  support  "0",  an  additional  equation 
is  required: 

e0  =  o, 

.  0  _  e  _  Vi  +  Vi  .  Aibi 

•#U  U0         3EI-L  6EI1        j5~ElT 

-M1         3A1b, 

Calculate:   Mn  =  — ??—  -  —* 

0     2    2*1 


where  A.  ,  b..  ,  &  1  ,  are  known  from  previous  calcula- 


tions . 


The  end  force  or  reaction  is,  then 

m   -M    -M  +M  ... 
r  =  R   +  p    +   n-1   n  +   n  n+1 
n    n   ~n     T^     3T^~ 

For  a  simple  supported  continuous  span,  the  discontinuity  mo- 
ments at  the  supports  can  be  represented  as 


M  H  M   ,A     A  a 

1     n  n    n-1  n    n  n 


Gn   TEI~    '      6EI    '  o„EI 
n       n    ^n  n 


••   Vn+1   «n+l*n+l   An+lbn+l 


+ 


+  I-    ..EI 


or 


n         3EIn+1  6EIn+1         -»+1--n+1 


t  I    I 

e    =  -  e 

n  n 


,      1    fyLS   +2  M     fe  +  >™±) 
n_1   Vn-l)  n  Vn        Vl/ 


+  M 


n+1 


h+1   I 


n+1 


6A  A  6A  x,b    ,- 

n  n  n+1  n+1 

n  n  n+1  n+1 


k    2(A 


M 


h-1. 


n+l'     '        n+1 


}  k 


M 


n+1 


"6Anan 
*nXn 

6An+lbn+l 


U 


n+1  n+1  J 


For   the  n  =   4   span  case, 


V  2  /*  2+"  3\ 

n=3,   _|_3  +  2M3y_ij 


m2    2 


v3 


4    4 


6Alal         6A2b2 


*1X1 


6A2a2 

y2 


6A3a3 
£3I3 


Aol 


2"2 


6A3b3 
I3I3 


6A4b4 


M3% 
n  =    4,    -|-i  +   2M4 


MA  .    +    2MC 
4    4  5 


felt) 


Vs 


+  Ve 


6A4a4 

K..I. 
4    4 


6A5b5 
I?7~ 


=    0 


Case    I.    Ends  w/variable   fixity 


6=0 
o 


If 


Vi 


Ml*l 

21. 


3Albl 
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2    /A    +* 


.  A 


2 

£     .A 
_2+_3\    _ 

,I2    I3/    I3 

H  Vl3    I4/    ^ 

A  2C*4> 


0 
0 


0  *3 


0 


0  rA 


M 


M. 


M, 


M. 


M 


=    -6 


/Ai*l  A2b2\ 


U2 

ft 


A3b3 


) 


3X3 


A,a  „ 
4    4 


4*4 


J£±\ 


V 


-M; 


Since  M.    =   0      if   simple   supported  end 

M       =    0 

o 

Solve  M    ,    Mx,    M2,    M3,    and  M4 . 


Case  II.  Simply  supported  ends 


If  simply  supported  at  ends  "0"  and  "4 


M  =  M„  =  0 
o     4 


3A3k3 
£3X4 


J 


A 

£1 

0 

l2 

0 

3 

0 

4 

A 

000 

///// 

///// 

///// 

///// 

000 

///// 

Then, 


x 


0       £_2      2/_?3+-^.\ 
X2       VJ3  V 


~M1~ 

M2 

M3 

„     _ 

=  -6E 


x 


Alal 

Z — T~ 

11 

A2a2 
^2*2 

A3a3 
^3*3 


A2b2 
2X2 

A3b3 


£3I3 

A4b4 
£4X4 


128 


Solve 


M. 


M. 


LM3 


=    -6 


7VV\  li        o 

h    \Z2    h)  h 

~2  V3    V 


~      1 


Aiai 

11 
A2a2 

rf 

2X2 

A3a3 

3X3 


+ 


A2b2 

12 

A3b3 

*~^ 

3  3 

4  4 
X4    4 


General   n   spans    continuous   Beam  Bridge  with   end 
constraints . 


IK 


h 
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9   A+^ 
Xl    J2 
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-^=-  0 


2    f2+il 
X2    *3 


I 


V*4 

2    J3    Z4 


.    0 


7    *4   +?J 
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n-1 


n-1 


n-1      n        n   Q 
I        I 
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n        n 
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M. 
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A    b       &     t 
W    ^1 
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=  -6 


r 


Alal  +  A2b2 


I    T    A 
11 


2I2 


A2a2  +  A3b3 


*2*2 


3*3 


A3a3    A4b4 


3  3 


4  4 


' 


A   -.a*      A  b 
n-1  n-1  ,   n  n 

I      71  T  +  I    1 
n-1  n-1    n  n 

A  a 
n  n 

n  n      _J 


Anisotropic  Composite  Deck 


i  ■  ' 

a  =  E  e  +E  e 
x    x   x   1  1  y 

a  =  E  '  e  +E  e 

y   y   y     x 


r  = 

xy 


d  = 

x 


Gy 


xy 


E   .3 
x  h 

12 


«   =   E   .  3 
Dy      Y   h 
y      12 

1  '  3 
D,  =   E   h 
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3 
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Neutral  Axis 


D 


ax' 


+   2(D.,+2D   ) 


3^00 


3x23y2     °Y 


By* 


=  q 


H  =  D1+2Dxy 


For  isotropy/ 


Ex  "  E"  = 


y     i-v 

Eh3 


.   E 


VE 


.  G  = 


1-v 


*   D   =  D   = 

"    x      y     i2d-v2) 


2(l+v) 


H  =  D,+2D    =  -^ 

1    ^     12  (1-v2) 


For  reinforced  concrete  beam, 


n  =  Es/Ec    , 


V    ~ 

c  ~ 


E 


E  'E  ' 
x  y 
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For  two-way  reinforced  concrete, 
E_  E 


E        r                    E      r       E 
Dv  =  -£ i  +(n-l)l    =  2 U  +  (_s  -i)i  I 

X     l-v  2    L  CX  SXJ       x_v  2    L  C*     Ec         S*J 


E       r 
D  =  _£ i   +(n-l)I 

Y   i-v  2   [  cy  v     syj 


1     c    x  y 


D    =  1-v   /n  \/D~D 
xy      c  '  2     v  x  y 


where   I   =  moment  of  inertia  of  slab  material 
ex 

I_„  =  moment  of  inertia  of  steel  reinforced  rods 
sx 

y 


/////////// 


D   = 
x 


D   = 

y 


E  a,h    1 

c  1 


/ 


12 (a  -t+a3t) 

V 


a  =  h/H 


Dl  " 


I  =  moment  of  inertia  of  a  Tee  section  of  width  a. 
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D   =  Bl 

X     r— 


H  =v6~D 
x  y 


B. 


D   = 


y   a. 


B..  ,  B-    are  flexural  rigidities 


a 

a-M 

"■ 

■ 

kJ 

1 

0 

1 

1 

— 



# — 
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Calculate  the  beam's  sectional  modulus 
Steel  Girder 

Calculate  the  geometric  centroid,  i.e.,  neutral  axis. 

Given, 

b  =  effective  width 

t,=  concrete  deck  thickness 

I  =  moment  of  inertia  for  WF 

s         .  i+ 

beam  — 

A  =  sum  of  the  cross  sectional  area 

s 

of  reinforced  steel  bars 
Bars 
A  =  I  Beam's  Area 

EA  C 
C,  =   i  i  (Sum  of  each  cross  section  Area  centroid  distance) 

Total  Cross  Sectional  Area 


A.  =  Area  of  Material i 


C.  =  Centroid  Distance  with  Reference  to  Base  Line  of 

i 

Material  i 

b,tn  .h,,   „  ,«  „    ,.     ^  ,t2 


;.  c    =    1  1,r^As-cs+Arci+b2t2(i-+cP) 


Area  of  Concrete  (b,t,  ) +b~t2  (coverplate+A,.  (WF=Beam)  ) 

st . 


Combined  moment  of  inertia, 

Jc  "  15-  blfcl3  +  ^'r'vi'2*   I.t.bar.+Vcrc.»2   + 

.     h  +  Wci)2  +I5-  b2fc23  +  Va-  (cc-ci)2 
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Z    '    =  —       Sectional  modulus   for  beam's   upper  element 
cl        cl 


Z_     =  : —       Lower  sectional  modulus   of  beam 
C2         c2 


3     5 


*    6   *     &     ♦     £>      u,    <  >    a       0  •    ^      »  •  Q    &  &#  <j 


i   vzznmnzzj 


Neutral  Axis 


\„ 


-A 


/ 


/ 


C>      &         O         V> 


;i 


V 


fejj 


-<> th- 


where  CT  =  distance  of  WF  beam  centroid  to  the  base 
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Concrete  Girder 


For  a  concrete  girder; 

b 


°c  °s/n 

led-    aPEd" 


n 


kd 


<lik>/k  ■    SV 


n 


7ec  =  7 


kd 


bkd  2-  =  nAs  (d"kd) 


P  - 


'bd  Area  ratio  of  steel  to  that  for  concrete 
above  the  steel  reinforcing  bars 


k  =   V2pn+(pn)    -  pn 

C  =  ^compressive  force 
bkd 


-w 


a   a 
T  =   s  s  =  tensile  force 


M      = 


n     = 


o   A   jd,  M      =   i>    [a   f   d(l-0.6p   ^Xr)j 

s    sJ    '  uYLsy  c     -* 


_f_     =      30x.CT 
Ec  4.3X106 


=   7 


0.9  mod.    factor 
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Mc  =  f-   b(kd)  Cjd)  -  |  ac  kjbdZ 


M  &  M  are  then  the  bending  moments  of  a  concrete 
s    c 


beam 

To  determine  n  =  — — 


a      the  compressive  concrete  stress  is  evaluated,  or 


is  assigned,  and  is  defined  as 

a. 


where 


a  =  J5. 
c    n 


a   &  n  are  known,  and 
c  ' 


a   is  computed. 


Find  the  neutral  axis  of  the  reinforced  concrete  beam,  i.e., 

k  =  [V2pn+(pn)2  -  pnj 

A 
where  P  =  gf 

calculate  kd,  i.e.,  neutral  axis  from  the  base  of 
composite  beam 

M  is  bending  moment  in  the  beam  at  any  point  along 

the  beam,  and 

M 
B 
Bending  stress  =  = — ,  where 

Zec 
Z   =  equivalent  bending  sectional  modulus  for  concrete 

Now         Z   =  equivalent  sectional  modulus  for  steel  section 
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■[' 


EI  =  E  I  ,   +  EI 
c  g/s    s  s 


I   =  moment  of  inertia  of  gross  concrete  section 
g 

about  the  controid  axis,  neglecting  the  reinforcement 
steel,  and 
R  =  ratio  of  maximum  design  dead  load  moment  to  maximum 
design  total  load  moment,  always  positive. 
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DESCRIPTION  OF  ANALYTIC  METHODOLOGY  DEVELOPED 

The  objectives  of  the  investigation  of  analytic  methods  was 
to  attempt  to  develop  and  provide  a  practical  means  for  pre- 
dicting the  long  term  effects  due  to  heavy  truck  loadings.   The 
methodology  which  is  based  upon  the  use  of  a  traffic  simulator 
is  inadequate  for  long  term  applications,  for  example,  50-year 
life  spans.   In  order  to  utilize  such  a  simulator  small  time 
samples,  for  example,  a  continuous  two-week  sample,  must  be 
taken  for  traffic  distributions  anticipated  for  a  given  period 
of  time.   It  is  then  necessary  to  run  several  such  samples, 
varying  the  time  dependent  parametric  traffic  data.   The  present 
version  of  BRIGLD1  can  be  satisfactorily  operated  at  up  to  a 
720  to  1  real-time  to  computer  use  time  compression  ratio,  as 
compared  to  the  2  to  1  of  the  original  program.   However,  even  at 
this  increased  rate  a  total  of  5  hours  of  CP  time  is  required 
merely  to  generate  the  loads,  without  any  consideration  of  the 
stress  calculations.   The  dynamic  stress  calculations  are  far 
more  complex  and  require  a  very  small  time  interval.   Hence, 
the  CP  time  required  to  calculate  the  response  of  a  bridge  to 
the  prescribed  loads  will  be  far  greater  than  that  necessary 
to  generate,  the  loads.   Present  range  from  6  to  1,  to  30  to  1 
for  real  time  compression  to  perform  the  stress  calculations. 
This  implies  the  need  for  an  additional  112  to  560  hours  of  CP 
time.   The  main  cause  is  the  integration  interval  size.   The 
load  generator  is  essentially  a  pseudo-truck  load  event  simu- 
lator and  will  integrate  motion  at  a  maximum  interval.   This 
interval  is  far  too  large  to  be  satisfactory  for  the  calculation 
of  a  bridge's  structural  response  to  a  load.   Consequently,  the 
structural  program  moves  the  axle  defined  loads  over  the  deck 
at  a  necessarily  small  interval,  independent  of  the  simulator's 
incremental  interval.   The  stresses  are  calculated  at  each  of 
these  time  points. 

The  two  programs  are  independent  of  each  other,  that  is, 
stand  alone.   They  can  be  linked  together  or  run  totally 
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independently.   The  output  of  the  load  generator  may  be  used' for 
other  purposes  and  the  input  to  the  structural  program  may  be 
derived  from  other  sources.   Additionally,  if  a  pseudo  bridge 
length  of  sufficient  magnitude,  which  would  include  all  bridge 
lengths  to  be  analyzed  structurally,  is  used  by  the  simulator, 
there  is  no  need  to  re-run  the  load  generator  for  shorter 
bridge  lengths.   The  structural  program  will  properly  use  the 
one  set  of  data,  recorded  on  magnetic  tape,  for  the  analysis 
of  all  bridge  lengths  equal  to  or  less  than  the  original 
pseudo  bridge  length  used  by  the  simulator,  for  the  highway 
defined  by  the  simulator's  data  base. 

While  the  above  defined  approach  has  increased  the  useability 
of  the  BRIGLDl  simulator  and  the  orginal  concept,  which 
caused  its  creation,  this  approach  is  not  economically  a 
feasible  approach  to  provide  insight  into  all  possible  forms 
of  significant  truck  platoon  events  or  into  long  term 
accumulated  effects. 

A  more  practical  approach  is  to  ignore  traffic  simulation, 

merely  generate  each  unique  platoon  event  and  evaluate  its 

effect  upon  a  bridge.   Subsequently,  generate  long  term 

effects  by  operating  on  each  single  event  effect  with  the  incidence 

of  the  event  over  the  period  of  interest.   However,  even  this 

approach  cannot  be  practically  applied  against  the  number 

of  potential  platoon  events,  as  is  discussed  below. 

Presently,  the  simulator  has  a  total  of  10  truck  types 
defined  in  its  present  highway  traffic  data  base.   This  is 
less  than  the  actual  total  number  of  types  in  existence, 
for  example,  the  set  of  truck  types  utilized  in  the  far  western 
states.   However,  using  the  10  types,  as  defined  below: 
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1. 

Type 

2D 

2: 

Type 

3 

3. 

Type 

2S1 

4. 

Type 

2S2 

5. 

Type 

3S2 

with  one  length.  =1 

with  one  length  =  1 

with  three  lengths  =  3 

with  three  lengths  =  3 

with  two  lengths  =  2 
We  have  a  total  of  10   elements 

Assuming,  that  10,000  lb  load  increments  will  normally 
provide  a  significant  variation  in  the  response  of  a  bridge, 
and  that  no  significant  response  can  be  expected  below  a 
gross  load  of  10,000  lb,  the  above  number  of  elements  increases 
to  the  following : 


1. 

Type  2D,  1  x  2  wt 

classes 

=  2 

2. 

Type  3,   1x3   " 

ii 

=  3 

3. 

Type  2S1,3  x  3   " 

H 

=  9 

4. 

Type  2S2,3  x  15  " 

it 

=15 

5. 

Type  3S2,2  x  8   " 

ii 

=16 
4T 

elements 

If  the  effects  due  to  speed  are  considered  in  an  essentially 
minimal  manner,  that  is,  that  20  mph  speed  increments  will 
provide  a  significance  threshold  in  the  dynamic  behavior 
of  a  bridge  and  that  a  lower  speed  limit  of  20  mph  is  used, 
the  above  number  of  elements  increase  to  the  following: 

1.  Type  2D,   2x3  speed  classes  =  6 

2.  Type  3,    3  x  3      "       "  =9 

3.  Type  2S1,  9x3      "       "  =27 

4.  Type  2S2,15  x  3      "       "  =45 

5.  Type  3S2,16  x  3      "       "  =48 

135   elements 

Simplistic  approximation  of  the  number  of  unique  single  truck 
platoon  load  events  may  be  determined  by  taking  all  combinations 
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of  the  indicated  elements  for  platoon  sizes  of  interest. 
T6  provide  a  reasonable  basis  for  platoon  definition,  it 
will  be  assumed  that  a  platoon  continues  to  exist  as  long 
as  the  bridge  is  excited  by  a  truck  load.   Hence,  a  truck 
platoon  size  of  six  trucks  appears  to  be  a  maximum 
possibility.   If  all  unique  platoon  events  were  to  be 
determined,  then,  it  would  be  necessary  to  consider  spatial 
configurations  of  the  component  trucks,  that  is,  spatial 
or  positional  permutations.   This  leads  to  an  immense  number. 
Merely  considering  combinations  indicates  the  following  number 
of  unique  platoons  (not  platoon  events) : 

1.  Single  Truck  Platoons:  135 

2.  Double  Truck  Platoons:  9,045 

3.  Triple  Truck  Platoons:  400,995 

4.  Quadruple  Truck  Platoons:        13,232,835 

5.  Quintuple  Truck  Platoons:        346,700,277 

6 .  Sextuple  Truck  Platoons :       6 ,356, 839,335 

6,717,182,522 
6,717,182,522  unique  truck  platoons  is  far  too  many  to  allow 
practical  evaluation  of  the  response  of  a  bridge  to  each 
platoon,  ignoring  positional  permutations,  especially  via  a 
traffic  simulator. 

Hence,  a  critical  analysis  of  the  definition  of  significant 
events,  that  is,  significant  truck  loading  effects,  must  be 
made  and  an  attempt  made  to  reduce  the  total  number  of  unique 
truck  platoons  load  events  to  a  manageable  size. 

A  more  realistic  approach  to  providing  a  long  term  cumulative 
damage  estimate  to  bridges,  is  to  start  with  a  definition 
of  significant  stress  class  intervals  and  dynamic  stress 
range  class  intervals  to  provide  a  basic  definition  of 
significant  events  to  a  bridge.   A  reasonable  basis  would 
be  to  use  10  percent  or  10  intervals  of  maximum  stress.   This 
leads  to  a  total  10  types  of  maximum  stress  events. 
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Having  defined  the  significant  event  it  is  then  necessary 
to  define  the  conditions  which  cause  the  event,  that  is , 
in  this  case  the  set  of  trucks,  in  terms  of 

1)  type  (including  configuration  and  length) 

2)  gross  weight 

3)  speed 

4)  platoon  configuration 

The  approach  to  this  is  to  predetermine  which  combinations 
are  equivalent  under  constant  nominal  structural  conditions, 
and  select  a  representative  truck  platoon  load  case  for  each  10 
percent  stress  range  interval.   Each  of  these  10  cases  can 
then  be  run  on  the  structural  analysis  program  to  generate 
the  specific  response  for  each  bridge  of  interest. 

The  total  incidence,  for  any  given  time  period,  would  be 
made  up  by  the  incidence  of  each  truck  platoon  element  in 
each  10  percent  equivalence  set  over  the  time  period  of 
interest.   This  reduces  the  structural  calculations  and  the 
load  generation.   It  imposes  the  need  to 

1.  Segregate  the  different  truck  platoon  events 
into  the  10  percent  stress  interval  equivalence  sets.   This, 
in  itself,  is  a  very  large  job,  although  of  a  non-recurring 
nature . 

2.  Calculate  the  incidence  of  each  separate  truck 
platoon  event  in  each  equivalence  set  and  sum  these  events 
over  each  equivalence  set  to  determine  the  total  incidence 
in  each  10  percent  stress  interval. 

While  the  foregoing  provides  a  better  basis  for  generating 
long  term  histograms,  of  stress  classes,  it  is  still  almost 
an  impossible  task  to  form  the  necessary  equivalence  classes 
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of  truck  platoon  load  events.   Hence,  another  approach  was 
conceived  that  appeared  to  have  significant  merit  and  was  of 
a  general  nature,  which  would  allow  broad  utilization. 

This  approach  is  predicated  upon  the  concept  that  each 
axle  has  a  stress  signature  and  that  a  total  truck  signature, 
or  platoon  signature,  is  actually  a  composite  of  the  individual 
axle  stress  signatures.   Further,  it  is  predicated  that  the 
amplitude  is  a  first  order  effect  of  the  axle  weight  and  the 
duration  is  a  first  order  effect  of  the  axle  speed. 

Using  such  an  hypothesis,  traffic  simulation  can  be  eliminated 

and  the  structural  analysis  program  need  only  be  utilized  to 

generate  a  minimum  set  of  single  axle  events,  for  example, 

at  increments  of  4000  lbs  in  weight  from  4000  lbs  to  40,000  lbs 

on  the  basis  of  a  speed  normalized  duration.   If  necessary, 

speeds  of  10  mph,  50  mph  and  90  mph  could  be  imposed  to  form 

an  interpolative  basis  for  speed,  rather  than  using  an  estimation 

from  the  single  normalized  set.   This  would  require  60  cases 

to  be  run  on  the  stress  program  for  each  bridge  to  be  analyzed, 

considering  two  lanes. 

On  the  basis  of  such  axle  weight  speed  stress  curves,  a 
composite  curve  may  be  formed  for  a  total  truck  from  the  set 
of  axle  stress  curves  which  represent  its  component  axles. 
These  individual  curves  become  essentially  additive  as  a 
function  of  the  axle  distance,  i.e.,  lag  time, 


TAB 

= 

V 
dAB+dBC 

TAC 

V 
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Hence,  the  second  axle  begins  to  effect  the  bridge  t__. 

......  Ac 

seconds  after  the  first  axle's  effect  occurs,  and  the 
third  axle  begins,  similarly,  t   seconds  after  the  first 
axle.   The  stress  imposed  at  tQ  +  xAB  is  then  the  reference 
for  axle  #2's  stress  curve,  that  is,  the  stress  curve  for 
the  second  axle  is  added  to  the  stress  at  time  tQ  +  t,r. 
Similarly,  the  stress  for  the  third  axle  is  added  to  the 
stress  remaining  on  the  bridge  beam  at  time  t.  +  t   seconds, 
that  is,  given  the  axle  curves: 


Axle   1,    W.,V 


Axle   2,    W9,    V 


i 


«wwlr 


/ 


/ 


02 


,VW\Vli*l 
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a     Axle  3,W3,V 


The  total  approximate  truck  stress  curve  is  then  formed  as 
a  composite  of  these  single  axle  curves,  i.e., 


I  I 

Axle  #1     Axle  #2 


Axle  #3 


This  eliminates  the  need  for  evaluating  each  truck's 
effect  on  a  bridge  through  the  structural  analysis  program, 
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Given  a  particular  truck,  i.e.,  distances  to  each  axle  from 
the  front  axle,  speed,  total  weight  and  fraction  of  weight 
at  each  axle,  a  composite  stress  curve  can  be  generated 
from  a  set  of  such  single  axle  stress  curves.   The  steps 
to  construct  such  a  composite  stress  curve  are  as  follows: 

1.   Select  four    single  axle  stress  curves,  for  each 
axle  of  the  given  truck,  i.e., 


o. . (t)  =  F(V. ,W. ,t) 

lj  1   j  ' 

a.  .  ..  (t)  =  F(V.  ,W.  .,  ,t) 
li+l         i   1+1 


ai+1.(t)  =  F(Vi+1/Wj,t) 
°i+l,j+l(t>  =  F<Vi+l'Wj+l't) 


where 


V.<V<V.,n,  W.<W<W.11,t  =  time 
l—  —  l+l    3  —  —  j+1 

and,  V  and  W  are  the  speed  and  weight  respectively  for  each 
of  the  given  truck's  axles. 

2.  Linearly  interpolate  for  a(V,W.,t)  and  o(V,W.  ,,t) 
from  the  selected  single  axle  stress  curves  for  each  axle 
of  the  given  truck. 

3.  Linearly  interpolate  between  a(V,W.,t)  and  a(V,W.  ,,t) 
for  a(V,W,t)  for  each  axle  of  the  given  truck. 

4.  Establish  a  base  stress  curve  with  aA(v,W,t)  for 
the  first  axle,  A. 
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5.   Calculate,  as  a  function  of  distance  between  the 
axles  and  the  given  truck's  speed,  the  lag  time  for  each 
rear  axle,  i.e., 

dAB 
Tfi  =  -y-   where  A  represents  the  front  axle,  B 

the  2nd  axle,  dAB  the  distance  between  A  and  B,  and  V  the 

speed  of  the  truck.   For  a  third  axle, 

'      ,   =  **£ 

C     V   * 

6.   The  second  axle's  stress  curve,  a_,(V,W,t)  is  then 
added  to  the  base  curve  a  (V,W,t)  translated  by  t  ,  i.e., 

"Truck  -  °A<V'W't)  +  oB(V,W,t+TB) 

Similarly,  the  third  axle,  if  it  exists,  is  added  to  the 
base  curve  to  form  the  final  truck  stress  curve,  i.e., 

°Truck  =  °A(V,W,t)+aB(V,W,t+xB)+a(V,W,Tc) 

Truck  platoons  can  further  be  approximated  in  the  same  manner 
as  the  individual  truck  stress  curves,  without  simulating  each 
possible  platoon  event,  in  a  similar  manner.   Lateral,  or  lane 
adjacent,  formed  platoons  can  be  approximated  by  considering 
the  load  imposed  on  each  beam  by  any  given  axle  load  in  the 
same  manner  as  described  above.   This  imposes  the  need  to  re- 
tain the  stress  curve  for  each  beam  in  a  bridge  for  each  of 
the  axle's  needed,  for  example,  the  60  at  4000#  increments. 
If  a  bridge  has  eight  beams  then  a  total  number  of  stress  curves 
required,  for  one  sample  point  per  beam,  in  the  "library"  of 
axle  stress  curves  for  that  bridge  would  be  480  curves.   This 
approach  provides  a  realistic  and  minimal  computation  method 
for  the  generation  of  unique  truck  platoon  events  on  a  bridge. 

The  statistics  of  truck  platoon  events,  per  unique  event, 

could  then  be  imposed  for  a  long  term  prediction  of  stress 

level  incidence. 
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For  further  simplicity  dual  axles  may  be  represented  as 
single  axles  with  the  combined  load  of  two  real  axles.   At 
88fps  (60  mph)  a  lag  of   ~. 05  sec  occurs  between  dual  axles 
at  4'  centers  and  at  44fps  (30  mph)  a  lag  of  ~.10  sec  occurs. 
Resolution  of  field  data  is  not  good  enough  to  completely 
verify  the  existence  of  a  composite  dual  axle  signature. 
Hence,  the  single  axle  signature  representation  must  vary, 
nominally  between  4000   and  40000#  rather  than  the  2000# 
and  30000#  which  would  normally  be  expected.   This  provides 
10  axle  signature  weights  at  constant  speed  for  a  single 
position  on  a  beam.   Two  nominal  speed  classes  of  30  mph 
(midpoint  of  10  to  50  mph  interval) and  70  mph  (midpoint  of 
50  to  90  mph  interval)  would  produce  a  total  of  20  stress 
signatures  per  beam  position  per  lane.   Since  the  structural 
program  is  designed  as  it  is,  it  will  simultaneously  produce 
the  stress  signatures  for  each  sample  point  on  each  beam. 

In  order  to  provide  realism  if  felt  necessary,  the  second  axle 
of  a  dual  set  of  axles  may  be  included  in  the  synthesis  of  the 
total  response  trace  in  the  same  manner  as  a  rear  axle 

as  a  function  of  speed. 

In  order  to  generate  the  stress  data  necessary  to  the 

method  described  in  the  foregoing,  a  simple  load  generator 

was  required.   It  must  output  load  data  that  is  consistent  with 

the  present  input  load  data  specification  of  the  structural 

analysis  program,  which  was  developed  for  BRIGLD1.   All  data 

is  sequential  as  shown  below,  and  in  binary  form,  on  magnetic 

tape, 

t  =  time  of  load  sampling  (sees) ,  a  single  value 

and  real 
At  =  time  increment  of  load  samples  (sees) ,  a  single 

value 
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NEV  =  Event  number.   This  number  will  also  act  as  a 
pointer  to  the  generated  stress  curve  for  a 
particular  parametric  triplet,  that  is, 
|v,W,NLf  ,  it  is  an  integer 

NAXLE  =  1  (integer) 

NTRUCKS  =  1  (integer) 

LAST  =  A  logical  indicator  where 

1.  True  implies  the  last  set  of  load  data 
for  a  given  event,  or 

2.  False  implies  all  load  data  sets  except 
the  last. 

Truck  Type  (20) :   Twenty  integer  values  are  furnished 

in  the  data  block.   However,  this  parameter 

has  no  meaning  for  this  program  and  all  values 

will  be  set  to  unity. 
Truck  Weight  (20) :   Twenty  real  values  are  furnished. 

However,  the  first  value  will  be  W  and  all  other 

values  will  be  set  to  zero. 
Truck  Speed  (20)  :   Twenty  real  values  are  furnished. 

However,  the  first  value  will  be  V,  and  all 

other  values  will  be  set  to  zero. 
Truck  Lane  (20):   Twenty  integer  values  are  furnished. 

However,  the  first  value  will  be  NNL  and  all 

other  values  will  be  set  to  zero. 
Truck  Entry  Time  (20) :   Twenty  real  values  are 

furnished,  all  zero. 
Axle  Position  (50) :   Fifty  real  values  are  furnished, 

the  first  value  is  X  and  all  other  values  are 

zero. 
Axle  Lane  (50) :   Fifty  integer  values  are  furnished, 

the  first  value  is  NNL  and  all  others  are  zero. 
Axle  Weight  (50) :   Fifty  real  values  are  furnished, 

the  first  value-  is  W  and  all  others  are  zero. 
Travel  Distance  (50) :   Fifty  real  values  are  furnished, 

the  first  is  AX  and  all  others  are  zero. 
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Acceleration  (50) :  Fifty  real  values  are  furnished, 
all  zero. 

In  addition  to  requiring  a  special  load  event  generator,  a 

minor  modification  of  the  structural  analysis  program  was 

also  required.   This  modification  allowed  the  generated 

stress  data  to  be  retained  in  an  auxiliary  storage  media,  magnetic 

tape,  for  processing  by  the  program  developed  in  this  portion 

of  the  investigation,  which  is  described  in  FHWA-RD-73-45  and 

referred  to  as  the  "Stress  Histogram"  program,  HISGEN. 

Consistent  with  the  present  form  of  the  output  of  the  structural 
analysis  program,  from  one  to  three  sample  points  per  beam  are 
allowed,  for  the  stress  trace  output  to  tape.   No  sophisticated 
processing  of  the  generated  digital  stress  trace  data  is 
performed  by  the  structural  program. 

The  following  data  is  passed,  as  a  header  record  for  each  stress 
event,  from  the  structural  program  to  the  histogram  program: 

NEV  =  the  event'  number,  integer 

W   =  the  axle  weight,  #,  real 

V   =  the  axle  velocity,  f ps ,  real 

NL  =  number  of  lanes,  integer 

Vt  =  time  increment  used  by  the  structural  program,  real 

B_   =  bridge  length,  ft,  real 

NP(J)=number  of  sample  points  per  Jth  beam,  integer 

NB  =  number  of  beams ,  integer 

NNL  =  Lane  Number 
After  a  header,  the  set  of  stresses,  at  a  given  time,  for  each 
sample  point,  1  to  3  per  beam,  are  output  as  a  record.   This 
is  repeated  for  all  time  points  contained  in  a  given  stress 
event.   The  next  stress  event  repeats  the  above  process,  starting 
with  a  new  header  record. 
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o{.Xir   1)  ,  a(x2,  1),   a(x3,  1) 
a(X1/  2),  a(X2,  2)  ,   a(X3,  2) 

•  •  • 

•  •  • 

a(x1,NB),  a(x2,NB),   a(X3,NB) 

where  a  (X  .  j)  is  defined  as  the  stress,  at  time  t,  for 
position  X  ,  1<  p<  3,  on  beam  j,  for  1<  j  <ND. 

Stress  Histogram  Program  (HISGEN) 

The  overall  program  developed  in  this  portion  of  the 
project  was  the  "Stress  Histogram  Program."   Its  function 
is  to 

1.  Construct  and  save,  on  magnetic  tape,  time 
dependent  digital  stress  traces  at  each  sample  point,  per 
beam,  for  each  load  event,  that  is,  each  prescribed  axle 
weight,  velocity  and  lane.   This  amounts  to  primarily  a 
sort  of  the  data  output  by  the  structural  analysis. 

2.  Construct  and  save,  on  magnetic  disk,  the  stress 
trace  of  each  truck  category  for  each  sample  point  on  each 
beam.   Each  truck  category  is  defined  in  terms  of  weight, 
velocity  and  lane  of  occupancy.   These  three  variables  must 
correlate  with  the  distribution  functions  prescribed  for  these 
variables  and  input  to  the  program,  as  discussed  later  under 
"Platoon  Construction." 

3.  As  a  function  of  the  given  distribution  functions 
of  truck  weight,  velocity,  lane  occupancy  and  platoon  size, 
construct  digital  platoon  stress  traces  for  each  sample  point, 
per  beam. 

4.  From  each  digital  platoon  stress  trace,  per  sample 
point,  extract  the  maximum  stress  values. 
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5.  Calculate  the  probability  of  existence  of  each 
platoon  construct  and,  as  a  function  of  ADT,  truck  fractions 
of  the  traffic  population,  the  predicted  growth/decay 
function,  and  the  life  span  desired  for  the  bridge,  estimate 
the  number  of  occurrences  of  each  platoon. 

6.  For  each  maxima  stress  value  extracted,  for  each 
platoon,  add  the  number  of  occurrences  to  each  stress  class 
interval  in  each  sample  point's  histogram  of  the  platoon. 

7.  Exhaust  the  set  of  platoons  defined  in  the  given 
distribution  data  and  output  the  histogram  data. 

This  program  has  been  completely  subroutined  in  order  to 
allow  ease  for  future  modifications.   From  a  logic, 
functional,  and  repetitive  use  point  of  view  there  was  no 
need  to  use  subroutines.   It  is  only  from  a  replacement 
or  upgrading  point  of  view  that  this  approach  was 
warranted. 

Truck  Trace  Construction 

The  formation  of  the  composite  axle  stress  trace  to  synthesize 
the  stress  trace  of  a  given  truck  was  accomplished  in  a  very 
straightforward  manner.   Linearity  of  response  is  assumed. 

The  delay  time  between  axles  was  estimated  as 
AAX 

where  AAX  is  the  distance  of  any  given  axle  from  the  front 
axle,  and  V  is  the  truck's  speed,  fps. 

The  integer  number  of  time  increments,  N,  contained  in  the 
delay  time  is  calculated  and  the   fractional  part  of  the 
delay  time,  that  is,  less  than  a  full  time  increment 
determined  from 
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At  =  T-N.At 

A. linear  interpolative  formula  was  used  to  approximate 
that  portion  of  the  axle  stresses  to  be  superimposed  on  the 
first  axle's  stress  trace,  that  is,  for 

t  z.t  <t„ 
m  —  p  mH 

V  V-tt    (°pH  "  °p> 


where  a     is  from  the  reference  on  base  stress  trace  and  a 
m  p 

is  from  the  stress  trace  of  the  axle  being  superimposed. 


The  above  process  is  performed  for  each  set  of  truck  parameters 
defined  in  the  input  data  at  each  sample  point  on  each  beam 
from  the  synthetically  generated  stress  data. 

The  efficient  retrieval  of  the  appropriate  point  trace 
data  for  each  axle  specified,  via  the  truck  parameters,  was 
a  larger  problem  than  forming  an  approximate  truck  stress 
trace.   This  is  discussed  in  the  next  subsection. 

Point  Trace  Retrieval  Algorithm 

To  determine  the  record  within  the  total  point  trace  file 
containing  a  trace  representative  of  an  axle  for  a  given 
truck,  with  speed,  weight  and  lane  occupancy  specified, 
the  event  number,  or  record  number,  is  determined  as  follows: 


NT   =   given  as  the  total  number  of  lanes 

w   -w 

N   =  jnax_o 
W       AW 

.  v   -v 

N   =   max   o 

V       AV     +  ± 

where  N  ,  w    ,    W  ,  AW,  V    ,  V   and  AV  are  the  same  values 
Jj   max    °       max   o 
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as  used  by  the  Synthetic  Load  Generator.   For  the 
given  truck 

V  =  speed  is  given 

W  =  weight  of  axle  of  interest  is  given,  and 

NNL  =  lane  of  occupancy  is  given 

The  record  number  NR  is  then  determined  from 
V-V 


N. 


i   _ 


V      AV 


NR0  =   NV  •  NW  •  NL 


N '       = 

W-W 
o 

AW 

NR1   = 

N'     .    N_ 
W           L 

where 


NB 

NPTS=   .2.   N  ,.x 
D=l    PCD) 


NR2  =  <NR0  +  NR1  +  NNL>  NPTS 

j-l 

NR3  =   NR2+  Z   <NP(P)>  +  i 
p=l 


NP(p)  is  the  number  of  sample  points  on  the  p   beam 


and  i  is  the  sample  point  index  at  the  sample  point  of 
interest  on  th  j   beam. 


Platoon  Construction 

In  order  to  apply  the  foregoing,  the  platoon  composition 
must  be  defined  and  be  constfuctable,  that  is,  in  terms  of 
constructing  its  truck  content  and  spatial  configuration,  its 
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stress  curve  and  the  incidence  of  occurrence.   A  basis  for 
providing  a  means  to  define  the  platoons  would  be  as 
follows,  as  an  example: 

1.  Assume  two  speed  class  intervals ,  as  indicated 
previously,  for  each  truck  type,  that  is,  lo-speed  (0-50) 
and  hi-speed  (50  and  above) , 

2.  Assume  nine  truck  configurations, 

3.  Assume  spatial  permutations  have  no  effect,  since 
stresses  are  in  the  elastic. range  and  are  linear  and  additive, 
and  that  a  worst  case  platoon  positional  configuration  will 
be  sufficiently  representative  of  anticipated  stresses  for 
spatial  permutations  of  the  given  platoon  components, 

4.  Assume  single  weight  for  each  truck  category, 
that  is,  maximum  anticipated, 

5.  Assume  a  maximum  platoon  size  of  six  trucks, 

This  constrains  the  set  of  platoon  elements  to  18  and  the 
total  number  of  load  events  is  determined  as: 


Platoon  Size 
1 
2 
3 
4 
5 
6 


Number  of 

Different  Platoons 

18 

153 

816 

3060 

8568 

18564 

31179 


Truck 

Population 

18 

306 

2448 

12240 

42840 

111384 

169236 

This  implies  31179  unique  load  events  involving  a  truck 
population  of  169236. 
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The  criteria  utilized  to  form  the  platoon  spatial  configura- 
tion should  be  realistic,  for  example, 

1.  Slowest  vehicle  (s)  in  the  right  lane 

2.  Minimum  headway  distance  between  vehicles 

For  this  case,  it  then  appears  that  31179  platoons  would 
need  to  be  evaluated,  that  is,  their  stress  effect  determined 
and  the  incidence  of  each  platoon  determined. 

The  construction  of  the  estimated  stress  response  of  each 
beam  would  be  determined  as  indicated  previously  for  trucks. 
Multiple  lane  occupancy  is  taken  care  of  by  the  addition  of 
each  sample  point's  response  for  each  truck  due  to  lateral 
effect.   Time  lag  effects  are  included  in  the  same  manner 
as  shown  previously  for  truck  trace  construction. 

In  order  to  determine  the  incidence  of  each  platoon  configura- 
tion's stress  response,  certain  statistics  regarding  the  traffic 
must  be  available  or  estimated  from  best  available  data, 
that  is, 

1.  Average  daily  traffic  (ADT) 

2.  Percent  of  traffic  that  is  truck  traffic,  P™ 

3.  Distribution  of  significant  truck  types  in  total 
truck  population, 

[PT1/  PT2'  PT3'  PT4'  PT5'  PT6*  PT7'  pTs'  PT9^ 

9 

where.!.  P_.  -  1,  for  this  example. 

1=1  Ti    '     „        r 

4.  Distribution  of  velocity  population  for  each  type, 
for  example,  for  two  class  intervals, 
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1  (PV1/PV2}T1'  (PV1/PV2)T2/  CPV1/PV2)t3'  CPyl,PV2)T4/ 

(PV1 '  PV2 }  T5  '  (PV1 '  PV2 J  T6  '  (PVl '  PV2 )  T7  '  (PV1 '  PV2 )  T8  ' 

tPVl'PV2)T9] 
2 

where  I         ^v^Ti  =  1* 
i=i        J 

5.  Distribution  of  platoon  population  by  configura- 
tion, that  is,  trailing  or  passing 

jPNLl'  PNL2(' 

2 

Where  E    PXTT  .  =  1 
.  ,    NLi 
i=l 

6.  Distribution  of  platoon  size  in  total  platoon 
population 


[P1'  P2'  P3'  P4'  P5'  P6J 

6 
where  Z'       P.  =  1. 
i=i    x 

7.   A  growth  or  decay  function  for  long  term  growth 
or  decline.   An  arbitrary  function  can  be  established  to 
provide  a  mathematical  basis  for  this  function.   However, 
a  more  realistic  approach  is  to  use,  for  now,  a  linear 
approximation  between  current  estimated  ADT  and  the  designed 
daily  service  volume  for  the  relevant  section  of  highway, 
at  the  estimated  point  in  time  that  the  highway  will  reach 
its  designed  daily  service  volume.   From  that  point  assume 
the  volume  to  remain  constant  for  the  life  span  of  the 
bridge,  for  example, 
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Designed  Service  Volume  of  Adjacent  Highway 


i 

i  Time 

t(yearsj 


The  total  estimated  traffic  is  then 

(t  -t  ) 
Tr  =  365  [-^V2-  (ADTq  +  ADT^)  +  (t^) AD^] 


where  t  =  beginning  of  time  in  years 

point ,  in  years,  service  volume  of 

highway  reached  (estimated) 

total,  in  years,  service  life  of  bridge 

desired 

average  daily  traffic  initially 

average  daily  traffic  at  design  service 

volume  of  highway 


t,  = 


t„  = 


ADT   = 


ADT 


SV 


Lacking  statistical  data  and  a  rational  basis  for  predicting 
the  distribution  of  a  given  truck  type  over  the  population 
of  each  platoon's  population,  it  will  be  assumed  that  the 
distribution  of  truck  types  over  each  platoon's  population 
is  the  same  as  over  the  total  truck  population. 

If  the  total  traffic  for  the  period  of  interest  is  T   and  the 

r  r 

percent  of  trucks  in  that  population  is  PT  then  the  total 
estimated  truck  population  over  the  period  of  interest,  T  , 
is 

T  =  Total  Truck  Population  =  P   .  T  . 

The  total  truck  population  in  each  platoon's  population  is  then 
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m      _   m        p 
PI      -P   *   ^1 


m      _   m        p 
P2      XP   *   ^2 


TP3    Tp  *  P3 


TP4  =  TP  •  P4 


T    =  T     P 
P5    1P  *  *5 


m      _   m        p 
P6        P   '     6 


The  total  of  each  type  truck  in  each  platoon's  population  is 


TPlj  =  Tpl  .  PTj  for  1<  j  49,  j  =  type 


m  m  p 

1P2j  1P2     *     *TJ 


TP3j         TP3    '    PTj       " 


rp  _      rn  p  H  II  II 

P4j         1P4    *       Tj 

m  _      m  p  II  II  II 

1P5j  1P5    *       Tj 

rp  _      m      C  p  II  II  II 

iP6j  P       *       Tj 

The  total  of  each  velocity  class  of  each  type  truck  in 
each  platoon's  population  is: 

T^,  .,  =  T^,  .  .  P,rl  for  1  <  k  <  2,  k   =  velocity  class 
Pljk    Plj     Vk       —  —  '  2 

m  _   rp  p        II  M  ll  II 

1P2jk    iP2j  *  *Vk 


m  _   m  p 

P3jk    P3j  *   Vk 
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T_.  .,     =   T^.  .    .    PTT1       for   l<k<2,    k     =      velocity   class 
P4jk  P4j  Vk  —     —  . 

m  _      m  p  II  II  II  II 

P5jk  P5j    *      Vk 

m  _      m  p  II  II  II  II 

P6jk  P6j    *      Vk 

The  determination  of  the  incidence  of  single  truck  events 
is  straight  forward,  that  is,  the  calculation  of  Tp, .,  given 
the  type;  and  velocity  class  k. 

While  the  above  provides  the  incidence  of  each  j  type  and  k 
velocity  class  in  each  platoon's  population,  it  does  not 
provide  for  the  incidence  of  each  platoon's  configuration 
of  size  greater  than  one,  that  is,  given  a  particular  set 
of  trucks  what  is  the  incidence  of  that  specific  platoon. 

For  two  truck  events  the  incidence  of  a  specific  event, 
that  is,  for  a  truck  of  type  p  and  velocity  class  a,  and 
for  a  truck  of  type  q  and  velocity  class  b,  whose  composite 
stress  curve  has  just  been  generated,  is  estimated  as 

=  (Tp2Pa)  T 
XP2  (p,a)  (q,b)    l  Tpa  ;  ±P2qb 

This  value  then  becomes  the  incidence  of  the  platoon  stress 
traces  found  in  the  composite  stress  curve  for  the  (p,a)  (q,b) , 
2  component,  platoon.   It  is  assumed  that  each  element,  that 
is,  type  and  velocity,  will  pair  with  itself  to  form  a  2 
element  platoon. 

It  becomes  easier  to  work  with  the  probabilities  for  platoons 
of  size  3  and  larger  than  with  raw  incidence  values,  for 
example, 

T 
P3jk  =  "T^1  for  a11  j  and  k 
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T 
=  ■  !/3k  for  all  j  and  k 


4jk    T 


P4 


p    =  TP5jk 
P5jk     Tp5 

P6jk     Tp6 


Utilizing  the  joint  probability  of  occurrence  for  determining 
the  probability  that  three  different  trucks  will  occur  in  a 
single  event,  then 


P3Cp/a) (q,b) (r,c)    P3pa  '  P3qb  *  P3rc 


Hence,  given  a  composite  3  truck  platoon  stress  response 
curve  for  trucks  defined  as 

1 .  j  =  p  and  k  =  a 

2 .  j  =  q  and  k  =  b 

3 .  j  =  r  and  k  =  c 

then  the  incidence  of  this  platoon  event  becomes 

T3(p,a)  (q,b)  (r,c)    TP3  (P3pa  *  P3qb  *  P3rc} 
Similarly,  for  a  four  element  platoon: 

1 .  j  =  p  and  k  =  a 

2 .  j  =  q  and  k  =  b 

3 .  j  =  r  and  k  =  c 

4 .  j  =  s  and  k  =  d 

and     T4(P/a) (q,b) (r,c) (s,d)    TP4(P4pa*  P4qb*  P4rc*  P4sd) 

for  a  five  element  platoon: 

1.  j  =  p  and  k  =  a 

2 .  j  =  q  and  k  =  b 

3 .  j  =  r  and  k  =  c 

4 .  j  =  s  and  k  =  d 

5 .  j  =  t  and  k  =  e 


and 


■5(p,a)  (q,b)  (r,c)  (s,d)  (t,e)    TP5(P5pa*  P5qb'  P5rc* 
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P5sd*  P5te}/ 

and  for  a  six  element  platoon: 

1 .  j  =  p  and  k  =  a 

2 .  j  =  q  and  k  =  b 

3 .  j  =  r  and  k  =  c 

4 .  j  =  s  and  k  =  d 

5 .  j  =  t  and  k  =  e 

6 .  j  =  u  and  k  =  f 
and 

T6(p,a)  (q,b)  (r,c)  (s,d)  Ct,e)  (u,f)=  TP6 (P6pa' P6qb*Prc * 

P6sd*P6te,P6uf} 
Stress  Range  Determination 

The  stress  maxima  are  then  determined,  for  each  platoon's 
stress  trace.   The  proper  class  interval  for  each  sample 
point's  histogram  are  then  determined  for  each  maximum 
stress  point  and  the  platoon's  incidence  is  then  added  into 
each  of  the  appropriate  class  intervals.   The  entire  process 
is  repeated  until  all  platoon  parameters  have  been  varied 
as  specified  in  the  input  data  and  all  combinations  of  these 
have  been  exhausted. 

Hence,  the  foregoing  provides  a  synthetic  means  of  generating 
long  term  stress  histograms  in  a  reasonably  efficient  manner 
and  on  a  sound  structural  analysis  basis.   Obviously,  as  such 
an  approach  is  utilized  increased  learning  will  occur  and 
improvements  in  it  can  be  made,  for  example,  increased  refine- 
ment of  truck  categories  and  platoon  definitions,  if  economics 
and  efficiency  allow  it. 

The  implemented  frequency  formulae,  to  determine  the  incidence 
of  a  given  platoon,  are  as  follows: 
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TT 


N 

I  nPp(n) 
n=l  F 

where  T  is  the  total  truck  traffic  as  defined  earlier, 

N  is  the  largest  platoon  size,  and  Pp(n)  is  the  probability 

of  occurrence  of  a  platoon  of  size  n  in  the  total  platoon 

population.   This  distribution  must  be  furnished  from  field 

data. 


The  probability  of  occurrence  of  a  truck,  out  of  the  total 
truck  population,  in  a  platoon  of  size  n  is 

n  PD(n) 
PT  (n)  = 


N 

I   nPp(n) 
n=l  * 

Given  the  probability  of  existence  for  each  truck  of  a 
platoon  of  size  n,  based  on  type/weight  and  velocity,  as 
previously  discussed,  that  is, 


\P1'  P2>  *  *  * '  Pn  } 


The  probability  of  each  element  of  the  set,  truck  in  the 
platoon,  existing  in  a  platoon  of  size  n,  out  of  the  total 
truck  population  is  determined  from 


Pln  -  Pi'Vn) 


P2n  =  P2.PT(n) 


P   =  P  ,P_(n) 
nn    n   T 


The  probability  of  all  the  specified  trucks  coexisting  in  a 
platoon  of  size  n  is 


gn    In   2n       nn 
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where  g  is  the  gth  set. 

and  the  incidence  of  occurrence  is  then 


I   =  P   .T 
gn     gn   p 


Verification  of  the  correctness  of  the  calculations  can  be 
made  by  comparing  the  given  distributions,  Pp  ,  ,,    against 
the  calculated  occurrences,  that  is, 

Q 

II    =  T_.P_(n) 
,  gn     P   P 
g=l 


The  class  intervals  of  the  sample  points  of  interest 
containing  the  stress  maxima  from  the  composite  traces,  for 
each  of  the  sample  points,  are  then  incremented  by  I   . 

To  retrieve  the  individual  truck  stress  traces,  for  forming 
the  platoon  stress  traces,  which  are  a  function  of  the 
sample  point  of  interest,  beam  of  interest,  lane  of 
occupancy,  speed  and  weight,  the  record  number  of  the  trace 
is  determined  from 


NR  =  i  +  NPS  +  NPTS  J[(kT-D  DNv+(kv-l)]NL  +  (^-D) 


where 

i  =  sample  point  index  on  given  beam,  j 

j-l 

NPS  =    I    NP(p) 

p=l 

NB 
NPTS  =    I         NP(j) 
j  =  l 

KL  is  the  given  lane  index, 
KV  is  the  given  speed  index, 
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KT  is  the  given  truck  category  index 
N  =  total  number  of  velocity  classes 
NT  =  total  number  of  lanes 
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Application  of  the  Approach 

In  order  to  apply  the  foregoing  approach,  it  is  necessary 
to  utilize  the  three  previously  described  computer  programs, 
i.e., 

1.  The  synthetic  load  generator  (SYNGEN)  to  generate 
single  axle  load  events, 

2.  The  dynamic  stress  analysis  program  (BRGSTRS)  to 
generate  single  axle  stress  signatures,  and 

3.  The  stress  histogram  program  (HISGEN)  to  calculate 
platoon  incidences,  platoon  stress  ranges  and  construct  the 
stress  histograms. 

While  BRGSTRS  had  been  tested  and  used  on  a  production  basis, 
as  described  in  the  next  section,  SYNGEN  and  HISGEN  had  only 
been  tested  with  debug  drivers  and  test  data.   A  realistic 
test  was  performed  using  the  truck  traffic  and  bridge  data 
for  the  southbound  span  of  the  bridge  at  1-83  and  Bunker 
Hill  Road  in  Maryland  (3).   This  same  bridge  was  also 
evaluated  using  BRIGLDl,  and  is  described  in  the  next  section 
of  this  report. 

An  interesting  comparison  of  computer  time  was  one  of  the  re- 
sults of  this  testing.   In  using  BRIGLDl  and  BRGSTRS  in 
combination,  a  five  hour  simulation  of  real  time  was  performed, 
generating  a  total  truck  population,  and  their  related  stress 
time  histories,  of  117  trucks.   The  computer  time  used  in  the 
"go"  step  for  each  program  was 

1.  For  BRIGLDl  approximately  .5  minutes,  and 

2.  For  BRGSTRS  approximately  11.0  minutes. 
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The  total  time  required  was,  then,  approximately  11.5  minutes 
using  the  traffic  simulation  approach  to  achieve  5  hours  of 
simulated  time  and  117  trucks. 

In  using  SYNGEN,  BRGSTRS  and  HISGEN  in  combination,  a  fifty 
year  life  span  period  was  used.   For  this  approach  the  period 
of  simulated  real  time  has  no  affect  upon  the  use  of  com- 
puting time.   The  period  can  be  5  minutes  or  a  thousand  years 
and  will  require  the  same  amount  of  computer  time.   The 
portion  of  this  approach  which  affects  computer  useage  is 
the  generation  of  the  single  axle  stress  traces  by  BRGSTRS. 
The  finer  the  mesh  of  speed  and  weight  used,  the  more  com- 
puting time  is  used.   However,  for  this  test  case  the 
computer  time  used  in  the  "go"  step  for  each  program  was 

1.  For  SYNGEN  approximately    .5  minutes, 

2.  For  BRGSTRS  approximately  8.1  minutes,  and 

3.  For  HISGEN  approximately   1.6  minutes. 

The  total  time  required  was,  then,  approximately  10.2  minutes 
using  the  stress  signature  approach.   This  generated  stress 
ranges  on  9  sample  points  on  5  beams  for  the  equivalent  of 
19399750  trucks  over  fifty  years.   In  this  instance  66  single 
axle  stress  traces  were  generated  by  BRGSTRS.   Eleven  weight 
values,  2  lanes,  and  three  speeds  were  used  by  SYNGEN  to 
create  the  6  6  load  events  for  BRGSTRS,  i.e., 

58.67  fps  <  V  £  102.67  varied  at  22  fps,  and 
2000  lbs  <  W  <  82000  lbs  varied  at  8000  lbs. 

The  same  bridge  structural  data  was  used  in  this  case  to 
generate  the  single  axle  stress  signatures  as  was  used  in  the 
traffic  simulation  case.   Nine  sample  points  on  beams  1,  2, 
3,  4  and  5  were  used  to  collect  the  stress  data. 
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Instead  of  using  weight  distribution  functions  per  truck  type, 
HISGEN  uses  only  the  maximum  anticipated  weight  per  type  to 
predict  on  a  worst  case  basis.   The  maximum  weight  and  the 
distribution  by  axle  per  type  was  the  same  as  used  in  the 
BRIGLD1  case.   The  distribution  by  type  over  the  total  truck, 
population  and  the  speed  distribution  within  each  type  was 
also  the  same  as  was  used  in  the  BRIGLD1  case.   However,  the 
speed  distribution  is  defined  in  a  different  manner  in 
HISGEN,  i.e.,  by  specific  speeds  and  the  probability  of  each 
speed.   The  mid-points  of  each  class  interval  defined  in 
BRIGLD1  were  chosen  as  the  specific  speeds  and  since  there 
were  ten  class  intervals,  a  probability  of  .1  was  assigned 
to  each  speed. 

Platoon  definition  was  restricted  to  size  1  and  2  platoons 
with  .80  and  .20  probabilities  respectively.   This  was  based 
upon 

1.  The  low  incidence  of  multi-truck  platoons  generated 
in  the  traffic  simulation, 

2.  The  uncertainty  about  the  validity  of  the  traffic 
simulated  data,- 

3.  No  evidence  of  multi-truck  platoons  in  the  field 
collected  data,  and 

4.  Intuitive  estimating  that  at  least  platoons  of  size 

2  would  occur.   However,  the  20%  probability  was  probably  too 
high  and  should  have  been  closer  to  1%. 

The  truck  traffic  load  data  used  was  based  upon  actual  field 
samples  but  no  long  term  growth  or  decay  was  used.   This  was 
primarily  due  to  the  desire  to  compare  the  stress  range 
histograms  generated  by  this  approach  with  the  field  collected 
sample  and  the  traffic  simulation  generated  sample.   Passing 
and  trailing  truck  platoon  formations  were  set  at  10%  and  90% 
respectively . 
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The  results  of  this  test  were  unverifiable  because  of  the 
stipulation  of  platoons.   Consequently,  a  second  case  was 
run  using  only  single  truck  platoons  with  a  10%  passing 
lane  occupancy.   This,  again,  was  unverifiable  but  only 
transferred  the  stress  maxima  between  beams.   The  results 
of  both  cases  are  shown  in  table  21.   The  incidence  used 
for  each  class  interval  of  stress  range  were  the  accumulated 
stress  peaks  for  the  nine  sample  points  on  five  beams. 
Consequently,  there  is  a  multiplicity  of  stress  peaks  for  a 
given  live  load.   For  the  case  with  platoons  of  size  1,  a 

L 

7 


7 
total  accumulation  of  23.59671  x  10   maxima  were  accumulated 


over  the  nine  sample  points  for  a  total  of  1.9399750  x  10 

truck  loads.   Using  a  very  simplistic  analysis,  the  9  sample 

7 
points  averaged  2.6218562  x  10   significant  stress  maxima. 

This  implies  that  significant  stress  maxima  were  distributed 

to  an  average  of  1.35  beams  per  truck  load. 


Significant  use,  manipulation  and  evaluation  of  this  metho- 
dology is  required  before  any  significant  conclusions  can 
be  drawn.   Applied  utilization  of  this  approach  during  the 
study  was  very  small  and  insufficient  to  provide  an  adequate 
learning  basis.   Results  from  methods  such  as  this  are  ex- 
tremely dependent  upon  how  a  user  applies  them  and  how  well 
he  understands  the  approach.   Futher,  as  more  field  data  is 
collected  and  more  learned  about  the  statistics  of  platoon- 
ing  trucks  and  the  configurations  of  the  platoons,  the  use 
of  this  method  will  improve. 
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GENERATED  HISTOGRAMS  AND  RESULTS 

As  a  result  of  the  review  of  a  suggested  set  of  criteria  for 
the  bridges  to  be  used  for  generating  stress  histograms,  a 
candidate  set  of  bridges  were  defined.   The  suggested  criteria 
were 

1.  Composite  deck,  reinforced  concrete,  with  steel 
main  girders  (wide  flange) . 

2.  Concrete  deck  with  box  main  beams. 

3.  Reinforced  concrete  deck  with  concrete  main  girders. 

4.  Three  cases  to  be  run  with  a  simple  span  configuration 

5.  One  case  to  be  run  with  a  3  span  continuous 
beam  (.8L-L-.8L) . 

6.  One  case  to  be  run  with  a  5  span  continuous 
beam  ( . 8L-L-L-L- . 8L) . 

The  identified  candidate  bridges,  which  were  assumed  to  have 
adequate  data  available  were 

1.  One  simple  span  on  Md.  301. 

2.  One  simple  span  in  Virginia. 

3.  The  Dumfries  bridge. 

4.  A  continuous  beam  in  Virginia. 

5.  A  continuous  beam  in  Md. 

6.  A  continuous  beam  in  Connecticut. 

7.  A  continuous  beam  in  Minnesota. 
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The  bridges,  and  their  related  parameters,  that  were  finally 
chosen  on  the  basis  of  acquiring  adequate  data,  within  time 
and  project  constraints,  were 

1.  1   simple  span  of  a  two  span,  both  identical,  bridge 
at  U.S.  301  and  Md.  Rt .  5,  fabricated  of  steel  beams  and  rein- 
forced concrete  deck. 

2.  The  north  span  of  a  bridge  at  U.S.  301  at  Western 
Branch  in  Maryland.   This  is  a  continuous  beam  of  three  spans. 

3.  The  fifth  span  of  a  5  span  simply  supported  bridge 
at  1-495  and  U.S.  1  in  Maryland. 

4.  The  center  span  of  a  3  span  simply  supported  bridge, 
southbound  on  1-83  at  Bunker  Hill  Road  in  Maryland. 

5.  The  northern  span  of  the  simply  supported  bridge  on 
1-95  at  Dumfries,  Virginia. 

Upon  detailed  investigation  of  the  "U.S.  301"  data  originally 
furnished  in  BRIGLDl,  as  contained  data,  it  was  determined 
that  the  truck  incidence  was  approximately  a  factor  of  twenty 
time  that  collected  in  actual  samples  for  the  two  bridges  on 
U.S.  301  used  in  this  project.   In  general  it  was  found  that  this 
data  was  unrealistic  and  totally  dependent  upon  the  "Headway" 
data  used. 

Additionally,  the  truck  data  contained  in  this  same  source 
data  was  also  unuseable.   Consequently,  each  of  the  five  bridges 
evaluated  required  a  completely  new  set  of  traffic  data,  as 
opposed  to  the  requirements  of  the  contract  for  this  project. 
The  resulting  load  simulations  evidenced  a  far  greater  time 
compression,  approximately  720  to  1,  of  real  time  than  was 
evidenced  during  the  sensitivity  testing.   This  was  almost 
totally  due  to  the  change  in  traffic  distribution  data. 
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The  truck  weight  and  length  data  contained  in  the  original 
BRIGLD1  program  was  modified  to  conform  to  more  recent 
findings  (2) . 

The  simulation  period,  for  the  purposes  of  this  project,  was 
selected  at  five  hours.   This  was  predicated  on  the  in- 
dications resulting  from  the  Sensitivity  Analysis. 

No  arbitrary  platoon  distribution  was  imposed  on  the  simulations. 

The  procedure  followed  in  producing  the  stress  range  incidence 
data  was 

1.  Development  of  the  necessary  truck  traffic  distribu- 
tion data  and  other  required  input  data  for  the  load  simulation 
program,  BRIGLDl  for  each  bridge  to  be  analyzed. 

2.  Preparation  of  the  data  for  input  to  the  computer  and 
running  of  the  load  simulation  program  on  the  computer  for  each 
bridge.   A  separate  magnetic  tape  was  utilized  for  the  generated 
load  data  for  each  bridge. 

3.  Analyzing  the  results  of  the  generated  load  data 
against  the  available  real  load  data  on  each  bridge. 

4.  Development  of  the  structural  input  data  for  each  bridge 
to  be  analyzed  as  required  by  the  dynamic  structural  analysis 
program,  BRGSTRS. 

5.  Preparation  of  card  input  data  for  each  bridge,  for 
BRGSTRS,  and  running  each  case  on  the  computer  using  BRGSTRS,  the 
card  input  data  and  the  bridge  corresponding  load  data  generated 
by  BRIGLDl  on  magnetic  tape. 

6.   Analyzing  the  generated  stress  range  data  against 
the  available  stress  range  data  on  each  bridge. 
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The  following  subsections  describe  the  data  used  for  each  bridge 
and  the  results  obtained,  in  comparison  with  available  real 
field  data. 

U.S.  301  and  Md.  Rt.  5 


Auto  traffic  was  assumed  to  be  an  arbitrary  constant  and  set 
at  83%  of  the  total  traffic  for  this  bridge.   The  truck  dis- 
tribution data  used  was  based  upon  actual  measurements  taken 
at  the  bridge  (3) .   This  data  only  allowed  for  five  truck 
types,  i.e.,  2D,  3,  2S1,  2S2  and  3S2.   A  comparison  of  the 
measured  incidence  of  these  truck  types  and  the  five  hours  of 
simulation  generated  by  BRIGLD1  is  shown  in  Table  22. 

In  this  case  no  means  of  comparing  total  or  average  truck 
weights  directly  was  possible  from  available  data. 

The  manner  in  which  the  truck  traffic  appeared  on  the  bridge 
deck  in  the  simulation  generated  84  load  units  of  84  single  truck 
events  in  the  five  hours  of  simulated  loading.   The  real  sample 
contained  5284  trucks  in  a  7  day  period  or  one  truck  every  1.91 
minutes.   The  simulated  sample  contained  approximately  1/2  of 
the  truck  traffic  actually  measured.   This  implies  that  the 
headway  distribution  data  used  in  BRIGLDl  was  skewed  to  the  low 
side.   This  is  possible  since  it  was  generated  from  2  and  4 
hour  averages  provided  by  the  real  data.   The  high  frequency 
or  short  interval  data  was  masked  by  the  longer  term  means. 
This  variance  could  have  been  corrected  in  BRIGLDl  by  adjusting 
the  headway  tables  to  provide  for  the  proper  mean  rate. 

The  input  data  used  to  generate  this  loading  data  is  shown  in 
Figure  6  as  is  the  five  hour  generated  data.   A  magnetic  tape 
containing  the  synthetically  generated  time  dependent  load 
data,  was  produced  simultaneously  for  input  to  the  dynamic 
bridge  stress  analysis  program  (BRGSTRS) » 
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TABLE  22.   Comparison  of  Truck  Distributions  for  the  U.S.  301 
and  Md.  Rt.  5  Bridge 


Truck  Types 

2D 

3 

2S1 

2S2 

3S2 

Measured  (%) 

22.5 

5.7 

8.5 

20.0 

43.3 

Simulated  (%) 

22.6     5.9 

5.9 

22.6 

42.9 
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Card  input  data  was  prepared  from  available  information 
on  the  bridge,  in  particular  the  south  span,  and  where  data 
was  insufficient  ,  standard  interstate  bridge  data  was 
utilized. 
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This  particular  case  was  the  first  production  run  attempted 
on  BRGSTRS  and  there  were  the  usual  difficulties  in  learning 
how  to  use  a  newly  constructed  program  and  debugging  the 
input  data.   The  data  interface  from  BRIGLD1  worked  perfectly. 
However,  it  was  a  well  planned  move  to  limit  the  simulation 
period  to  five  hours.   The  dynamic  stress  analysis  program 
required  45  minutes  of  CPU  time  (360/65)  in  the  "go"  step. 
This  implies  a  6.67  to  1  compression  of  real  time  by  the 
structural  analysis  program.   However,  this  is  a  direct  func- 
tion of  truck  density  and  cannot  be  relied  upon  as  a  rule. 
This  case,  as  indicated  earlier,  simulated  84  single  trucks 
in  the  five  hour  period  or  one  truck  every  3.57  minutes. 
This  span  was  also  a  simple  span  which  did  not  require  re- 
petitive inversion  of  the  modified  flexibility  matrix. 

The  gross  results  of  the  generated  stress  data  for  the  simula- 
ted truck  traffic  over  the  five  hour  period  is  shown  in 
comparison  with  field  measured  data  on  this  span  in  Table  23. 

The  results  compare  very  favorably  with  the  data  measured 
at  the  center  span  position  on  the  cover  plate.   This  was  for 
the  3rd  beam  in  both  cases.   The  input  data  defined  a  cover 
plate  on  the  beam  and  the  resultant  calculations  should 
correlate  with  the  cover  plate  data.   The  variances  in  the 
400  -  600  and  800  -  1000  class  intervals  are  apparently  due  to 
a  variance  in  the  weights  of  the  trucks  of  the  same  types 
between  the  real  and  the  simulated.   As  indicated  earlier,  there 
was  no  direct  means  of  correlating  the  weights  of  the  trucks 
within  a  given  type.   The  traffic  simulator  program,  BRIGLD1, 
appears  to  have  generated  heavier  trucks,  by  type,  than  the 
real  sample.   This  tendency  was  verified  in  the  fifth  test  case, 
the  Dumfries,  Virginia,  bridge,  which  is  described  later.   The 
truck  weight  distribution  data  utilized  was  based  upon  a 
generalized  data  covering  several  of  the  mid-Atlantic  states 
(2) .   The  local  utilization  of  these  truck  types  in  the  area 
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surrounding  this  bridge  may  be  toward  lighter  loading  than 
that  found  in  the  total  mid-Atlantic  region. 

The  input  data  used  to  describe  the  structual  characteristics 
of  this  bridge  for  BRGSTRS  is  shown  in  Figure  7.   A  sample 
of  the  stresses  calculated  are  also  shown  in  Figure  7  as  is 
an  example  of  the  graphical  output.   The  scaling  on  the 
graphical  output  was  changed  in  the  later  cases  to  one  half 
that  used  in  this  case.   The  histogram  generated  for  this  case 
is  also  shown  in  Figure  7.   As  can  be  seen,  the  selection  of 
the  class  interval  size  and  maximum  range  was  poorly  defined 
in  the  input  data  to  the  program.   This  input  was  modified 
in  the  other  cases. 
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U.S.  ,301  and  Western  Branch 

Auto  traffic  was  assumed  to  be  an  arbitrary  constant  and  set 
at  83%  of  the  total  traffic  for  this  bridge.   The  truck  dis- 
tribution data  used  was  based  upon  actual  measurements  taken 
at  the  bridge  (4).   This  data  only  allowed  for  five  truck 
types,  i.e.,  2D,  3,  2S1,  2S2  and  3S2.   A  comparison  of  the 
measured  incidence  of  these  truck  types  and  the  five  hours  of 
simulation  generated  by  BRIGLDl  is  shown  in  Table  24. 

In  this  case  no  means  of  comparing  total  or  average  truck 
weights  directly  was  possible  from  available  data. 

The  manner  in  which  the  truck  traffic  appeared  on  the  bridge 
deck  in  the  simulation  generated  65  load  events  of  65  single 
truck  events  in  the  five  hours  of  simulated  loading.   The 
output  from  BRIGLDl  of  this  data  is  shown  in  Figure  8. 
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Table  24.   Comparison  of  Truck  Distribution  for  the  U.S.  301/ 
Western  Branch  Bridge 


Truck  Types 

, 

2D 

■  ' 

3 

2S1 

2S2 

3S2 

Measured  (%) 

27.0 

6.0 

5.0 

21.0 

41.0 

Simulated  (%) 

30.8 

4.6 

0.0 

20.0 

44.6 
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As  in  the  previous  case,  the  truck  traffic  generated  by 
BRIGI^Dl  was  less  than  the  real  sample.   However,  in  this 
case  the  real  sample  was  taken  over  a  period  of  four  days 
in  nominally  4  hour  periods  during  daytime  and  may  not  be 
truly  representative  of  long  term  behavior.   The  simulator 
generated  at  the  rate  of  13  trucks  per  hour  while  the  real 
sample  evidenced  a  rate  of  20  per  hour.   This  case  actually 
came  closer  to  generating  the  proper  amount  of  truck  traffic 
for  the  five  hour  simulation  period  than  any  of  the  other 
cases  at  65%  of  the  real  sample's  rate.   This  undergeneration 
is  again  due  to  the  unavailability  of  the  high  frequency  truck 
occurrence  data  which  would  have  allowed  development  of  a 
better  set  of  headway  input  data.   This  can  be  easily  adjusted 
by  shifting  the  headway  distribution  function  into  the  shorter 
time  interval  region. 

In  terms  of  the  stress  range  portion  of  the  analysis ,of  this 
bridge,  this  was  the  only  continuous  multi-span  beam  bridge 
tested  on  the  program  BRGSTRS.   All  of  the  other  cases  were 
of  a  simple  span  design. 

The  run  with  BRGSTRS  was  terminated  due  to  excessive  CP  time. 
Eight  load  events  were  processed  in  58  minutes  of  CP  time. 
These  eight  load  events  were  equivalent  to  approximately  40 
minutes  of  simulated  times.   The  extremely  slow  calculation 
of  the  stresses  is  due  to  the  repetitive  solution  of  the 
flexibility  matrices  and  the  inversion  of  their  modified  form. 
However,  the  calculated  dynamic  stresses  provided  13  points  of 
maxima  from  4  specified  sample  points,  i.e.,  one  sample  point  on 
beams  2,  3,  4,  and  5.   It  was  decided  that  due  to  the  extreme 
slowness  of  BRGSTRS  in  handling  a  continuous  beam  case  that  the 
results  obtained  from  this  one  run  would  be  used  for  comparative 
analyses  with  real  data  to  provide  some  insight  into  the  validity 
of  the  BRGSTRS  output  for  continuous  beams.   This  comparison 
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is  shown  in  Table  25.   The  predicted  data  was  based  upon  an 
extremely  small  sample.   Further,  the  load  data  generated 
by  BRIGLD1  during  the  first  half  hour  of  simulated  time  was 
skewed  toward  low  truck  weights.   The  heaviest  truck  generated 
during  this  period  was  in  the  32,000  to  40,000  lb.  range  and 
it  caused  a  maximum  stress  point  of  2076  psi,  or  about  1/2 
of  the  total  measured  stress  range.   Correlated  with  this  is 
that  the  mean  measured  stress  range  was  approximately  1750  psi 
while  the  predicted  mean  stress  range  was  820  psi,  or  about 
1/2  the  measured  mean.   The  heavier  trucks  did  not  appear 
in  the  simulated  load  data  until  approximately  1-1/2  hours 
of  simulation  had  occurred. 

Hence,  the  predicted  results,  while  small  in  guantity  appear 
to  relate  well  to  the  collected  field  data  for  lighter  truck 
loadings . 

1-495  and  U.S.  1 

Auto  traffic  was  assumed  to  be  an  arbitrary  constant  and  set 
at  83%  of  the  total  traffic  for  this  bridge.  The  truck  dis- 
tribution data  used  was  based  upon  actual  measurements  taken 
at  the  bridge  (3).  This  data  only  allowed  for  five  truck  types, 
i.e.,  2D,  3,  2S1,  2S2  and  3S2.  A  comparison  of  the  measured 
incidence  of  these  truck  types  and  the  five  hours  of  simulation 
generated  by  BRIGLD1  is  shown  in  Table  26 . 

In  this  case  no  means  of  comparing  total  or  average  truck  weights 
directly  was  possible  from  available  data. 
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TABLE  26.   Comparison  of  Truck  Distributions  for  the 
I-495/U.S.  1  Bridge 


Truck  Types 

2D 

3 

2S1 

2S2 

3S2 

Measured  (%) 

37.7 

19.5 

8.8 

15.3 

18.7 

Simultated  (%) 

37.9 

20.4 

8.7 

13.6 

19.4 
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The  manner  in  which  the  truck  traffic  appeared  on  the  bridge 
deck  in  the  simulation  generated  102  load  events  of  101 
single  truck  events  and  one  two  truck  platoon,  in  the  five 
hours  of  simulated  loading. 

This  case  also  under  generated  truck  traffic  in  comparison  to 
the  available  real  sample.   The  real  sample  acquired  637 
trucks  in  a  14.5  hours  of  sampling.   However,  the  simulated 
sample  only  produced  103  trucks  in  a  five  hour  period,  or 
only  about  1/2  the  rate  indicated  by  the  real  data.   Two 
factors  must  be  considered  in  this  instance,  i.e., 

1.  The  low  frequency  data  was  not  available  from  the  real 
sample  to  generate  the  headway  data  for  the  simulation  which 
was  generated  from  2  and  4  hour  means. 

2.  The  14-1/2  hours  of  real  sampling  was  heavily  skewed 
to  the  high  intensity  period  of  the  day,  i.e.,  from  12:00 
midnight  to  2:30  p.m.,  and  does  not  include  the  low  intensity 
periods  during  a  24  hour  cycle  nor  the  weekly  cycle. 

The  input  data  used  to  generate  this  loading  data  is  shown  in 
Figure  9  and  the  five  hour  generated  data  is  shown  in  Figure  9. 
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The  results  from  the  stress  program  did  not  appear  as 
favorable  in  this  case,  when  compared  to  field  collected  data 
per  Table  27,  as  did  the  other  simple  span  cases.   Two 
significiant  factors  enter  into  the  evaluation  of  the  com- 
parison of  real  stress  range  data  and  synthesized  data  which 
can  cause  variances  between  these  data,  i.e., 

1.  The  content  of  the  simulated  truck  traffic,  the 
weights  of  the  trucks,  their  speeds  and  the  size  of  the  sample 
as  compared  to  the  real  truck  traffic  that  occurred  in  the 
generation  of  the  real  stress  range  data. 

2.  The  theoretical  specification  of  the  structural 
characteristics  of  the  bridge  as  compared  to  its  actual 
structural  characteristics. 

While  the  stress  range  data,  as  shown  in  Table  27,  does  not 
distribute  over  the  range  of  stresses  collected  in  the  field 
it  is  contained  within  the  ranges  collected  in  the  field. 
The  synthetic  results  are  heavily  skewed  toward  the  lower 
stress  ranges.   This  tendency  could  be  caused  by  one  of  two, 
or  both,  factors  i.e., 

1.  The  generated  truck  traffic  from  BRIGLDl  was 
composed  of  lighter  trucks  on  the  whole  than  the  real  sample 
contained. 

2.  The  specified  structural  properties  in  the  input  data 
to  BRGSTRS  was  weighted  toward  a  stiffer  response  than  was 
true  of  the  real  bridge. 

Partial  output  of  this  case  is  shown  in  Figure  10. 
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1-83  and  Bunker  Hill  Road 

Auto  traffic  was  assumed  to  be  an  arbitrary  constant  and  set 
at  83%  of  the  total  traffic  for  this  bridge.   The  truck  dis- 
tribution data  used  was  based  upon  actual  measurements  taken  at 
the  bridge  (3).   This  data  only  allowed  for  five  truck  types, 
i.e.,  2D,  3,  2S1,  2S2  and  3S2.   A  comparison  of  the  measured 
incidence  of  these  truck  types  and  the  five  hours  of  simulation 
generated  by  BRIGLD1  is  shown  in  Table  28. 

In  this  case  no  means  of  comparing  total  or  average  truck 
weights  directly  was  possible  from  available  data. 

The  manner  in  which  the  truck  traffic  appeared  on  the  bridge 
deck  in  the  simulation  generated  117  load  events  of  117  single 
truck  events  in  the  five  hours  of  simulated  loading. 

The  real  sample  contained  7444  trucks  in  a  7  day  period  or  one 
truck  every  1.35  minutes.   The  simulated  sample  contained  only 
117  trucks  in  a  five  hour  period  or,  again,  approximately  1/2 
the  rate  indicated  by  the  real  data.   This  can  easily  be 
accounted  for  and  corrected  in  the  BRIGLD1  input  data.   The 
headway  tables  require  adjustment  due  to  the  lack  of  the  high 
frequency  data  masked  in  the  4  hour  means  available  from  the 
real  sample  data. 

The  input  data  used  to  generate  this  loading  data  is  shown  in 
Figure  11,  and  the  five  hour  generated  data  is  also  shown  in 
Figure  11. 
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TABLE  28.   Comparison  of  Truck  Distributions  for  the  1-83 
and  Bunker  Hill  Road  Bridge 


Truck  Types 

2D 

3 

2S1 

2S2 

3S2 

Measured  (%) 

18.4 

2.6 

7.2 

27.4 

44.4 

Simulated  (%) 

16.2 

4.3 

10.2 

28.2 

41.0 
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The  input  data  for  the  stress  calculations  is  shown  in  Figure  12. 
The  results  produced  by  BRGSTRS  as  compared  to  measured  field 
data  are  also  shown  in  Figure  12.   This  case,  of  all  the  simple 
span  cases.,  required  more  effort  to  get  results  than  all  of  the 
other  simple  span  cases  together.   Input  errors  caused  serious 
problems  and  several  reruns.   Fortunately,  this  case  only 
utilized  eight  minutes  of  CPU  time  for  the  five  hour  simu- 
lated load  period,  for  a  compression  ratio   of  37  to  1  of  real 
time.   This  bridge  was  a  three  lane  bridge  with  8  beams  and 
sampled  with  strain  gages  on  beams  4  and  5.   Initially,  the 
sample  point  input  data  to  BRGSTRS  was  established  as  the  same  as 
the  field  test.   However,  this  generated  no  output  from 
BRGSTRS.   The  reason  for  this  was  the  manner  in  which  the 
synthetic  load  data  had  been  generated  by  BRIGLD1.   There  was, 
in  this  synthetic  sample  of  truck  load  data,  a  concentration  of 
trucks  in  the  first  lane  with  little  truck  traffic  in  the 
second  lane.   This  caused  low  stresses  to  appear,  via  the 
coupling  effect,  in  beams  4  and  5  of  the  synthesized  test  case. 
Consequently,  no  stress  output  was  generated  by  BRGSTRS.   All 
of  the  stresses  were  calculated  and  the  program  proceeded 
correctly  but  the  selected  sample  points,  the  steady  state 
threshold  value  input  and  the  change  of  stress  required  by  the 
control  on  determining  stress  points  of  maxima  prevented  the 
output  of  data.   A  sample  point  at  midspan  of  beam  2 ..was  also 
specified  in  the  final  run.   It  was  selected  because  of  its 
location  under  lane  1.   The  cause  of  the  problem  was  antici- 
pated but  not  known  for  certain.   The  midspan  point  was  chosen 
for  the  obvious  reason  of  having  the  best  chance  of  acquiring 
significant  stress  values.   The  results  obtained  were  those 
expected  and  are  shown  in  comparison  with  the  field  collected 
on  and  off  the  cover  plate  on  beams  4  and  5  in  Table  29. 
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While  the  synthetic  sample  point  on  beam  2  was  on  the  cover 
plate,  it  was  also  subjected  to  higher  loadings,  due  to  the 
tendency  of  truck  traffic  to  travel  in  the  right  most  lane, 
than  the  real  data  reflected  for  beams  4  and  5  on  the  cover 
plate.   Further,  the  real  sample  points  on  beams  4  and  5  were 
at  the  ends  of  the  cover  plate  and  not  at  midspan.   A  more 
valid  comparison  perhaps  would  have  been  to  set  the  synthetic 
sample  points  as  they  were  for  beams  4  and  5,  at  the  ends  of 
the  cover  plate  on  beam  2.   Additional  variation  possibly 
occurs  because  of  weight  per  truck  variances  between  the 
synthetic  sample  generated  and  the  real  set  of  trucks  causing 
the  collected  set  of  stress  range  data. 

In  general,  the  results  of  the  synthesis  measured  at  mid  span 
of  bean  2  are  essentially  bracketed,  as  seen  in  Table  29,  by 
the  collected  data  on  and  off  the  cover   plate  on  beams  4  and 
5  and  compare  favorably. 

Dumfries  Bridge 

Auto  traffic  density  was  considered  as  an  arbitrary  constant 
and  set  at  80%  of  the  total  traffic,  leaving  the  truck  traffic 
at  20%  of  the  total.   The  truck  distribution  was  based  upon 
actual  measurements  taken  on  the  bridge  (5).   This  data  only 
allowed  for  4  truck  types,  i.e.,  2,  2S1,  2S2,  and  3S2.   A 
comparison  of  the  measured  incidence  of  these  truck  types  and 
the  five  hours  of  simulation  generated  by  BRIGLD1  is  shown  in 
Table  30. 
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TABLE  30.   Comparison  of  Truck  Distributions  for  the  Dumfries 
Bridge 


Truck  Types 

2 

2S1 

2S2 

3S2 

Measured  (%) 

19.9 

14.0 

29.7 

36.4 

Simulated  (%) 

23.35 

14.2 

26.9 

35.5 
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In  this  case,  data  was  available  to  allow  a  comparison  of  the 
truck  weights  which  caused  the  stress  ranges  of  the  real 
sample  with  the  truck  weights  generated  by  BRIGLDl  for  use 
in  predicting  the  synthetic  stress  ranges.   Table  31  presents 
a  comparison  of  these  truck  weights. 

The  real  data  also  provided  a  measure  of  truck  weights  on 
1119  trucks  over  this  bridge  for  a  total  weight  of  38,306,200 
pounds.   The  average  truck  weight  for  this  sample  was  34,232 
pounds  per  truck.   The  average  simulated  truck  weight  was 
42,538  pounds  per  truck.   The  loading  of  the  synthetic  bridge 
was,  then,  approximately  8372  pounds  per  truck  heavier  than 
this  particular  real  sample  indicated  for  the  actual  bridge. 

The  manner  in  which  the  truck  traffic  appeared  on  the  bridge 
deck  in  the  simulation  generated  196  load  events  of  195 
single  truck  events  and  one  two  truck  platoon,  in  the  five 
hours  of  simulated  loading. 

The  annual  weighing  data  over  a  two  year  period  indicates 
an  hourly  mean  truck  rate  of  67.18.   The  simulated  data  only 
provided  for  39.4  trucks  per  hour,  or  approxiamtely  only  60% 
of  the  real  rate.   The  headway  data  utilized  for  this  case 
was  estimated  from  the  given  gross  annual  rates.   However, 
correction  of  the  headway  tables  to  adjust  for  this  variance 
can  be  easily  accomplished.   Similarly,  the  variance  in  weight 
per  truck  can  be  adjusted  in  the  weight  distribution  input 
tables.   However,  the  weight  variance  would  merely  tend  to- 
ward the  conservative.   The  rate  correction  would  be  far  more 
important,  in  terms  of  the  objective  of  these  computer  programs. 

The  input  data  used  to  generate  this  loading  is  shown  in 

Figure  13,  and  the  five  hour  generated  data  is  shown  in  Figure  13 
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TABLE  31.   Comparison  of  Measured  and  Predicted  Truck  Weights 
(1-95  Dumfries,  Virginia) 


Weight  Range 


0 
10000 


1000 
20000 


0200003 


0001 
3000i 


04 


0000 
0000 


40001 
5000i 


05 


06 


0000 
0000 


60000 


70000 


70000!  80000 


Measured  (%) 
(859  Trucks) 


7.0 


15.0 


26.2 


14.9 


13.2  I  9.0    8.3 


5.4 


Weight  Range 


0 
8000 


18000 
6000 


16000J32000  4000CJ48000   64000   72000 
3200040000  14800064000  I  72000   88000 


Predicted  (%) 
(197  Trucks) 


0.0 


1.0 


29.9   21.3 


4.0 


23.9 


14.7 


1.5 
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The  gross  results  of  the  stress  maxima  generated  by  BRGSTRS, 
from  the  five  hour  simulated  traffic  sample,  are  shown  in 
Table  32.   The  total  incidence  of  maxima  over  .41  to  ap- 
proximately 1.2  KIPS  is  82.8%  for  the  real  sample,  which 
compares  to  78.4%  for  the  simulated  sample.   The  distribution 
within  this  interval  to  the  .4  to  .8  and  .8  to  1.2  intervals 
varies  by  10%  in  both  intervals. 

The  results  shown  in  Table  32  indicate  a  lower  stress  maxima 
in  the  simulated  data  than  in  the  real  data.   This  is  con- 
trary to  the  previously  mentioned  tendency  of  the  simulated 
trucks  being  an  average  of  8000  pounds  per  truck  heavier  than 
that  shown  in  the  real  data.   The  values  used  for  the  effective 
beam  moments  of  inertia  in  BRGSTRS  appear  to  have  been  about 
10%  to  high  which  would  cause  a  stiffer  bridge  and  reflect 
lower  stress  values  for  the  same  or  heavier  loads. 

At  the  midspan  point  on  beams  2  and  3,  there  was  a  total  of 
436  incidents  for  beam  2  and  102  for  beam  3.   The  ranges 
on  these  two  beams,  as  compared  to  the  field  measured  data 
are  shown  in  Table  33.   The  stress  ranges  with  the  greatest 
intensity,  on  beam  2,  i.e.,  between  .41  and  1.65,  vary 
approximately  from  2%  to  4%.   The  lower  end,  0.  to  .4,  does 
not  correlate  well  and  relates  to  the  cut  off  threshold 
stress  value  used,  and  possibly  too  stiff  a  bridge  used  in 
the  simulation.   This  incidence  variance  roughly  equals  the 
variance  on  the  higher  end  of  the  stress  ranges  between  the 
real  and  predicted  occurrences,  which  further  evidences  a 
stiffer  bridge  in  the  simulation.   Beam  3  did  not  correlate  well 
between  the  simulated  and  reat  data.   This  could  be  due  to 
either  or  both  of  two  factors,  i.e., 

1.   An  insufficient  traffic  sample  and  traffic  rate  in  the 
2nd  lane  in  the  simulation  than  were  encountered  in  the  real 
sample. 
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2.      Too  stiff  a  coupling' effect  in  the  lateral  direction 
used  in  BRGSTRS . 

In  evaluating  these  results,  it  should  be  kept  in  mind  that  the 
real  sample  represented  34,448  stress  range  occurrences,  while 
the  simulated  sample  only  represented  700  stress  range  occur- 
rences.  In  order  to  simulate  an  equivalent  amount  of 
incidents,  it  would  require  approximately  one  half  hour  of 
computer  time  to  simulate  the  traffic  and  thirty  hours  of 
computer  time  to  calculate  the  stresses. 

One  additional  note  on  this  case  is  that  it  utilized  36 
minutes  of  CPU  time  in  the  "go"  step  for  a  compression  ratio 
of  8.3  to  1  of  real  time. 

Figure  14  presents  samples  of  the  output  from  BRGSTRS  on  this 
case. 
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Overall  Results 

The  results  of  the  analysis  of  the  four  simple  span  cases 
indicates  that  both  BRIGLD1  and  BRGSTRS  can  be  used  to 
predict  stress  range  occurrences  with  at  least  ball  park 
accuracy.   However,  the  results  of  the  four  initial  attempts 
to  predict  stress  range  occurrences  can  not  be  used  as  the 
sole  measure  of  the  validity  or  accuracy  of  these  programs. 
As  was  indicated  earlier,  the  cases  presented  could  be 
easily  improved  by  especially  rerunning  BRIGLD1  with  better 
headway  data  to  improve  the  generated  truck  rate.   Further, 
the  input  data  to  BRGSTRS,  representing  the  structural 
characteristics  of  each  bridge,  could  be  further  refined  to 
improve  the  correlation  of  each  bridge's  structural 
characteristics.   No  attempt  was  made  to  improve  the  input 
data,  which  should  be  the  normal  procedure  when  performing 
design  analysis. 

One  significant  output  of  this  set  of  cases  was  some  insight 
into  the  computer  running  time  of  the  dynamic  stress  analysis 
program  BRGSTRS.   The  original  test  case  used  during  the  de- 
velopment of  the  program  utilized  truck  load  data  generated 
by  BRIGLD1,  which  in  turn,  used  the  erroneous  U.S.  301  (Md.) 
data  that  caused  traffic  densities  to  be  exceedingly  high  and 
unrealistic.   This  test  case  implied  a  1  to  1  relationship  of 
computer  time  to  real  time.   However,  the  cases  run  in  this 
effort  varied,  for  the  simple  span  cases,  from  6  to  1  up  to 
37  to  1.   The  utilization  of  computer  time  by  BRGSTRS  is  a 
direct  function  of  truck  traffic  density. 

Another  significant  output  of  the  work  performed  on  these  test 
cases  was  the  determination  of  the  extremely  slow  soltuion  of 
stresses  for  continuous  beam  cases  by  BRGSTRS  and  the  need  to 
improve  this  situation.   Such  improvement  should  be  made  before 
practical  and  frequent  use  of  the  program  is  made  on  con- 
tinuous beam  bridges. 
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A  limitation  exists  on  the  use  of  this  approach  due  to 
available  disk  storage  on  the  computer  being  used.   The 
number  of  loads  generated  by  SYNGEN  is  only  limited  by  the 
disk  storage  available  for  use  in  HISGEN  to  store  each 
single  axle  stress  time  curve,  generated  by  BRGSTRS.   Each 
curve  is  allowed  to  have  up  to  100  time  points.   A  maximum 
of  30  sample  points,  3  per  each  of  a  maximum  of  10  beams, 
is  established  in  BRGSTRS.   It  is  anticipated  that  a  nominal 
maximum  case  would  exist  for 

1.  Eleven  axle  weights 

2.  Six  speeds 

3.  Two  lanes 

4.  Thirty  sample  points 

This  requires  approximately  4000  records  of  100  words  or  400 
bytes  of  disc  storage  availability  for  use  by  HISGEN  to 
utilize  single  axle  stress  traces  generated  by  BRGSTRS  from 
the  single  axle  loads  generated  by  SYNGEN. 

Similarly,  the  number  of  truck  stress  time  traces  generated 
by  the  subroutine  TRKTRC  in  HISGEN  is  limited  by  the  available 
disc  storage  on  which  to  store  the  truck  traces.   It  is  an- 
ticipated in  this  instance  that  a  nominal  maximum  case  would 
exist  for 

1.  Six  truck  types 

2.  Ten  speeds 

3 .  Two  lanes 

4.  Thirty  sample  points 

In  this  case  each  truck  stress  trace  is  allowed  to  have  up  to  200 

time  points.   This  requires  approximately  3600  records  of  200 

words  each  or  800  bytes  each.   However,  the  program,  HISGEN, 

does  allow  up  to  twenty  different  truck  types  to  be  specified. 

This  would  increase  the  required  disk  storage  to  12,000  records 

of  2  00  words  each. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  and  recommendations  presented  here  will  be  of 
four  forms ;  i.e., 

1.  General  Comments 

2.  Identified  deficiencies  in  the  existing  programs  which 
require  resolution  and  the  recommended  resolution. 

3.  Identified  limitations  in  the  existing  programs  which 
require  extension  or  generalization  and  the  recommended  action. 

4.  Major  considerations  in  the  research  and  development 
of  the  techniques  for  attacking  the  long-term  fatique  problem 
and  recommendations  of  actions . 

5.  Other  avenues  that  should  be  investigated. 

However,  the  discussion  of  the  above  categories  will  first  be 
within  the  context  of  each  program  followed  by  a  general 
discussion.   The  programs  are: 

1.  The  BRIGLDl  Simulator 

2.  The  Structural  Analysis  Program  (BRGSTRS) 

3.  The  Stress  Signature  System,  which  is  made  up  of 

(1)  The  Synthetic  Load  Generator  (SYNGEN) , 

(2)  The  Structural  Analysis  Program  (BRGSTRS),  and 

(3)  The  Histogram  Program  (HISGEN) . 
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General  Comments 

There  is  an  inherent  flexibility  in  the  programs  which  compose 
the  two  systems  defined  in  Figures  1  and  2.   The  keystone  program 
is  BRGSTRS.   It  can  interface  with 

1.  Arbitrary  bridge  load  sources,  provided  the  data  is 
in  a  compatible  format,  on  the  input  side, 

2.  Arbitrary  analytical  programs  on  the  output  side, 
or 

3.  Provide  a  limited  analytical  output  of  its  own  as  the 
ending  program. 

This  capability  lends  itself  to  the  research  and  investi- 
gation of  other  forms  for  generating  bridge  loads  and 
predicting  long  term  stress  ranges. 

In  addition  to  theoretical  use,  data  from  other  sources, 
configured  in  the  proper  format  such  as  field  data,  may  be 
input  to  BRGSTRS  and  HISGEN,  e.g., 

1.  In  the  case  of  BRGSTRS  real  truck  data  collected 
on-site,  weight,  speed  and  lane,  from  a  controlled  experi- 
ment or  actual  data  can  be  substituted  for  the  synthetic  data 
generated  by  BRIGLD1  or  SYNGEN  on  magnetic  tape. 

2.  In  the  case  of  HISGEN  real  stress  data,  reduced  from 
collected  strain  gage  data,  correlated  with  truck  data,  can 
be  used  as  input,  in  the  proper  format,  to  HISGEN  as  a  sub- 
stitute for  the  output  of  BRGSTRS. 

In  using  BRIGLD1  and  BRGSTRS,  care  should  be  taken  in  the 
development  of  all  input  data.   In  the  case  of  BRIGLD1,  truck 
statistics  should  be  carefully  generated  and  a  simulation 
performed  and  analyzed  before  processing  the  load  data  on 
BRGSTRS.   If  the  results,  generated  truck  statistics,  do  not 
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agree  well  enough  with  known  data,  then,  the  input  data  should 
be  modified  and  the  case  rerun  on  BRIGLD1,  prior  to  using 
BRGSTRS.   BRIGLD1  is  substantially  more  economical  to  utilize 
than  BRGSTRS  and  will  be  the  basic  source  of  variance  between 
real  and  simulated  stress  data.   Providing  sufficiently  accurate 
truck  load  data  from  BRIGLD1  will  minimize  reruns  on  the 
long  running  BRGSTRS  program. 

In  the  event  the  direct  input  data  for  BRIGLD1  is  not  available, 
estimates  can  be  used  initially  and  varied  until  the  desired 
truck  statistics  are  generated  by  the  program,  at  which  time 
BRGSTRS  may  be  used. 

The  BRIGLD1  simulator 

The  identified  deficiencies  in  the  BRIGLDl  simulator  are: 

1.  Inadequate  weight-horsepower  influence  in  the 
simulation. 

2.  Lack  of  influence  of  downgrades  in  the  simulation. 
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3.   Lack  of  influence  of  the  passing  restriction  in  the 
simulation. 

'4.   In  bidirectional  traffic  flow,  crashes  occur. 

5.  Simulation  algorithm  converges  to  zero  or  specified 
minimum  speed. 

6.  Degradation  of  the  preservation  of  prescribed  traf- 
fic statistics  over  long  simulation  time  periods. 

In  the  case  of  the  weight-horsepower  effect ,  a  better  repre- 
sentation could  be  included  with  a  small  effort. 

The  downgrade  and  passing  restriction  does  not  appear  to  work 
in  the  simulation.   This  would  require  thorough  checking  of  the 
coding  and  modification  of  the  validity  of  this  conclusion. 
Correcting  it  will  be  of  a  minor  nature. 

The  bidirectional  flow  simulation  appears  to  have  problems. 
Inadvertent  "crashes"  occur  which  are  due  to  an  inadequate  re- 
presentation of  two-way  traffic.   This  requires  detailed  ana- 
lysis of  the  coding  to  determine  the  cause  and  a  small  effort 
to  correct  the  problem. 

Far  more  significant  is  the  tendency  of  the  algorithm  to  con- 
verge on  zero  speed.   It  should  oscillate  about  the  steady 
state  speed.   Logical  or  mathematical  errors,  or  an  inadequate 
representation,  apparently  are  the  cause  of  this  problem.   This 
would  require  a  moderate  effort  to  investigate.   However,  it 
might  not  be  correctable  except  by  a  complete  replacement  of 
the  motion  simulation  algorithm. 

Another  extremely  significant  apparent  deficiency  is  the  simu- 
lator's tendency  to  fail  in  the  preservation  of  the  prescribed 
traffic  statistics  over  long  simulation  periods.   The  implica- 
tion, as  indicated  earlier,  is  that  the  pseudo-random  number 
generator  is  being  used  improperly  for  the  necessary  functions. 
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Correction  of  this  is  straightforward  and  of  a  minor  nature. 
However,  there  is  an  uncertainty  that  such  a  correction  will 
eliminate  the  problem  and  adequate  testing  after  such  a  cor- 
rection would  be  necessary. 

The  identified  limitations  in  the  BRIGLDl  Simulator  are: 

1.  Limited  to  two  lanes  of  traffic,  and 

2.  Limited  to  initial  generating  of  traffic  in  the 
right  lane  only. 

In  order  to  provide  more  than  two  lanes  of  traffic  simulation, 
a  moderate  effort  is  required.   Substantial  rework  of  logic, 
data  handling  and  buffering,  and  motion  synthesis  would  be 
required.   It  is  recommended  that  if  such  an  extension  is  at- 
tempted, a  complete  replacement  of  the  traffic  simulation 
algorithm  be  made,  which  would  also  correct  all  deficiencies 
and  limitations. 

The  restriction  to  one-lane  generation  is  nonsensical  and 
should  be  replaced.   A  different  basis  of  generation  would 
be  required  but  the  effort  required  would  be  small  in  nature. 

Perhaps  the  most  significant  conclusion  is  the  basic  ineptness 
ponderousness  and  economic  infeasibility  of  the  approach  for 
handling  the  long-term  fatigue  problem  upon  which  the  BRIGLDl 
simulator  is  based.   Further,  the  non-rigorousness  of  the  approach 
places  a  large  uncertainty  upon  the  value  of  any  fatigue 
generated  by  this  technique  because  the  load  content,  trucks, 
is  uncontrollable.   This  is  regardless  of  the  validity  and 
accuracy  of  the  empirical  data  fed  into  the  simulator.   Also, 
the  above  defined  remaining  deficiencies  in  the  simulator,  even 
after  the  significant  effort  expended  in  correcting  and  im- 
proving it,  cause  a  great  deal  of  doubt  about  the  value  of  its 
generated  load  data. 

In  order  to  analyze  one  bridge  for  a  100-year  life  span,  let's 
assume  a  2-week  simulation  period,  which  will  be  representative 
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of  each  5-year  interval  in  the  100  years.   This  amounts  to 
6,720  hours  of  real  time.   At  a  compression  of  720  to  one,' 
for  the  simulator  alone,  this  requires  9.3hours  of  computer 
time.   Each  two-week  sample  could  be  multiplied  by  130,  i.e., 
for  each  stress  maxima  occurrence,  to  represent  a  five-year 
period.   The  best  estimate  at  present  on  the  stress  analysis 
calculational  time  indicates  at  best  about  a  37  to  1  compres- 
sion on  time.   This  implies  the  need  for  a  total  of  about  234 
hours  of  computer  time,  which  is  not  a  practical  approach. 

The  only  solution  to  this  problem  is  to  use  a  very  simplified 
stress  analysis  calculation,  which  will  compress  real  time  on 
the  order  of  the  simulator.   This  would  still  require  approxi- 
mately 20  hours  of  computer  time. 

It  is  concluded  that  the  concept  of  simulating  traffic  to 
generate  time  dependent  load  data  and  the  subsequent  genera- 
tion of  dynamic  stress,  for  the  prediction  of  stress  maxima, 
in  turn  for  predicting  long-term  fatigue  is  not  practical  and 
has  a  large  uncertainty  regarding  the  proper  representation  of 
truck  platoon  events/  their  configuration  and  most  importantly 
their  weight. 
The  Dynamic  Structural  Analysis  Program 

This  program  is  of  a  utilitarian  nature  to  the  overall  problem 
of  predicting  long-term  stress  maxima  histograms  for  fatigue 
analysis.   It  was  not  developed  in  the  sense  of  advancing  the 
state-of-the-art . 

If  the  traffic  simulation  approach  is  utilized  to  generate 
load  events,  an  extremely  simplified  method  will  need  to  be 
implemented  as  a  replacement  for  this  program. 

In  general,  it  is  felt  that  the  program  developed  is  adequate 
for  the  purpose  it  is  to  serve.   However,  it  requires  improvements 
from  a  utility  point  of  view,  but  these  should  be  established 
after  some  significant  use  of  the  program  has  been  made. 
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While  substantial  testing  of  this  program  was  performed 
and  several  production  type  runs  were  made  with  it  which 
provided  very  good  comparisons  with  field  collected  data, 
it  is  felt  that  substantial  further  testing,  or  very  care- 
ful production  use,  is  still  required  to  determine  its 
limitations  and  deficiencies. 

A  serious  deficiency  exists  in  BRGSTR  in  the  calculation  of 
stresses  for  multiple  span  continuous  beams.   Because  of 
the  restriction  on  core  memory,  it  was  not  feasible  to  cal- 
culate all  of  the  flexibility  matrices  and  their  modified 
inverses  for  each  span  and  save  them  for  use  during  the  entire 
run.   Consequently  they  are  currently  being  calculated  for 
each  time  point  which  causes  excessive  use  of  computer  time. 
So  much  so  that  it  is  not  recommended  that  the  program  be 
used  for  more  than  a  few  load  events,  at  present,  on  con- 
tinuous beam  cases.   This  problem,  i.e.,  use  of  large  amounts 
of  computer  time  on  continuous  beam  problems,  can  be 
minimized  by  initially  generating  these  matrices,  as  is  done 
for  the  simple  span  case,  saving  them  in  temporary  disk  files 
during  the  run  and  retrieving  them  from  disk,  at  each  time 
point  for  use. 

At  present,  the  following  suggestions  are  made  concerning  this 
problem: 

1.  Analysis  of  the  effects  of  varying  At  on  a  highly 
responsive  bridge  should  be  performed. 

2.  Comparison  of  dynamic  stresses  for  specific  loads  on 
a  given  bridge  should  be  made  against  real  data  to  provide  a 
basis  for  measuring  its  accuracy  in  predicting  stresses.   Com- 
parative testing  performed  in  this  study  was  limited  to  stress 
range  incidence  of  random  truck  traffic. 
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The  Stress  Signature  Program  System 

The  analytical  technique  developed  in  this  investigation,  based 
upon  the  idea  of  single  axle  stress  signatures  for  any  given 
bridge,  is  certainly  a  more  efficient  and  more  rigorous  ap- 
proach toward  predicting  long-term  fatigue  of  bridges  than 
the  traffic  simulation  approach.   However,  insufficient 
testing  and  evaluation  of  this  approach  does  not  allow  any 
valid  conclusions  to  be  drawn. 

Continued  improvement  and  extension  of  the  model  and  the  com- 
puter program  is  also  required  before  a  viable  tool  will  be 
forthcoming.   Computer  time  requirements  for  this  approach  are 
quite  small  when  compared  to  the  simulation  approach. 

Particular  limitations  of  the  HISGEN  program  relate  to  the 
truck  weight  distributions,  by  type,  platoon  distributions  and 
platoon  configuration   definitions  and  distributions.   The 
primary  deficiency  lies  in  the  lack  of  data  to  provide  realistic 
choices  of  these  values.   The  truck  weight  distribution  approach 
was  to  use  a  maximum  for  each  type  which  would  always  provide 
conservative  stress  range  estimates.   This  is  a  practical  ap- 
proach but  does  not  tend  toward  the  concept  of  optimum  bridge 
design.   Platoon  distribution  data  is  literally  nonexistent 
and  this  creates  a  definite  problem  in  the  use  of  HISGEN.   Two 
techniques  can  be  used  to  develop  basic  data,  i.e., 

1.  The  collection  and  analysis  of  actual  field  data,  or 

2.  The  use  of  a  good  validated  traffic  simulation  program 
and  analysis  of  its  output. 
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3.  Extension  of  the  discontinuity  moment  contribution 
due  to  a  live  load  on  one  span  to  generating  the  resultant" 
dynamic  stresses  on  the  adjacent  unloaded  spans  of  a  con- 
tinuous beam  type  bridge  should  be  made. 

4.  Other  forms  of  beams,  not  presently  synthesized 
in  BRGSTRS ,  should  be  included.   This  can  be  accomplished 
by  replacing  or  providing  an  alternative  option  switch  to 
the  subroutine  MODULI.   It  is  this  subroutine  which 
characterizes  a  beam  or  diaphram. 

5.  Extending  the  number  of  lanes  appears  to  be  essential. 
Of  the  comparative  tests  run  using  real  data,  only  one  was  a 
bridge  with  two  lanes.   All  others  were  three  lanes. 

6 .  Providing  for  user  or  automatic  specification  of 
the  scaling  of  the  stress  graphical  output. 

7.  Modifying  the  program  to  calculate  the  flexibility 
matrix  and  the  inverse  of  its  modified  form  for  each  span  of 
a  continuous  beam  type  bridge  only  once,  i.e.,  in  the 
utilization  block,  storing  them  in  temporary  disc  files, 

and  retrieving  each  matrix  as  necessary  at  each  time  point 
instead  of  calculating  them  at  each  time  point. 

In  general,  a  great  deal  of  learning  must  be  accomplished  about 

using  this  program  before  significant 'criticism  can  be  made. 

One  significant  feature  of  the  program,  not  a  deficiency, 

is  its  ability  to  essentially  allow  a  user  to  outsmart  himself. 

Careful  selection  of  output  options  and  sample  points  is 

mandatory.   Careless  specification  can  result  in  no  stress  output 

even  though  it  was  calculated  for  all  elements.   Conversely, 

a  user  can  be  overwhelmed  with  output  by  careless  specification 

of  output  options,  especially  the  debug  options. 
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Probably  the  single  weakest  point  in  HISGEN,  and  again  its 
due  to  the  lack  of  data  primarily  but  also  due  to  limiting 
complexity,  are  the  definitions  of  the  platoon  configurations 
used  in  the  program  and  the  establishing  of  their  distribution 
by  a  user  in  the  input  data.   To  provide  improvements  on  these 
two  factors,  the  same  two  approaches  specified  above  for  the 
platoon  distribution  problem,  are  again  the  practical  means. 

This  program,  and  SYNGEN,  like  the  others  previously  discussed 
should  be  extended  to  more  than  2  lanes.   The  total  time 
nominally  necessary  is  then  approximately  .75  hours.   This  is 
in  comparison  to  the  234  hours  of  computer  time  necessary  to 
the  simulation  approach  previously  discussed. 

Merely  from  the  pragmatic  consideration  of  computer  time,  this 
approach  warrants  serious  further  investigation  and  development. 
Controlled  field  experiments  and  the  collection  of  validation 
data  andr  parametric  data  is  necessary  in  order  to  evaluate  this 
method.   Also,  from  the  rigorous  point  of  view,  it  has  a  lot  of 
merit.   It  certainly  does  not  have  the  extent  of  uncertainties 
that  the  BRIGLDk  simulator  contains,  or  the  general  simulation 
approach. 

This  approach  also  allows  the  practical  use  of  a  fairly  realistic 
dynamic  structural  analysis  program.   Whereas,  the  simulator 
approach  must  be  constrained  to  use  an  extremely  simplified 
stress  analysis  approach  to  significantly  reduce  computer 
requirements . 

Additionally,  continued  research  to  develop  such  approaches 
is  required  and  should  be  continued.   It  is  felt  that  the  tech- 
nique developed  is  of  an  innovative  or  novel  nature  and  could 
lead  to  improved  insight  into  better  approaches,  as  compared 
to  the  classical  engineering  and  simulation  approaches. 
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Nominally,  about  twenty  minutes  of  computer  time  will  be  more 
than' adequate  to  generate  the  single  axle  stress  signatures 
using  BRGSTRS .   Twenty  minutes  allows  the  generation  of  axles 
for  eleven  weights  for  each  of  six  speeds  in  each  of  two  lanes 
for  approximately  24  sample  points,  i.e.,  3168  stress  curves. 
This  provides  a  relatively  fine  grid  which,  in  turn,  should 
provide  rather  good  interpolative  results.   The  histogram 
program  will  nominally  require  about  40  minutes  for  any  length 
of  bridge  life  span. 

Other  Considerations 

An  alternative  approach  to  the  simulation  of  traffic  down  a 
synthetic  highway  to  a  synthetic  bridge  undergoing  analysis 
is  to  stochastically  generate  the  load  events  at  the  entrance 
to  the  deck  in  all  lanes  simultaneously.   A  load  event  would  be 
defined  by  the  existence  of  a  truck  in  the  set  of  vehicles  en- 
tering the  bridge.   An  analytical  approach  can  be  established 
for  implementing  such  a  model.   Further,  a  conceptual  variation 
of  the  manner  in  which  the  dynamic  stress  analysis  program  is 
applied,  somewhat  analogous  to  the  stress  signature  concept, 
could  reduce  the  computer  time  to  a  few  hours ,  as  shown  pre- 
viously for  the  stress  signature  method.   This  overall  approach 
retains  a  rigorous  stress  analysis  capability,  generates  re- 
alistic load  events  without  simulating  traffic  and  would  use  a 
practical  amount  of  computer  time  for  long  life  span  periods. 

Summary 

In  general,  it  is  recommended  that  continued  research  or 
analytical  methods  such  as  the  Stress  Signature  Method,  in- 
cluding further  development  of  it,  and  stochastic  load  event 
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generation  techniques  which  allow  the  use  of  fairly  rigorous- 
dynamic  stress  analysis  programs,  and  utilize  practical 
amounts  of  computer  time  be  performed.   It  is  also  concluded 
that  the  use  of  highway  traffic  simulation  techniques  are 
economically  impractical  for  long-term  fatigue  analysis  and 
that  they  should  be  abandoned. 
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